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PREFACE. 



This book was planned as the third book in the series of three textbooks for 
three basic geometric disciphnes of the university education. These are 

- «Coursc of analytical gconictry^»; 

- «Course of linear algebra and multidimensional geometry*; 

- «Course of differential geometry*. 

This book is devoted to the first acquaintance with the differential geometry. 

Therefore it begins with the theory of curves in three-dimensional Euclidean space 
E. Then the vectorial analysis in E is stated both in Cartesian and curvilinear 
coordinates, afterward the theory of surfaces in the space E is considered. 

The newly fashionable approach starting with the concept of a differentiable 
manifold, to my opinion, is not suitable for the introduction to the subject. In 
this way too many efforts are spent for to assimilate this rather abstract notion 
and the rather special methods associated with it, while the the essential content 
of the subject is postponed for a later time. I think it is more important to make 
faster acquaintance with other elements of modern geometry such as the vectorial 
and tensorial analysis, covariant differentiation, and the theory of Riemannian 
curvature. The restriction of the dimension to the cases n = 2 and n = 3 is 
not an essential obstacle for this purpose. The further passage from surfaces to 
higher-dimensional manifolds becomes more natural and simple. 

I am grateful to D. N. Karbushev, R. R. Bakhitov, S. Yu. Ubiyko, D. I. Borisov, 
and Yu. N. Polyakov for reading and correcting the manuscript of the Russian 
edition of this book. 

November, 1996; 

December, 2004. R. A. Sharipov. 



^ Russian versions of the second and the third boolfs were written in 1096, but the first book 
is not yet written. I understand it as my duty to complete the series, but I had not enough time 
all these years since 1996. 



CHAPTER I 



CURVES IN THREE-DIMENSIONAL SPACE. 



1. Curves. Methods of defining a curve. 
Regular and singular points of a curve. 



Let E be a three-dimensional Euclidean point space. The strict mathematical 
definition of such a space can be found in [1]. However, knowing this definition 
is not so urgent. The matter is that E can be understood as the regular 
three-dimensional space (that in which we live). The properties of the space E 
are studied in elementary mathematics and in analytical geometry on the base 
intuitively clear visual forms. The concept of a line or a curve is also related to 
some visual form. A curve in the space E is a spatially extended one-dimcnsional 
geometric form. The one-dimensionality of a curve reveals when we use the 
vectorial-parametric method of defining it: 



We have one degree of freedom when choosing a point on the curve (1.1), our 
choice is determined by the value of the numeric parameter t taken from some 
interval, e. g. from the unit interval [0, 1] on the real axis M. Points of the curve 
(1.1) are given by their radius-vectors^ r = r{t) whose components x^{t), x'^{t), 
x^{t) are functions of the parameter t. 

The continuity of the curve (1.1) means that the functions x^{t), x'^{t), x^{t) 
should be continuous. However, this condition is too weak. Among continuous 
curves there are some instances which do not agree with our intuitive understand- 
ing of a curve. In the course of mathematical analysis the Peano curve is often 
considered as an example (see [2]). This is a continuous parametric curve on a 
plane such that it is enclosed within a unit square, has no self intersections, and 
passes through each point of this square. In order to avoid such unusual curves 
the functions x'^{t) in (1.1) are assumed to be continuously differentiable (C^ class) 
functions or, at least, piecewise continuously differentiable functions. 

Now let's consider another method of defining a curve. An arbitrary point of 
the space E is given by three arbitrary parameters — its coordinates. 

We can restrict the degree of arbitrariness by considering a set of points whose 
coordinates cc^, a;^, x^ satisfy an equation of the form 



r = r(t) 



x\t) 
x'{t) 
x^t) 



(LI) 



F{x\x^,x^) = 0, 



(1.2) 



Here we assume that some Cartesian coordinate system in E is taken. 
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where F is some continuously dilFerentiable lunction of three variables. In a 
typical situation formula (1.2) still admits two-parametric arbitrariness: choosing 
arbitrarily two coordinates of a point, we can determine its third coordinate by 
solving the equation (1.2). Therefore, (1.2) is an equation of a surface. In 
the intersection of two surfaces usually a curve arises. Hence, a system of two 
equations of the form (1.2) defines a curve in E: 



G{x\x^,x^) = 0. 



(1.3) 



If a curve lies on a plane, we say that it is a plane curve. For a plane 
curve one of the equations (1.3) can be replaced by the equation of a plane: 
Ax^ + Bx^ + Cx^ + D = 0. 

Suppose that a curve is given by the equations (1.3). Let's choose one of the 
variables a;^, x^, or x^ for a parameter, e.g. we can take a;^ = t to make certain. 
Then, writing the system of the equations (1.3) as 

F{t,x'^,x^) = 0, 
G{t,x'^,x^) = 0, 

and solving them with respect to x^ and x'^, we get two functions x'^{t) and x^{t). 
Hence, the same curve can be given in vectorial-parametric form: 



r(t) 



t 

x'{t) 
x^t) 



Conversely, assume that a curve is initially given in vectorial-parametric form 
by means of vector-functions (1.1). Then, using the functions x^{t), x^{t), x^{t), 
we construct the following two systems of equations: 



(x^-x^{t) = 0, (x^~x^{t) = 0, 

[x^-x^{t) = 0, \x^-x^t)^Q. ^^'^^ 

Excluding the parameter t from the first system of equations (1.4), we obtain some 
functional relation for two variable x^ and x^. We can write it as F{x'^,x'^) = 0. 
Similarly, the second system reduces to the equation G{x^,x^) = 0. Both these 
equations constitute a system, which is a special instance of (1.3): 

r F{x\x^) = 0, 
I G{x\x^) =0. 

This means that the vectorial-parametric representation of a curve can be trans- 
formed to the form of a system of equations (1.3). 

None of the above two methods of defining a curve in E is absolutely preferable. 
In some cases the first method is better, in other cases the second one is used. 
However, for constructing the theory of curves the vectorial-parametric method 
is more suitable. Suppose that we have a parametric curve 7 of the smoothness 
class C^. This is a curve with the coordinate functions x^(t), x'^{t), x^{t) being 
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continuously differentiable. Let's choose two different values of the parameter: 
t and t = t + At, where At is an increment of the parameter. Let A and B 
be two points on the curve corresponding to that two values of the parameter 
t. We draw the straight line passing through these points A and B; this is a 

secant for the curve 7. Directing vectors 
of th is secant are coUinear to the vector 
AB . We choose one of them: 




a = 



AB _ r(t + A<)-r(i) 



At 



At 



(1-5) 



Tending At to zero, we find that the 
point B moves toward the point A. Then 
the secant tends to its limit position and 
becomes the tangent line of the curve at 
the point A. Therefore limit value of the 
vector (1.5) is a tangent vector of the 
curve 7 at the point A: 



rlt) = lim a 

At^oo 



dr(t) 
dt 



= r(t). (1.6) 



Fig. 1.1 



The components of the tangent vector 
(1.6) are evaluated by differentiating the 
components of the radius-vector r(t) with respect to the variable t. 

The tangent vector r(t) determines the direction of the instantaneous displace- 
ment of the point r(t) for the given value of the parameter t. Those points, at 
which the derivative r(t) vanishes, are special ones. They are «stopping points*. 
Upon stopping, the point can begin moving in quite different direction. For 
example, let's consider the following two plane curves: 



r(t) 



r(t) 



(1.7) 



At t = both curves (1.7) pass through the origin and the tangent vectors of both 
curves at the origin are equal to zero. However, the behavior of these curves near 
the origin is quite different: the first curve has a beak-like fracture at the origin. 




Fig. 1.2 



Fig. 1.3 
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while the second one is smooth. Therefore, vanishing of the derivative 

T{t) - Ht) = 



(1.8) 



is only the necessary, but not sufficient condition for a parametric curve to have a 
singularity at the point r{t). The opposite condition 



T{t) = rit) ^ 



(1.9) 



guaranties that the point r{t) is free of singularities. Therefore, those points of a 
parametric curve, where the condition (1.9) is fulfilled, are called regular points. 

Let's study the problem of separating regular and singular points on a curve 
given by a system of equations (1.3). Let A = {a^,a^,a^) be a point of such a 



curve. The functions F{x^,x'^,x^) and G{x^, 
continuously differcntiablc. The matrix 



dF 


dF 


dF 




dx^ 


dx^ 


dG 


dG 


dG 


dx^ 


dx^ 


dx^ 



,x ) in (1.3) are assumed to be 



(1.10) 



composed of partial derivatives of F and G at the point A is called the Jacobi 
matrix or the Jacobian of the system of equations (1.3). If the minor 



Ml = det 



dF 


dF 


dx"^ 


d^ 


dG 


dG 


dx"^ 


dx^ 



^0 



in Jacobi matrix is nonzero, the equations (1.3) can be resolved with respect to 
x'^ and x^ in some neighborhood of the point A. Then we have three functions 
x^ = x^ = x^(t), x^ = x'^it) which determine the parametric representation of 
our curve. This fact follows from the theorem on implicit functions (see [2]). Note 
that the tangent vector of the curve in this parametrization 



1 



is nonzero because of its first component. This means that the condition il/i ^ 
is sufficient for the point ^ to be a regular point of a curve given by the system of 
equations (1.3). Remember that the Jacobi matrix (1.10) has two other minors: 



A'h = det 



dF 


dF 


dx^ 


dx^ 


dG 


dG 


dx^ 


dx^ 



Ah = det 



dF 


dF 


dx^ 


dx"^ 


dG 


dG 


dx^ 


dx"^ 
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For both of them the similar propositions are fulfilled. Therefore, we can formulate 
the following theorem. 

Theorem 1.1. A curve given by a system of equations (1.3) is regular at all 
points, where the rank of its Jacobi matrix (1.10) is equal to 2. 

A plane curve lying on the plane = can be defined by one equation 
F{x^,x'^) = 0. The second equation here reduces to = 0. Therefore, 
G{x^,x'^,x^) = x^. The Jacoby matrix for the system (1.3) in this case is 



J = 



dF OF ^ 
dx^ dx'^ 
1 



(1.11) 



If rank J ~ 2, this means that at least one of two partial derivatives in the matrix 
(1.11) is nonzero. These derivatives form the gradient vector for the function F: 



OF 

dx"^ 



Theorem 1.2. A plane curve given by an equation F{x^, x'^) 
all points where gradF ^ 0. 



is regular at 



This theorem 1.2 is a simple corollary from the theorem 1.1 and the relationship 
(1.11). Note that the theorems 1.1 and 1.2 yield only sufficient conditions for 
regularity of curve points. Therefore, some points where these theorems are not 
applicable can also be regular points of a curve. 



§ 2. The length integral 
and the natural parametrization of a curve. 

Let r = r{t) be a parametric curve of smoothness class C^, where the parameter 
t runs over the interval [a, b]. Let's consider a monotonic increasing continuously 
differentiable function ip(t ) on a segment [a, b] such that ip{a) = a and (p{b) = b. 
Then it takes each value from the segment [a, b] exactly once. Substituting 
t = (fit) into r(i), we define the new vector-function f(t) = r((^(i)), it describes 
the same curve as the original vector- function v{t). This procedure is called the 
reparametrization of a curve. We can calculate the tangent vector in the new 
parametrization by means of the chain rule: 

T{i) = ip'{i)-r{^{i)). (2.1) 

Here ^p'it) is the derivative of the function </j(t). The formula (2.1) is known 
as the transformation rule for the tangent vector of a curve under a change of 
parametrization. 

A monotonic decreasing function (p{t) can also be used for the reparametrization 
of curves. In this case (p{a) = b and 99(6) = a, i. e. the beginning point and the 
ending point of a curve are exchanged. Such reparametrizations are called changing 
the orientation of a curve. 

From the formula (2.1), we see that the tangent vector f{i) can vanish at some 
points of the curve due to the derivative 'fi'{t) even when T{if{t)) is nonzero. 
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Certainly, such points are not actually the singular points of a curve. In order 
to exclude such formal singularities, only those reparametrizations of a curve 
are admitted for which the function (pit) is a strictly monotonic function, i.e. 
ip'{i) > or ip'{i) < 0. 

The formula (2.1) means that the tangent vector of a curve at its regular point 
depends not only on the geometry of the curve, but also on its parametrization. 
However, the effect of parametrization is not so big, it can yield a numeric factor to 
the vector t only. Therefore, the natural question arises: is there some preferable 
parametrization on a curve ? The answer to this question is given by the length 
integral. 

Let's consider a segment of a parametric curve of the smoothness class with 
the parameter t running over the segment [a, 6] of real numbers. Let 

a = to <ti < ... <tn=b (2.2) 

be a series of points breaking this segment into n parts. The points r(io), • ■ • , r(in) 

on the curve define a polygonal line with 
n segments. Denote Atk = tk — tk~i and 
let £ be the maximum of At^: 




Fig. 2.1 



£ — max Atfc. 

A:— 1, ... , n 

The quantity £ is the fineness of the par- 
tition (2.2). The length of fc-th segment 
of the polygonal line AB is calculated by 
the formula Lk = \^{tk) — ^{tk~i)\- Us- 
ing the continuous differentiability of the 
vector- function r{t), from the Taylor ex- 
pansion of r{t) at the point tk-i we get 
Lk = \'''{tk-i)\ ■ Atk + o{e). Therefore, as 
the fineness e of the partition (2.2) tends 
to zero, the length of the polygonal line 
AB has the limit equal to the integral of 
the modulus of tangent vector T{t) along 
the curve: 



L = lim y Lk = / |x(t)| dt. (2.3) 

a 

It is natural to take the quantity L in (2.3) for the length of the curve AB. Note 
that if we reparametrize a curve according to the formula (2.1), this leads to a 
change of variable in the integral. Nevertheless, the value of the integral L remains 
unchanged. Hence, the length of a curve is its geometric invariant which does not 
depend on the way how it is parameterized. 

The length integral (2.3) defines the preferable way for parameterizing a curve 
in the Euclidean space E. Let's denote by s{t) an antiderivative of the function 
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ipit) = \T{t)\ being under integration in the formula (2.3): 



t 



sit)^ \Tit)\dt. (2.4) 



to 

Definition 2.1. The quantity s determined by the integral (2.4) is called the 
natural parameter of a curve in the Euclidean space E. 

Note that once the reference point r(to) and some direction (orientation) on a 
curve have been chosen, the value of natural parameter depends on the point of 
the curve only. Then the change of s for — s means the change of orientation of 
the curve for the opposite one. 

Let's differentiate the integral (2.4) with respect to its upper limit t. As a result 
we obtain the following relationship: 

§ = k(t)l- (2.5) 

Now, using the formula (2.5), we can calculate the tangent vector of a curve in its 
natural parametrization, i. c. when s is used instead of t as a parameter: 

dr dr dt dr / ds T 
ds ~ (ft ' ds ~ dt/ " |t|' 

From the formula (2.6), we see that in the tangent vector of a curve in natural 
parametrization is a unit vector at all regular points. In singular points this vector 
is not defined at all. 

§ 3. Prenet frame. The dynamics of Prenet 
frame. Curvature and torsion of a spacial curve. 

Let's consider a smooth parametric curve r(s) in natural parametrization. The 
components of the radius-vector r(s) for such a curve are smooth functions of 
s (smoothness class C°°). They are differentiable unlimitedly many times with 
respect to s. The unit vector t(s) is obtained as the derivative of r(s): 

xW^f. (3^1) 

Let's differentiate the vector t(s) with respect to s and then apply the following 
lemma to its derivative t'(s). 

Lemma 3.1. The derivative of a vector of a constant length is a vector perpen- 
dicular to the original one. 

Proof. In order to prove the lemma we choose some standard rectangular 
Cartesian coordinate system in E. Then 

|r(s)p = (r(s) I t(s)) = {t^)^ + (r^)^ + (r^)^ = const . 
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Let's differentiate this expression with respect to s. As a result we get the 
following relationship: 

= 2 ri (ri)' + 2 (r^)' + 2 r'' (r^)' = 0. 

One can easily see that this relationship is equivalent to (t(s) | t'(s)) — 0. Hence, 
x(s) _L 'r'(s). The lemma is proved. □ 

Due to the above lemma the vector r'(s) is perpendicular to the unit vector 
t{s). If the length of t'(s) is nonzero, one can represent it as 

r'(s) = fc(s) • n(s), (3.2) 

where k{s) = \t'(s)\ and |n(s)| — 1. The scalar quantity fc(s) = [''"'(•s)! in formula 
(3.2) is called the curvature of a curve, while the unit vector n(s) is called its 
primary normal vector or simply the normal vector of a curve at the point r(s). 
The unit vectors t(s) and n(s) arc orthogonal to each other. We can complement 
them by the third unit vector b(s) so that t, n, b become a right triple^: 

b(s) = [r(s), n(s)]. (3.3) 

The vector b(s) defined by the formula (3.3) is called the secondary normal 
vector or the binormal vector of a curve. Vectors t(s), n(s), b(s) compose an 
orthonormal right basis attached to the point r(s). 

Bases, which are attached to some points, are usually called frames. One should 
distinguish frames from coordinate systems. Cartesian coordinate systems are also 
defined by choosing some point (an origin) and some basis. However, coordinate 
systems are used for describing the points of the space through their coordinates. 
The purpose of frames is different. They are used for to expand the vectors which, 
by their nature, are attached to the same points as the vectors of the frame. 

The isolated frames are rarely considered, frames usually arise within families 
of frames: typically at each point of some set (a curve, a surface, or even the whole 
space) there arises some frame attached to this point. The frame r(s), n(s), b(s) 
is an example of such frame. It is called the Frenet frame of a curve. This is the 
moving frame: in typical situation the vectors of this frame change when we move 
the attachment point along the curve. 

Let's consider the derivative n'(s). This vector attached to the point r(s) can 
be expanded in the Frenet frame at that point. Due to the lemma 3.1 the vector 
n'(s) is orthogonal to the vector n(s). Therefore its expansion has the form 

n'(s) ==a-r(s) + >fb(s). (3.4) 

The quantity a in formula (3.4) can be expressed through the curvature of the 



^ A non-coplanar ordered triple of vectors ai , a2 , aa is called a right triple if, upon moving 
these vectors to a common origin, when looking from the end of the third vector as, we see the 
shortest rotation from ai to a2 as a counterclockwise rotation. 



14 



CHAPTER I. CURVES IN THREE-DIMENSIONAL SPACE. 



curve. Indeed, as a result of the following calculations we derive 

a(.) = (r(.)|n'(.)) = (r(.)|n(.))'- 

- {t'{s) I n(.)) ^ -(kis) ■ nis) \ n{s)) = -k{s). 

The quantity x = >c{s) cannot be expressed through the curvature. This is an 
additional parameter characterizing a curve in the space E. It is called the torsion 
of the curve at the point r — r(s). The above expansion (3.4) of the vector n'(s) 
now is written in the following form: 

n'(.s) = -k{s) ■ t{s) + x{s) ■ b(s). (3.6) 

Let's consider the derivative of the binormal vector b'(s). It is perpendicular 
to b(s). This derivative can also be expanded in the Frenet frame. Due to 
b'(s) _L b(s) we have b'(s) = /? • n(s) + 7 • r(s). The coefficients (3 and 7 in this 
expansion can be found by means of the calculations similar to (3.5): 

I3{s) = (n(.) I b'(.)) ^ (n(.) I b(s))' - (n'(s) | b(.)) = 
= -(-fc(.s) • t{s) + x{s) ■ b(s) I b(.s)) = ^h{s). 

7(s) = (r(s) I b'(.)) = {t{s) I b(,s))' - (r'(s) | h{s)) = 
= -(fc(s)-n(s)|b(s)) = 0. 

Hence, for the expansion of the vector b'(s) in the Frenet frame we get 

h'{s) ^ ->i{s) - nis). (3.7) 

Let's gather the equations (3.2), (3.6), and (3.7) into a system: 

t'{s) = k{s) ■ n{s), 

n'(s) = -k{s) ■ r(s) + >f(s) • b(s), (3.8) 
b'(s) = ->f(s) • n(s). 

The equations (3.8) relate the vectors t(s), n(s), b(s) and their derivatives with 
respect to s. These differential equations describe the dynamics of the Frenet 
frame. They are called the Frenet equations. The equations (3.8) should be 
complemented with the equation (3.1) which describes the dynamics of the point 
r(s) (the point to which the vectors of the Frenet frame are attached). 

§ 4. The curvature center and the curvature radius 
of a spacial curve. The evolute and the evolvent of a curve. 

In the case of a planar curve the vectors t{s) and n(s) lie in the same plane as 
the curve itself. Therefore, binormal vector (3.3) in this case coincides with the 
unit normal vector of the plane. Its derivative b'(s) is equal to zero. Hence, due 
to the third Frenet equation (3.7) we find that for a planar curve >c{s) = 0. The 
Frenet equations (3.8) then are reduced to 

T'(s) = fc(s)-n(s), 

n'{s) = -k{s)-r{s). ^ ■ ' 
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Let's consider the circle of the radius R with the center at the origin lying in the 
coordinate plane — 0. It is convenient to define this circle as follows: 



r(s) 



R cos{s/R) 
R sm{s/R) 



(4.2) 



here s is the natural parameter. Substituting (4.2) into (3.1) and then into (3.2), 
we find the unit tangent vector t(s) and the primary normal vector n(s): 



r(,5) 



- sin(s/i?) 
cos(s/i?) 



- cos(s/i?) 

- sin(s/i?) 



(4.3) 



Now, substituting (4.3) into the formula (4.1), we calculate the curvature of a 
circle k{s) = 1/R = const. The curvature fc of a circle is constant, the inverse 
curvature 1/fc coincides with its radius. 

Let's make a step from the point r(s) on a circle to the distance l//c in the 
direction of its primary normal vector n(s). It is easy to see that we come to the 
center of a circle. Let's make the same step for an arbitrary spacial curve. As a 
result of this step we come from the initial point r(s) on the curve to the point 
with the following radius-vector: 



pis) = i-(s) + 



n(s) 



(4.4) 



Certainly, this can be done only for that points of a curve, where fc(s) 7^ 0. The 
analogy with a circle induces the following terminology: the quantity R{s) = 
l/k{s) is called the curvature radius, the point with the radius- vector (4.4) is 
called the curvature center of a curve at the point r(s). 

In the case of an arbitrary curve its curvature center is not a fixed point. When 
parameter s is varied, the curvature center of the curve moves in the space drawing 
another curve, which is called the evolute of the original curve. The formula (4.4) 
is a vectorial-parametric equation of the evolute. However, note that the natural 
parameter s of the original curve is not a natural parameter for its evolute. 

Suppose that some spacial curve r(t) is given. A curve f (s) whose evolute p{s) 
coincides with the curve v{t) is called an evolvent of the curve v{t). The problem 
of constructing the evolute of a given curve is solved by the formula (4.4). The 
inverse problem of constructing an evolvent for a given curve appears to be more 
complicated. It is effectively solved only in the case of a planar curve. 

Let r(s) be a vector-function defining some planar curve in natural parametri- 
zation and let f(s) be the evolvent in its own natural parametrization. Two 
natural parameters s and s are related to each other by some function in form 
of the relationship s = (p{s). Let "0 = be the inverse function for then 
s = ?/'(s). Using the formula (4.4), now we obtain 



r(V'(s)) = f(s) 



n(£) 



(4.5) 



Let's differentiate the relationship (4.5) with respect to s and then let's apply the 
formula (3.1) and the Frenct equations written in form of (4.1): 



^'(s) • r(^(S)) = - 



d~s \k{s) 



n{S). 
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Here t('0(s)) and n(s) both are unit vectors which are coUinear due to the above 
relationship. Hence, we have the foUowing two equaUties: 



fi(s) = ±t(V'(s)), = ±: 

The second equahty (4.6) can be integrated: 

^ ±W:s)-c). 



1 



d§ Vfc(s) 



k{S) 



(4.6) 



(4.7) 



Here C is a constant of integration. Let's combine (4.7) with the first relationship 
(4.6) and substitute it into the formula (4.5): 

f (S) = r(^.(S)) + {C~ ■ rim)- 

Then wc substitute s ~ V'(^) iiito the above formula and denote p(s) = r((p(s)). 
As a result we obtain the following equality: 



pis) = ri.s) + {C-s)-T{s). 



(4.8) 



The formula (4.8) is a parametric equation for the evolvent of a planar curve r(s). 
The entry of an arbitrary constant in the equation (4.8) means the evolvent is not 
unique. Each curve has the family of evolvents. This fact is valid for non-planar 
curves either. However, we should emphasize that the formula (4.8) cannot be 
applied to general spacial curves. 

§ 5. Curves as trajectories of material points in mechanics. 

The presentation of classical mechanics traditionally begins with considering 
the motion of material points. Saying material point, we understand any material 
object whose sizes are much smaller than its displacement in the space. The 
position of such an object can be characterized by its radius-vector in some 
Cartesian coordinate system, while its motion is described by a vector-function 
r{t). The curve r{t) is called the trajectory of a material point. Unlike to purely 
geometric curves, the trajectories of material points possess preferable parameter t, 
which is usually distinct from the natural parameter s. This preferable parameter 
is the time variable t. 

The tangent vector of a trajectory, when computed in the time parametrization, 
is called the velocity of a material point: 



v(t) = ^ = rit) 



v\t) 
v\t) 



(5.1) 



The time derivative of the velocity vector is called the acceleration vector. 



- Tt 



v(i) 



a\t) 
a\t) 
a\t) 



(5.2) 
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The motion of a material point in mechanics is described by Newton's second law: 

ma = F(r,v). (5.3) 

Here m is the mass of a material point. This is a constant characterizing the 
amount of matter enclosed in this material object. The vector F is the force 
vector. By means of the force vector in mechanics one describes the action of 
ambient objects (which are sometimes very far apart) upon the material point 
under consideration. The magnitude of this action usually depends on the position 
of a point relative to the ambient objects, but sometimes it can also depend on 
the velocity of the point itself. Newton's second law in form of (5.3) shows that 
the external action immediately affects the acceleration of a material point, but 
neither the velocity nor the coordinates of a point. 

Let s = s[t) be the natural parameter on the trajectory of a material point 
expressed through the time variable. Then the formula (2.5) yields 

s{t)^W{t)\^v{t). (5.4) 

Through v(t) in (5.4) we denote the modulus of the velocity vector. 

Let's consider a trajectory of a material point in natural parametrization: 
r = r(s). Then for the velocity vector (5.1) and for the acceleration vector (5.2) 
we get the following expressions: 

v(t)=,s(t)-r(,s(i)). 

Taking into account the formula (5.4) and the first Frcnet equation, these expres- 
sions can be rewritten as 

v(t) = «(t) • r(s(t)), 

a(0 = v{t) ■ Tis{t)) + (kisit)) vitf) ■ n(.(0). ^ ' ' 

The second formula (5.5) determines the expansion of the acceleration vector into 
two components. The first component is tangent to the trajectory, it is called the 
tangential acceleration. The second component is perpendicular to the trajectory 
and directed toward the curvature center. It is called the centripetal acceleration. 
It is important to note that the centripetal acceleration is determined by the 
modulus of the velocity and by the geometry of the trajectory (by its curvature). 
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§ 1. Vectorial and tensorial fields in the space. 

Let again E be a three-dimensional Euclidean point space. We say that in 
E a vectorial field is given if at each point of the space E some vector attached 
to this point is given. Let's choose some Cartesian coordinate system in E; in 
general, this system is skew-angular. Then we can define the points of the space 
by their coordinates x^, , , and, simultaneously, we get the basis ei, e2, e-^ for 
expanding the vectors attached to these points. In this case we can present any 
vector field F by three numeric functions 



F = 



Fi(x) 
F2(x) 
F3(x) 



(1.1) 



where x — [x^,x'^,x'^) are the components of the radius- vector of an arbitrary 
point of the space E. Writing F(x) instead of F(a;^, x^, x^), we make all formulas 
more compact. 

The vectorial nature of the field F reveals when we replace one coordinate 
system by another. Let (1.1) be the coordinates of a vector field in some 
coordinate system O, ei, e2, e3 and let O, ei, §2, 63 be some other coordinate 
system. The transformation rule for the components of a vectorial field under a 
change of a Cartesian coordinate system is written as follows: 



i=i 



(1.2) 



Here 5^ are the components of the transition matrix relating the basis ei, e2, 
with the new basis ei, 62, 63, while a^,a^,a'^ are the components of the vector 

00 in the basis ei, e2, e3. 

The formula (1.2) combines the transformation rule for the components of a 
vector under a change of a basis and the transformation rule for the coordinates of 
a point under a change of a Cartesian coordinate system (see [1] ) . The arguments 
X and X beside the vector components F* and in (1.2) is an important novelty 
as compared to [1]. It is due to the fact that here we deal with vector fields, not 
with separate vectors. 

Not only vectors can be associated with the points of the space E. In linear 
algebra along with vectors one considers covectors, linear operators, bilinear forms 
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and quadratic forms. Associating some covector with each point of E, we get a 
covector field. If we associate some hnear operator with each point of the space, 
we get an operator field. An finally, associating a bilinear (quadratic) form with 
each point of E, we obtain a field of bilinear (quadratic) forms. Any choice of a 
Cartesian coordinate system O, ei, 82, 63 assumes the choice of a basis 61,62,63, 
while the basis defines the numeric representations for all of the above objects: 
for a covector this is the list of its components, for linear operators, bilinear 
and quadratic forms these are their matrices. Therefore defining a covector field 
F is equivalent to defining three functions -F'i(x), F2(x), -F3(x) that transform 
according to the following rule under a change of a coordinate system: 

i^.(x)=^7^f F,(x), 

(1.3) 



In the case of operator field F the transformation formula for the components of 
its matrix under a change of a coordinate system has the following form: 



3 3 



x^=i:s:^ 

p=l 



P=i,=i (^4) 



For a field of bilinear (quadratic) forms F the transformation rule for its compo- 
nents under a change of Cartesian coordinates looks like 



3 3 



3 

p=i 



p=lg=l ^^j.^ 



Each of the relationships (1.2), (1.3), (1.4), and (1.5) consists of two formulas. 
The first formula relates the components of a field, which are the functions of 
two different sets of arguments x = {x^,x^,x'^) and x = {x^,x'^,x^). The second 
formula establishes the functional dependence of these two sets of arguments. 

The first formulas in (1.2), (1.3), and (1.4) are different. However, one can see 
some regular pattern in them. The number of summation signs and the number 
of summation indices in their right hand sides are determined by the number of 
indices in the components of a field F. The total number of transition matrices 
used in the right hand sides of these formulas is also determined by the number of 
indices in the components of F. Thus, each upper index of F implies the usage of 
the transition matrix S, while each lower index of F means that the inverse matrix 
T = 5'-Ms used. 
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The number of indices of the field F in the above examples doesn't exceed 
two. However, the regular pattern detected in the transformation rules for the 
components of F can be generalized for the case of an arbitrary number of indices: 

^E^^l---^^---^ (1-6) 

91 •••9s 

The formula (1.6) comprises the multiple summation with respect to (r + s) indices 
pi, ■ ■ ■ ,Pr and qi, ■ ■ ■ ,qs each of which runs from 1 to 3. 

Definition 1.1. A tensor of the type (r, s) is a geometric object F whose 
components in each basis are enumerated by (r + s) indices and obey the transfor- 
mation rule (1.6) under a change of basis. 

Lower indices in the components of a tensor are called covariant indices, upper 
indices are called contravariant indices respectively. Generalizing the concept of 
a vector field, we can attach some tensor of the type (r, s), to each point of the 
space. As a result we get the concept of a tensor field. This concept is convenient 
because it describes in the unified way any vectorial and covectorial fields, operator 
fields, and arbitrary fields of bilinear (quadratic) forms. Vectorial fields are fields 
of the type (1,0), covectorial fields have the type (0, 1), operator fields are of the 
type (1, 1), and finally, any field of bilinear (quadratic) forms are of the type (0, 2). 
Tensor fields of some other types are also meaningful. In Chapter IV we consider 
the curvature field with four indices. 

Passing from separate tensors to tensor fields, we acquire the arguments in for- 
mula (1.6). Now this formula should be written as the couple of two relationships 
similar to (1.2), (1.3), (1.4), or (1.5): 



Pl...Pr 



(1.7) 



The formula (1.7) expresses the transformation rule for the components of a 
tensorial field of the type (r, s) under a change of Cartesian coordinates. 

The most simple type of tensorial fields is the type (0,0). Such fields are 
called scalar fields. Their components have no indices at all, i. e. they are numeric 
fimctions in the space E. 

§ 2. Tensor product and contraction. 

Let's consider two covectorial fields a and b. In some Cartesian coordinate 
system they are given by their components ai(x) and &j (x). These are two sets of 
functions with three functions in each set. Let's form a new set of nine functions 
by multiplying the functions of initial sets: 

Cy(x) = a,(x)&j(x). (2.1) 

Applying the formula (1.3) we can express the right hand side of (2.1) through the 
components of the fields a and b in the other coordinate system: 

c.,(x)= (^Tfa,) (e^'O =EE^f^/(«^M- 

\p=l / \9=1 / p=lq=l 
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If we denote by Cpg(x) the product of ai(x) and 6j(x), then we find that the quan- 
tities Cij(x) and Cpg(x) are related by the formula (1-5). This means that taking 
two covectorial fields one can compose a field of bilinear forms by multiplying the 
components of these two covectorial fields in an arbitrary Cartesian coordinate 
system. This operation is called the tensor product of the fields a and b. Its result 
is denoted as c = a (g) b. 

The above trick of multiplying components can be applied to an arbitrary pair 
of tensor fields. Suppose we have a tensorial field A of the type (r, s) and another 
tensorial field B of the type (m, n). Denote 

- (2-2) 

Definition 2.1. The tensor field C of the type (r+7Ti, s+n) whose components 
are determined by the formula (2.2) is called the tensor product of the fields A 
and B. It is denoted C = A (g) B. 

This definition should be checked for correctness. We should make sure that 
the components of the field C are transformed according to the rule (1.7) when we 
pass from one Cartesian coordinate system to another. The transformation rule 
(1.7), when applied to the fields A and B, yields 

A.il...ir _ \ ^ qil qir rpqi rpqs lpi...pr 



Pr + r, 

qs + 1... qs+„ 



The summation in right hand sides of this formulas is carried out with respect 
to each double index which enters the formula twice — once as an upper index 
and once as a lower index. Multiplying these two formulas, we get exactly the 
transformation rule (1.7) for the components of C. 

Theorem 2.1. The operation of tensor product is associative, this means that 
(A(gB)(gC = A®(B(gC). 

Proof. Let A be a tensor of the type (r, s), let B be a tensor of the type 
{m,n), and let C be a tensor of the type (j>,q). Then one can write the following 
obvious numeric equality for their components: 



I All. ..It rylr+l---lr + m \ ^ Ir + m+ 1 • • • Jr+ 7n+p 



Aii...ir f pir+l---ir+m ^ 1 ■ ■ ■ V+m+p A 



(2.3) 



As we see in (2.3), the associativity of the tensor product follows from the 
associativity of the multiplication of numbers. □ 

The tensor product is not commutative. One can easily construct an example 
illustrating this fact. Let's consider two covectorial fields a and b with the 
following components in some coordinate system: a = (1,0,0) and b = (0,1,0). 
Denote c = a (g) b and d = b (g a. Then for ci2 and di2 with the use of the formula 
(2.2) we derive: Ci2 = 1 and di2 = 0. Hence, c 7^ d and a (g b 7^ b g) a. 
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Let's consider an operator field F. Its components Fj{x.) are the components of 
the operator F(x) in the basis ei, 62, 63. It is known that the trace of the matrix 
F*(x) is a scalar invariant of the operator F(x) (see [1]). Therefore, the formula 

3 

/(x)=trF(x)=5^^^/(x) (2.4) 

determines a scalar field /(x) in the space E. The sum similar to (2.4) can be 
written for an arbitrary tcnsorial field F with at least one upper index and at least 
one lower index in its components: 

(x) = E f;i::;::^;:;:::-;z' (-)• (2-5) 

In the formula (2.5) the summation index k is placed to m-th upper position 
and to 71-th lower position. The succeeding indices im, ■ ■ ■ V-i and j„, . . . js-i 
in writing the components of the field F are shifted one position to the right as 
compared to their positions in left hand side of the equality (2.5): 



I I 



h 




-1 


k 






ir-l 



h 




]n-l 


k 


in 




7s -1 



T f 



Definition 2.2. The tensor field H whose components are calculated according 
to the formula (2.5) from the components of the tensor field F is called the 
contraction of the field F with respect to m-th and n-th indices. 

Like the definition 2.1, this definition should be tested for correctness. Let's 
verify that the components of the field H are transformed according to the 
formula (1.7). For this purpose we write the transformation rule (1.7) applied to 
the components of the field F in right hand side of the formula (2.5): 

api---Pr-l 
Pqi-.-Qs-i 

rriljl rpqn-1 rpl3 rpq„ rpQe-l pPl ■ ■ ■ Pm - 1 " Pm ■ • • Pr - 1 

^ -^ji ■ ■ ■ ^J^-ik ■ ■ ■ ^gi...g„_i/39„...g,_i ' 

In order to derive this formula from (1.7) we substitute the index k into the m-th 
and n-tli positions, then we shift all succeeding indices one position to the right. 
In order to have more similarity of left and right hand sides of this formula we 
shift summation indices as well. It is clear that such redesignation of summation 
indices does not change the value of the sum. 

Now in order to complete the contraction procedure we should produce the 
summation with respect to the index k. In the right hand side of the formula the 
sum over k can be calculated explicitly due to the formula 

E^arf = 5f, (2.6) 

fe=l 
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which means T — S ^. Due to (2.6) upon calculatmg the sum over k one can 
calculate the sums over P and a. Therein we take into account that 



E 



pPl...Pm_l QPm...Pr-l _ ffPl - Pr 



As a result we get the equality 

ir-l \ ^ oil qir-l 7^9l rpQa-l frpi...pr_i 

~ '^pi ■ ■ ■ '^Pr-i ^Ji ■ ■ ■ ^qi...qs-l ' 

Pl---Pr-1 

qi...qs-l 

which exactly coincides with the transformation rule (1.7) written with respect to 
components of the field H. The correctness of the definition 2.2 is proved. 

The operation of contraction introduced by the definition 2.2 implies that the 
positions of two indices are specified. One of these indices should be an upper 
index, the other index should be a lower index. The letter C is used as a 
contraction sign. The formula (2.5) then is abbreviated as follows: 

H = C„,„(F) = C(F). 

The numbers m and n are often omitted since they are usually known from the 
context. 

A tensorial field of the type (1, 1) can be contracted in the unique way. For a 
tensorial field F of the type (2, 2) we have two ways of contracting. As a result 
of these two contractions, in general, wc obtain two different tensorial fields of the 
type (1, 1). These tensorial fields can be contracted again. As a result wc obtain 
the complete contractions of the field F, they are scalar fields. A field of the type 
(2, 2) can have two complete contractions. In general case a field of the type (n, n) 
has n! complete contractions. 

The operations of tensor product and contraction often arise in a natural way 
without any special intension. For example, suppose that we are given a vector 
field V and a covector field w in the space E. This means that at each point we 
have a vector and a covector attached to this point. By calculating the scalar 
products of these vectors and covectors we get a scalar field / = (w | v) . In 
coordinate form such a scalar field is calculated by means of the formula 

3 

/ = (2.7) 

k=l 

From the formula (2.7), it is clear that / = C(w (g) v). The scalar product 
/ = (w I v) is the contraction of the tensor product of the fields w and v. In a 
similar way, if an operator field F and a vector field v arc given, then applying F 
to V we get another vector field u = F v, where 



3 



In this case we can write: u = C(F (X) v); although this writing cannot be uniquely 
interpreted. Apart from u = F v, it can mean the product of v by the trace of the 
operator field F. 
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§ 3. The algebra of tensor fields. 

Let V and w be two vectorial fields. Then at each point of the space E we have 
two vectors v(x) and w(x). We can add them. As a result we get a new vector 
field u = V + w. In a similar way one can define the addition of tensor fields. Let 
A and B be two tensor fields of the type (r, s). Let's consider the sum of their 
components in some Cartesian coordinate system: 

c;-::;:-^::X + B-^-::X. o.i) 

Definition 3.1. The tensor field C of the type (r, s) whose components are 
calculated according to the formula (3.1) is called the sum of the fields A and B 
of the type (r, s). 

One can easily check up the transformation rule (1.7) for the components of the 
field C. It is sufficient to write this rule (1.7) for the components of A and B then 
add these two formulas. Therefore, the definition 3.1 is consistent. 

The sum of tensor fields is commutative and associative. This fact follows 
from the commutativity and associativity of the addition of numbers due to the 
following obvious relationships: 

, r>ii...i,. _ Tjii...i,. . Aii...ir 



Let's denote by T^^r.s) the set of tensor fields of the type {r,s). The tensor 
multiplication introduced by the definition 2.1 is the following binary operation: 

The operations of tensor addition and tensor multiplication (3.2) are related to 
each other by the distributivity laws: 

(A + B)(8)C = A®C + B(8)C, 

C«)(A + B) = C®A + C®B. ^ ' ' 

The distributivity laws (3.3) follow from the distributivity of the multiplication of 
numbers. Their proof is given by the following obvious formulas: 

/ Ail...ir I nil... ,^4t-+1--- ir+m _ 

_ Ail...ir ^ir+l...ir + m _|_ r>ii...ir ^ ir+1 ■ ■ ■ ir+m 



^■ll...lr /lJr+1--- «r+m .^ll...lr ri'r + 1 ■ ■ ■ «r+ m 



Due to (3.2) the set of scalar fields K = 2^(o,o) (which is simply the set 
of numeric functions) is closed with respect to tensor multiplication (g), which 
coincides here with the regular multiplication of numeric functions. The set K is 
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a commutative ring (see [3]) with the unity. The constant function equal to 1 at 
each point of the space E plays the role of the unit element in this ring. 

Let's set TO = 71 = in the formula (3.2). In this case it describes the 
multiplication of tensor fields from T(^r,s) by numeric functions from the ring 
K. The tensor product of a field A and a scalar filed £ K is commutative: 
A (g) ^ = ^ (g) A. Therefore, the multiplication of tensor fields by numeric functions 
is denoted by standard sign of multiplication: ^ (g) A = ^ • A. The operation of 
addition and the operation of multiplication by scalar fields in the set T(^r,s) possess 
the following properties: 

(1) A + B = B + A; 

(2) (A + B) + C = A+(B + C); 

(3) there exists a field G T(^r,s) such that A + = A for an arbitrary tensor 
field A e T(^_^); 

(4) for any tensor field A G Tf^r.s) there exists an opposite field A' such that 
A + A' = 0; 

(5) ^ • (A + B) = ^ • A + ^ • B for any function ^ from the ring K and for any 
two fields A,B G T^^r.s)'-, 

(6) {£, + Q ■ A — £^ ■ A + ( ■ A ior any tensor field A G T/^r.s) a-nd for any two 
functions ^, C G K; 

(7) (^ C) ■ = ^ ■ (C ■ -^) fo'' a-iiy tensor field A G T(r s) and for any two functions 

(8) 1 • A = A for any field A G Ti^r,s)- 

The tensor field with identically zero components plays the role of zero element 
in the property (3). The field A' in the property (4) is defined as a field whose 
components are obtained from the components of A by changing the sign. 

The properties (l)-(8) listed above almost literally coincide with the axioms of 
a linear vector space (see [1]). The only discrepancy is that the set of functions K 
is a ring, not a numeric field as it should be in the case of a linear vector space. 
The sets defined by the axioms (l)-(8) for some ring K are called modules over 
the ring K or K -modules. Thus, each of the sets r(r.s) is a module over the ring of 
scalar functions K = T(^o,o)- 

The ring K ~ ^(o.o) comprises the subset of constant functions which is 
naturally identified with the set of real numbers M. Therefore the set of tensor 
fields T(r.s) ill the space E is a linear vector space over the field of real numbers M. 

If r ^ 1 and s ^ 1, then in the set T(^r,s) the operation of contraction with 
respect to various pairs of indices are defined. These operations are linear, i. e. the 
following relationships are fulfilled: 

C(A + B) = C(A) + C(B), 

c(^A) = ^c•(A). 

The relationships (3.4) are proved by direct calculations in coordinates. For the 
field C = A + B from (2.5) we derive 



Yjii... i, — 1 X ^ ^ 11 ... im-i fe i,„ ... v-i 



k=l 



kijn...ir—i , \ j-fil . . . ijn — i k . . 

fc=l fe=l 



26 



CHAPTER II. ELEMENTS OF TENSORIAL ANALYSIS. 



This equality proves the first relationship (3.4). In order to prove the second one 
we take C = ^ • A. Then the second relationship (3.4) is derived as a result of the 
following calculations: 



3 



jl---js-l ~ jl---jn-lkj„...j. 



3 3 
fe=l k=l 



The tensor product of two tensors from TJ-^ .j) belongs to 7(^,8) only if r = s = 
(see formula (3.2)). In all other cases one cannot perform the tensor multiplication 
staying within one A'-module T(^r.s)- Iii order to avoid this restriction the following 
direct sum is usually considered: 

CXD OO 

T = ^^T^r,sy (3.5) 

r=0 s=0 

The set (3.5) consists of finite formal sums A^^' + . . . + A'''") , where each summand 
belongs to some of the A'-modules T(^r,s)- The operation of tensor product is 
extended to the if-module T by means of the formula: 

k q 

(A(i^ + . . . + A^^)) ® (A(i^ + . . . + A^")) - ^ ^ aW (g) A(^'). 

i=l 3 = 1 

This extension of the operation of tensor product is a bilinear binary operation in 
the set T. It possesses the following additional properties: 

(9) (A + B)®C = A®C + B®C; 

(10) (e-A)®C = e-(A®C); 

(11) C(g)(A + B) = C®A + C®B; 

(12) C®(^B)^e-(C®B). 

These properties of the operation of tensor product in T are easily derived from 
(3.3). Note that a A'-module equipped with an additional bilinear binary operation 
of multiplication is called an algebra over the ring K or a K-algebra. Therefore 
the set T is called the algebra of tensor fields. 

The algebra T is a direct sum of separate iiT-modules T^^^s) in (3.5). The 
operation of multiplication is concordant with this expansion into a direct sum; 
this fact is expressed by the relationship (3.2). Such structures in algebras are 
called gradings, while algebras with gradings are called graded algebras. 

§ 4. Symmetrization and alternation. 

Let A be a tensor filed of the type (r, s) and let r ^ 2. The number of upper 
indices in the components of the field A is greater than two. Therefore, we can 
perform the permutation of some pair of them. Let's denote 

^;;:::.!r;.::;;;.}: = ^};::::r.^:::::.}:- (4.1) 

The quantities "j^ in (4.1) are produced from the components of the tensor 
field A by the transposition of the pair of upper indices im and i„. 
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Theorem 4.1. The quantities B^^ '^ J produced from the components of a ten- 
sor field A by the transposition of any pair of upper indices define another tensor 
field B of the same type as the original field A. 

Proof. In order to prove the theorem let's check up that the quantities (4.1) 
obey the transformation rule (1.7) under a change of a coordinate system: 

nil-.-V _ \ ^ oil qin, aim qir rpqi rpQs ^pi-.-Pr 

ji---js Pi ■ ■ ■ Pm ■ ■ ■ '-'pn ■ ■ ■ '-'pt ^ji ■ ■ ■ -^js ^qi---q3 ■ 

qi---q3 

Let's rename the summation indices pm and p„ in this formula: let's denote Pm 
by Pn and vice versa. As a result the S matrices will be arranged in the order 
of increasing numbers of their upper and lower indices. However, the indices Pm 
and Pn in Ap^ - ^'j will exchange their positions. It is clear that the procedure of 
renaming summation indices does not change the value of the sum: 



r>il---ir \ ^ cil Qir rpqi rpqs Ap 

^n -js - -^pi • ■ • -^p, ^ii • ■ • ^j, ^q_ 



l...Pn. 

qi--qs 



Pi...Pr 
qi---qs 



Due to the equality B^^-'-p^ = ylPi - Pn - Pm. -PT- ^j^j; above formula is exactly the 
transformation rule (1-7) written for the quantities (4.1). Hence, they define a 
tensor field B. The theorem is proved. □ 

There is a similar theorem for transpositions of lower indices. Let again A be a 
tensor field of the type (r, s) and let s ^ 2. Denote 

uii ir Ail ir (A 2) 

Theorem 4.2. The quantities B^^ ' produced from the components of a ten- 
sor field A by the transposition of any pair of lower indices define another tensor 
held B of the same type as the original field A. 

The proof of the theorem 4.2 is completely analogous to the proof of the 
theorem 4.1. Therefore we do not give it here. Note that one cannot transpose 
an upper index and a lower index. The set of quantities obtained by such a 
transposition docs not obey the transformation rule (1.7). 

Combining various pairwise transpositions of indices (4.1) and (4.2) we can 
get any transposition from the symmetric group S,- in upper indices and any 
transposition from the symmetric group 6^ in lower indices. This is a well-known 
fact from the algebra (see [3]). Thus the theorems 4.1 and 4.2 define the action of 
the groups and 6s on the A'-module T(^r,s) composed of the tensor fields of the 
type (r, s). This is the action by linear operators, i. e. 

a o t(A + B) = cr o r(A) + a = t(B), 

a.r(C-A)=e-(aor(A)) ^^'^^ 

for any two transpositions a £ 6^ and r e S^. When written in coordinate form, 
the relationship B — ct<=t(A) looks like 

where the umbers (t(1), . . . ,cr(r) and r(l), . . . , r(s) are obtained by applying a 
and T to the numbers 1, . . . , r and 1, . . . , s. 
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Definition 4.1. A tensorial field A of tlie type (r, s) is said to be symmetric 
in TO-tli and n-th upper (or lower) indices if cr(A) = A, where a is the permutation 
of the indices given by the formula (4.1) (or the formula (4.2)). 

Definition 4.2. A tensorial field A of the type (r, s) is said to be skew- 
symmetric in TO-th and n-th upper (or lower) indices if cr(A) = —A, where a is 
the permutation of the indices given by the formula (4.1) (or the formula (4.2)). 

The concepts of symmetry and skew-symmetry can be extended to the case 
of arbitrary (not necessarily pairwise) transpositions. Let e = a°T be some 
transposition of upper and lower indices from (4.4). It is natural to treat it as an 
element of direct product of two symmetric groups: e G 6,- x &s (see [3]). 

Definition 4.3. A tensorial field A of the type {r,s) is symmetric or skew- 
symmetric with respect to the transposition e G &r x &.n if one of the following 
relationships is fulfilled: e(A) = A or e(A) = (—1)" ■ A. 

If the field A is symmetric with respect to the transpositions ei and £2, then it 
is symmetric with respect to the composite transposition ei ° £2 and with respect 
to the inverse transpositions e^^ and £^^. Therefore the symmetry always takes 
place for some subgroup 25 € &r x S^. The same is true for the skew-symmetry. 

Let (5 C &r X &s be a subgroup in the direct product of symmetric groups and 
let A be a tensor field from T^^r.s)- The passage from A to the field 



is called the symmetrization of the tensor field A by the subgroup C S^. x 6^. 
Similarly, the passage from A to the field 



is called the alternation of the tensor field A by the subgroup C Sr x 6,. 

The operations of symmetrization and alternation are linear operations, this 
fact follows from (4.3). As a result of symmetrization (4.5) one gets a field B 
symmetric with respect to ©. As a result of alternation (4.6) one gets a field skew- 
symmetric with respect to 0. If © = ©r x Ss then the operation (4.5) is called 
the complete symmetrization^ while the (4.6) is called the complete alternation. 



The smoothness class of a tensor field A in the space E is determined by the 
smoothness of its components. 

Definition 5.1. A tensor field A is called an in-times continuously differen- 
tiable field or a field of the class C™ if all its components in some Cartesian system 
are m-times continuously differentiablc functions. 

Tensor fields of the class are often called differentiablc tensor fields, while 
fields of the class C°° are called smooth tensor fields. Due to the formula (1.7) the 
choice of a Cartesian coordinate system does not affect the smoothness class of a 




(4.5) 




(4.6) 



§ 5. Differentiation of tensor fields. 
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field A in tlie definition 5.1. Tlie components of a field of tlie class C™ are the 
functions of the class C™ in any Cartesian coordinate system. This fact proves 
that the definition 5.1 is consistent. 

Let's consider a differentiable tensor field of the type (r, s) and let's consider all 
of the partial derivatives of its components: 

The number of such partial derivatives (5.1) is the same in all Cartesian coordinate 
systems. This number coincides with the number of components of a tensor field 
of the type (r, s + 1). This coincidence is not accidental. 

Theorem 5.1. The partial derivatives of the components of a differentiable 
tensor field A of the type (r, s) calculated in an arbitrary Cartesian coordinate 
system according to the formula (5.1) arc the components of another tensor Bled 
B of the type {r,s + 1). 

Proof. The proof consists in checking up the transformation rule (1.7) for the 
quantities BjlZjlj^ ^ in (5.1). Let 0, 61,62,63 and O', 61,62, §3 be two Cartesian 
coordinate systems. By tradition we denote by S and T the direct and inverse 
transition matrices. Let's write the first relationship (1.7) for the field A and let's 
differentiate both sides of it with respect to the variable x^'+^ : 



Pl...Pr 

qi---q. 

In order to calculate the derivative in the right hand side we apply the chain rule 
that determines the derivatives of a composite function: 

93+1=1 

The variables x = {x^,x'^,x^) and x = (a;^,i^,i^) are related as follows: 

3 3 

One of these two relationships is included into (1.7), the second being the inversion 
of the first one. The components of the transition matrices S and T in these 
formulas arc constants, therefore, we have 

Let's substitute (5.3) into (5.2), then substitute the result into the above expression 
for the derivatives dA^^'^'^jJOx^'^'^. This yields the equality 

r>il...ir _ \ ^ oil Oir T^?! rpqs + 1 rtPl-.-Pr 

Pl...Pr 
9l---9s + l 
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which coincides exactly with the transformation rule (1-7) applied to the quantities 
(5.1). The theorem is proved. □ 

The passage from A to B in (5.1) adds one covariant index js+i- This is the 
reason why the tensor field B is called the covariant differential of the field A. 
The covariant differential is denoted as B = VA. The upside-down triangle V is 
a special symbol, it is called nabla. In writing the components of B the additional 
covariant index is written beside the nabla sign: 

Due to (5.1) the sign Vfc in the formula (5.4) replaces the differentiation operator: 
Vfc = d/dx^. However, for Vfe the special name is reserved, it is called the operator 
of covariant differentiation or the covariant derivative. Below (in Chapter III) we 
shall see that the concept of covariant derivative can be extended so that it will 
not coincide with the partial derivative any more. 

Let A be a differentiable tensor field of the type (r, s) and let X be some 
arbitrary vector field. Let's consider the tensor product VA ® X. This is the 
tensor field of the type (r + 1, s + 1). The covariant differentiation adds one 
covariant index, while the tensor multiplication add one contravariant index. We 
denote by VxA = C(VA ® X) the contraction of the field VA ® X with respect 
to these two additional indices. The field B = VxA has the same type (r, s) as 
the original field A. Upon choosing some Cartesian coordinate system we can 
write the relationship B = VxA in coordinate form: 

^};:::}:-E^^v,A-;;-. (5.5) 

q=l 

The tensor field B = VxA with components (5.5) is called the covariant derivative 
of the field A along the vector field X. 

Theorem 5.2. The operation of covariant differentiation of tensor fields posses- 
ses the following properties 

(1) Vx(A + B) = VxA + VxB; 

(2) Vx+yA = VxA + VyA; 

(3) Vf.xA = ^ • VxA; 

(4) Vx(A (g) B) = VxA (g) B + A (g) VxB; 

(5) VxC(A) = C(VxA); 

where A and B are arbitrary differentiable tensor fields, while X and Y are arbi- 
trary vector fields and ^ is an arbitrary scalar field. 

Proof. It is convenient to carry out the proof of the theorem in some Cartesian 
coordinate system. Let C = A + B. The property (1) follows from the relationship 

^ dx'i ^ dxi ^ dxi 

q=l g=l q=l 

Denote Z = X + Y and then we derive the property (2) from the relationship 

q=l 9=1 9=1 
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In order to prove the property (3) we set Z = ^ • X. Then 

q=l q=l 

This relationship is equivalent to the property (3) in the statement of the theorem. 

In order to prove the fourth property in the theorem one should carry out the 
following calculations with the components of A, B and X: 

3 \ / ^ BA^^- - \ 

9=1 \ 9=1 / 

.„+i...^„+„ ^.i...^, [Yx". ^=+1-^=+" I 



dxi 

. 9=1 



And finally, the following series of calculations 

3 ^ / 3 

o 



EVq ___ I lil--- im-l fe«m--- ir-l j _ 

Qxl \ ^ jl---j^-lkj„...]s-l I ~ 

9=1 \fc=l / 

3 3 kim-.-ir-l 



k=l 9=1 

proves the fifth property. This completes the proof of the theorem in whole. □ 

§ 6. The metric tensor and the volume pseudotensor. 

Let 0,ei,e2,e3 be some Cartesian coordinate system in the space E. The 
space E is equipped with the scalar product. Therefore, the basis 61,62,63 of any 
Cartesian coordinate system has its Gram matrix 

g,j = {ei\ej). (6.1) 

The gram matrix g is positive and non-degenerate: 

detg>0. (6.2) 

The inequality (6.2) follows from the Silvester criterion (see [1]). Under a change 
of a coordinate system the quantities (6.1) are transformed as the components 
of a tensor of the type (0,2). Therefore, we can define the tensor field g 
whose components in any Cartesian coordinate system are the constant functions 
coinciding with the components of the Gram matrix: 

5y (x) = g^j = const . 

The tensor field g with such components is called the metric tensor. The metric 
tensor is a special tensor field. One should not define it. Its existence is 
providentially built into the geometry of the space E. 

Since the Gram matrix g is non-degenerate, one can determine the inverse 
matrix g = g^^. The components of such matrix are denoted by g*-' , the indices i 
and j are written in the upper position. Then 

3 

T.9''9ik^5]. (6.3) 
j=i 
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Theorem 6.1. The components of the inverse Gram matrix g are transformed 
as the components of a tensor held of the type (2, 0) under a change of coordinates. 

Proof. Let's write the transformation rule (1.7) for the components of the 
metric tensor g: 

3 3 

In matrix form this relationship is written as 

g = T"gr. (6.4) 
Since g, g, and T are non-degenerate, we can pass to the inverse matrices: 

g-' = {T'^iTr'^Si-'S'\ (6.5) 
Now we can write (6.5) back in coordinate form. This yields 

5''=EE^;^^'5''- (6-6) 

The relationship (6.6) is exactly the transformation rule (1.7) written for the 
components of a tensor field of the type (2,0). Thus, the theorem is proved. □ 

The tensor field g = g^^ with the components g*-' is called the inverse metric 
tensor or the dual metric tensor. The existence of the inverse metric tensor also 
follows from the nature of the space E which has the pre-built scalar product. 

Both tensor fields g and g are symmetric. The symmetruy of gij with respect to 
the indices i and j follows from (6.1) and from the properties of a scalar product. 
The matrix inverse to the symmetric one is a symmetric matrix too. Therefore, 
the components of the inverse metric tensor g^^ are also symmetric with respect to 
the indices i and j. 

The components of the tensors g and g in any Cartesian coordinate system are 
constants. Therefore, we have 

Vg = 0, Vg = 0. (6.7) 

These relationships follow from the formula (5.1) for the components of the 
covariant differential in Cartesian coordinates. 

In the course of analytical geometry (see, for instance, [4]) the indexed object 
Eijk is usually considered, which is called the Levi-Civita symbol. Its nonzero 
components are determined by the parity of the transposition of indices: 

if i ~ j, i = k, or j ~ fc, 

1 if (i j/c) is even, i. e. sign{ijk) = l, (6.8) 
1 if (ijk) is odd, i.e. sign(ijA:) = —1. 

Recall that the Levi-Civita symbol (6.8) is used for calculating the vectorial prod- 
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uct^ and the mixed product^ through the coordinates of vectors in a rectangular 
Cartesian coordinate system with a right orthonormal basis ei, 62, 63: 

[X, Y]=^eJ^^e.,,X^r'= 

2 = 1 \j"=lfc = l 

3 3 3 

(X, Y, z) = Y.T.T.'^^''^'^'^'- 

i=l j=l k=l 

The usage of upper or lower indices in writing the components of the Levi- 
Civita symbol in (6.8) and (6.9) makes no difference since they do not define 
a tensor. However, there is a tensorial object associated with the Levi-Civita 
symbol. In order to construct such an object we apply the relationship which is 
usually proved in analytical geometry: 

333 

^ ^ ^ Spqi Mip Mjq Mki = det M • e^Jk (6.10) 

p=l q=l 1=1 

(see proof in [4]). Here M is some square 3x3 matrix. The matrix M can be 
the matrix of the components for some tensorial field of the type (2,0), (1, 1), or 
(0,2). However, it can be a matrix without any tensorial interpretation as well. 
The relationship (6.10) is valid for any square 3x3 matrix. 

Using the Levi-Civita symbol and the matrix of the metric tensor g in some 
Cartesian coordinate system, we construct the following quantities: 

Wijfc = Vdctg e^Jk■ (6.11) 

Then we study how the quantities ujijk and ojpqi constructed in two different 
Cartesian coordinate systems O, ei, 62, 63 and O', ei, §2, §3 are related to each 
other. From the identity (6.10) we derive 

333 

E E E ^/ ^P,i = det T (6.12) 

p=l q=l 1=1 

In order to transform further the sum (6.12) we use the relationship (6.4), as 
an immediate consequence of it we obtain the formula dot g = (det T)^ det g. 
Applying this formula to (6.12), we derive 

333 

E E E ^/ = ^^Sn{det T) v/d^ (6.13) 

p=l g=l 1 = 1 

Note that the right hand side of the relationship (6.13) differs from ujijk in (6.11) 
only by the sign of the determinant: sign(detr) = sign(detS') = ±1. Therefore, 
we can write the relationship (6.13) as 

333 

oj,,k = sign(det ^) E E E ^P,i ■ (6.14) 

p=l q=l 1=1 

^ It is also called the cross product of vectors. 

2 The mixed product is defined as (X, Y, Z) = (X | [Y, Z]). 
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Though the difference is only in sign, the relationship (6.14) differs from the 
transformation rule (1.6) for the components of a tensor of the type (0,3). The 
formula (6.14) gives the cause for modifying the transformation rule (1.6): 

F;i:t- E (-i)^^;...^Tr---Jt (e-is) 

Pl-..Pr 
<3l-- - 9a 

Here (— 1)"^ = sign(det 5) = ±1. The corresponding modification for the concept 
of a tensor is given by the following definition. 

Definition 6.1. A pseudotensor F of the type (r, s) is a geometric object 
whose components in an arbitrary basis are enumerated by (r + s) indices and 
obey the transformation rule (6.15) under a change of basis. 

Once some pseudotensor of the type (r, s) is given at each point of the space E, 
we have a pseudotcnsorial field of the type (r,s). Due to the above definition 6.1 
and due to (6.14) the quantities from (6.11) define a pseudotcnsorial field uj 
of the type (0, 3). This field is called the volume pseudotensor. Like metric tensors 
g and g, the volume pseudotensor is a special field pre-built into the space E. Its 
existence is due to the existence of the pre-built scalar product in E. 

§ 7. The properties of pseudotensors. 

Pseudotensors and pseudotcnsorial fields arc closely relative objects for tensors 
and tcnsorial fields. In this section wc repeat most of the results of previous 
sections as applied to pseudotensors. The proofs of these results are practically 
the same as in purely tensorial case. Therefore, below we do not give the proofs. 

Let A and B be two pseudotcnsorial fields of the type (r, s). Then the formula 
(3.1) determines a third field C = A + B which appears to be a pseudotcnsorial 
field of the type (r, s). It is important to note that (3.1) is not a correct procedure 
if one tries to add a tensorial field A and a pseudotcnsorial field B. The sum 
A + B of such fields can be understood only as a formal sum like in (3.5). 

The formula (2.2) for the tensor product appears to be more universal. It 
defines the product of a field A of the type (r, s) and a field B of the type (to, n). 
Therein each of the fields can be either a tensorial or a pseudotcnsorial field. The 
tensor product possesses the following properties: 

(1) the tensor product of two tensor fields is a tensor field; 

(2) the tensor product of two pseudotcnsorial fields is a tensor field; 

(3) the tensor product of a tensorial field and a pseudotcnsorial field is a 
pseudotcnsorial field. 

Let's denote by P(r.s) the set of pseudotcnsorial fields of the type {r,s). Due 
to the properties (l)-(3) and due to the distributivity relationships (3.3), which 
remain valid for pseudotcnsorial fields too, the set P{r,s) is a module over the ring 
of scaral fields K = T(o.o)- As for the properties (l)-(3), they can be expressed in 
form of the relationships 

^{r,s) ^ P(m,7i) ^ -^(r+m,.s+n) : P{r,s) ^ ^{m^n) ^ -^(r+m,s+n) 



that extend the relationship (3.2) from the section 3. 
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The formula (2.5) defines the operation of contraction for a field F of the type 
(r, s), where r ^ 1 and s ^ 1. The operation of contraction (2.5) is applicable to 
tensorial and pseudotensorial fields. It has the following properties: 

(1) the contraction of a tensorial field is a tensorial field; 

(2) the contraction of a pseudotensorial field is a pseudotensorial field. 

The operation of contraction extended to the case of pseudotensorial fields preserve 
its linearity given by the equalities (3.4). 

The covariant differentiation of pseudotensorial fields in a Cartesian coordinate 
system is determined by the formula (5.1). The covariant differential VA of a 
tensorial field is a tensorial field; the covariant differential of a pseudotensorial 
field is a pseudotensorial field. It is convenient to express the properties of the 
covariant differential through the properties of the covariant derivative Vx in the 
direction of a field X. Now X is either a vectorial or a pseudovectorial field. All 
propositions of the theorem 5.2 for Vx remain valid. 

§ 8. A note on the orientation. 

Pseudoscalar fields form a particular case of pseudotensorial fields. Scalar 
fields can be interpreted as functions whose argument is a point of the space 
E. In this interpretation they do not depend on the choice of a coordinate 
system. Pseudoscalar fields even in such interpretation preserve some dependence 
on a coordinate system, though this dependence is rather weak. Let ^ be a 
pseudoscalar field. In a fixed Cartesian coordinate system the field ^ is represented 
by a scalar function £,{P), where P S E. The value of this function ^ at a point P 
does not change if we pass to another coordinate system of the same orientation, 
i. e. if the determinant of the transition matrix S is positive. When passing to 
a coordinate system of the opposite orientation the function ^ changes the sign: 
£,{P) —£,{P). Let's consider a nonzero constant pseudoscalar field ^. In some 
coordinate systems ^ = c = const, in others ^ = — c = const. Without loss of 
generality we can take c — 1. Then such a pseudoscalar field ^ can be used 
to distinguish the coordinate systems where ^ = 1 from those of the opposite 
orientation where ^ = — 1. 

Proposition. Defining a unitary constant pseudoscalar Geld ^ is equivalent to 
choosing some preferable orientation in the space E. 

From purely mathematical point of view the space E, which is a three- 
dimensional Euclidean point space (see definition in [1]), has no preferable 
orientation. However, the real physical space E (where we all live) has such 
an orientation. Therein we can distinguish the left hand from the right hand. This 
difference in the nature is not formal and purely terminological: the left hemi- 
sphere of a human brain is somewhat different from the right hemisphere in its 
functionality, in many substances of the organic origin some isomers prevail over 
the mirror symmetric isomers. The number of left-handed people and the number 
of right-handed people in the mankind is not fifty-fifty as well. The asymmetry of 
the left and right is observed even in basic forms of the matter: it is reflected in 
modern theories of elementary particles. Thus, we can assume the space E to be 
canonically equipped with some pseudoscalar field £,e whose values are given by 

CopyRight © Shaiipov R.A., 1996, 2004. 



36 



CHAPTER II. ELEMENTS OF TENSORIAL ANALYSIS. 



the formula 

r 1 in right-oriented coordinate systems, 
^ \ — 1 in left-oriented coordinate systems. 

Multiplying by S,e, we transform a tensorial field A into the pseudotensorial 
field (g) A = • A. Multiplying by once again, we transform ■ A back to 
A. Therefore, in the space E equipped with the preferable orientation in form of 
the field S,e one can not to consider pseudotensors at all, considering only tensor 
fields. The components of the volume tensor u: in this space should be defined as 

(^i]k = S.E Vdetg £yfe. (8.1) 

Let X and Y be two vectorial fields. Then we can define the vectorial field Z 
with the following components: 

3 3 3 

^"^EEEs'^^u-fe^^'^'- (8.2) 

j=l j=l k=l 

From (8.2), it is easy to sec that Z is derived as the contraction of the field 
g(8)Ci;(g)a(8)b. In a rectangular Cartesian coordinate system with right-oriented 
orthonormal basis the formula (8.2) takes the form of the well-known formula for 
the components of the vector product Z = [X, Y] (see [4] and the formula (6.9) 
above). In a space without preferable orientation, where ujijk is given by the 
formula (6.fl), the vector product of two vectors is a pseudovector. 

Now let's consider three vectorial fields X, Y, and Z and let's construct the 
scalar field u by means of the following formula: 

3 3 3 

1=1 j=i k=i 

When (8.3) is written in a rectangular Cartesian coordinate system with right- 
oriented orthonormal basis, one easily sees that the field (8.3) coincides with the 
mixed product (X, Y Z) of three vectors (see [4] and the formula (6.9) above). 
In a space without preferable orientation the mixed product of three vector fields 
determined by the volume pscudotcnsor (6. II) appears to be a pseudoscalar field. 

§ 9. Raising and lovirering indices. 

Let A be a tensor field or a pseudotensor field of the type (r, s) in the space E 
and let r ^ I. Let's construct the tensor product A (g) g of the field A and the 
metric tensor g, then define the field B of the type (r — I, s + I) as follows: 

5;i:::;:;: = E41:::r^^:::^;\5.,„. O-d 

fe=i 

The passage from the field A to the field B according to the formula (9.1) is called 
the index lowering procedure of the m-th upper index to the n-th lower position. 
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Using the inverse metric tensor, one can invert the operation (9.1). Let B be a 
tensorial or a pseudotensorial field of the type (r, s) and let s ^ 1. Then we define 
the field A = C(B (g) g) of the type (r + 1, s — 1) according to the formula: 

3 

^n-j.-i- 2^^n-j„-iqi„-js-i 9 ■ (9.2) 

9=1 

The passage from the field B to the field A according to the formula (9.2) is called 
the index raising procedure of the n-th lower index to the m-th upper position. 

The operations of lowering and raising indices are inverse to each other. Indeed, 
we can perform the following calculations: 



I . . . im. — l k ijjiJ^l ... ir 



Jrl•■■^r Ql, 

9kq 



3-n.-l]n---3s 

g=l k=l 
3 

E,ii...i™_ifci™ + i...V sim _ Ail.. A,. 

k=l 



The above calculations show that applying the index lowering and the index 
raising procedures successively A ^ B ^ C, we get the field C = A. Applying 
the same procedures in the reverse order yields the same result. This follows from 
the calculations just below: 

3 3 

- 2^ 2^^h-3^-iq3^ + i-3s 9 9k]„ - 
k=l q=l 

3 

E,Jl...i„_li„...i, eg ^ 

jl... q + js "in 31... Jb' 

9=1 

The existence of the index raising and index lowering procedures follows from 
the very nature of the space E which is equipped with the scalar product and, 
hence, with the metric tensors g and g. Therefore, any tensorial (or pseudotenso- 
rial) field of the type (r, s) in such a space can be understood as originated from 
some purely covariant field of the type (0, r + s) as a result of raising a part of its 
indices. Therein a little bit different way of setting indices is used. Let's consider 
a field A of the type (0,4) as an example. We denote by Ai-^^i^i^i^ its components 
in some Cartesian coordinate system. By raising one of the four indices in A one 
can get the four fields of the type (1, 3). Their components are denoted as 

Ail 4 *2 4 *3 /I '4 (Q "i) 

i2i3U' il isU' ili2 U' 11^213 ' W'^J 

Raising one of the indices in (9.3), we get an empty place underneath it in the 
list of lower indices, while the numbering of the indices at that place remains 
unbroken. In this way of writing indices, each index has «its fixed position* no 
matter what index it is — a lower or an upper index. Therefore, in the writing 
below we easily guess the way in which the components of tensors are derived: 

^^jkq^ ^\q^ ^^3^q^ ^/ q- (9-4) 
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Despite to some advantages of the above form of index setting in (9.3) and (9.4), 
it is not commonly admitted. The matter is that it has a number of disadvantages 
either. For example, the writing of general formulas (1.6), (2.2), (2.5), and some 
others becomes huge and inconvenient for perception. In what follows we shall not 
change the previous way of index setting. 

§ 10. Gradient, divergency, and rotor. 
Some identities of the vectorial analysis. 

Let's consider a scalar field or, in other words, a function /. Then apply the 
operator of covariant differentiation V, as defined in §5, to the field /. The 
covariant differential V/ is a covectorial field (a field of the type (0, 1)). Applying 
the index raising procedure (9.2) to the covector field V/, we get the vector field 
F. Its components are given by the following formula: 

fe=i 

Definition 10.1. The vector field F in the space E whose components are 
calculated by the formula (10.1) is called the gradient of a function /. The 
gradient is denoted as F = grad/. 

Let X be a vectorial field in E. Let's consider the scalar product of the vectorial 
fields X and grad/. Due to the formula (10.1) such scalar product of two vectors 
is reduced to the scalar product of the vector X and the covector V/: 

(X| grad/)==^X'^|4 = (W|X). (10.2) 

fc=i 

The quantity (10.2) is a scalar quantity. It does not depend on the choice of a 
coordinate system where the components of X and V/ are given. Another form of 
writing the formula (10.2) is due to the covariant differentiation along the vector 
field X, it was introduced by the formula (5.5) above: 

(X|grad/) = Vx/. (10.3) 

By analogy with the formula (10.3), the covariant differential VF of an arbitrary 
tensorial field F is sometimes called the covariant gradient of the field F. 

Let F be a vector field. Then its covariant differential VF is an operator field, 
i. e. a field of the type (1, 1). Let's denote by ip the contraction of the field VF: 

v' = ^(VF) = Ea^- (10-4) 

k=l 

Definition 10.2. The scalar field (p in the space E determined by the formula 
(10.4) is called the divergency of a vector field F. It is denoted (f = divF. 

Apart from the scalar field divF, one can use VF in order to build a vector 
field. Indeed, let's consider the quantities 

z'" = E E E E 3™ ^^^^ v,F^ (10.5) 

i=l j = l k=l q=l 
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where Uijk are the components of the volume tensor given by the formula (8.1). 

Definition 10.3. The vector field p in the space E determined by the formula 
(10.5) is called the rotor^ of a vector field F. It is denoted p = rotF. 

Due to (10.5) the rotor or a vector field F is the contraction of the tensor field 
g (g) w (g) g (g) VF with respect to four pairs of indices: rot F = C(gCS)u;Cg)gCg) VF). 

Remark. If uiijk in (10.5) are understood as components of the volume pseu- 
dotcnsor (6.11), then the rotor of a vector field should be understood as a 
pseudovectorial field. 

Suppose that O, ei, 62, 63 is a rectangular Cartesian coordinate system in E 
with orthonormal right-oriented basis ei, 62, 63. The Gram matrix of the basis 
ei, 62, 63 is the unit matrix. Therefore, we have 



9ij 



1 for i = j, 
for i 7^ j. 



The pseudoscalar field defining the orientation in E is equal to unity in a right- 
oriented coordinate system: = ^- Due to these circumstances the formulas 
(10.1) and (10.5) for the components of grad/ and rotF simplifies substantially: 



3 3 

(rotF)'=^^£,,, 



(10.6) 
(10.7) 



j=i fc=i 

The formula (10.4) for the divergency remains unchanged: 

dF^ 



divF = 



k=l 



(10.8) 



The formula (10.7) for the rotor has an elegant representation in form of the 
determinant of a 3 x 3 matrix: 



rotF = 



ei 


62 


63 


d 


d 


d 




dx"^ 


dx^ 




F^ 


F^ 



(10.9) 



The formula (8.2) for the vector product in right-oriented rectangular Cartesian 
coordinate system takes the form of (6.9). It can also be represented in the form 
of the formal determinant of a 3 x 3 matrix: 



[X, Y] 



61 62 63 
yl y2 ^3 



(10.10) 



^ The term «curl» is also used for the rotor. 
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Due to the similarity of (10.9) and (10.10) one can formally represent the operator 
of covariant differentiation V as a vector with components d/dx^, d/dx^, d/dx^. 
Then the divergency and rotor are represented as the scalar and vectorial products: 

divF = (V|F), rotF=[V, F]. 

Theorem 10.1. For any scalar field Lp of the smoothness class the equality 
rot grad = is identically fulfilled. 

Proof. Let's choose some right-oriented rectangular Cartesian coordinate sys- 
tem and then use the formulas (10.6) and(10.7). Let F = rot grad 93. Then 

j=i k=i 

Let's rename the summation indices in (10.11). The index j is replaced by the 
index k and vice versa. Such a swap of indices does not change the value of the 
sum in (10.11). Therefore, we have 

3 3 3 3 ^9 

j=l k=l j=l k=l 

Here we used the skew-symmetry of the Levi-Civuta symbol with respect to the 
pair of indices j and k and the symmetry of the second order partial derivatives of 
the function (p with respect to the same pair of indices: 

- . . (10.12) 

dx^ dx^ dx^ dx^ 

For class functions the value of second order partial derivatives (10.12) do not 
depend on the order of differentiation. The equality i^' = — now immediately 
yields F' = 0. The theorem is proved. □ 

Theorem 10.2. For any vector field F of the smoothness class the equality 
div rot F = is identically fulfilled. 

Proof. Here, as in the case of the theorem 10.1, we choose a right-oriented 
rectangular Cartesian coordinate system, then we use the formulas (10.7) and 
(10.8). For the scalar field = div rot F from these formulas we derive 

333 F^ 

^ = EEE^.-a^- (10-13) 

i=i j=i k=i 

Using the relationship analogous to (10.12) for the partial derivatives 

q2 pk q2 pk 



dxi dx^ dx"^ dxi 

and using the skew-symmetry of £ijk with respect to indices i and j, from the 
formula (10.13) we easily derive ip = —tf. Hence, ip = Q. □ 
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Let (yC be a scalar field of tlie smoothness class . The quantity divgradiyS in 
general case is nonzero. It is denoted /\Lp — div grades. The sign A denotes the 
differential operator of the second order that transforms a scalar field Lp to another 
scalar field div grad (/?. It is called the Laplace operator or the laplacian. In a 
rectangular Cartesian coordinate system it is given by the formula 

Using the formulas (10.6) and (10.8) one can calculate the Laplace operator in a 
skew-angular Cartesian coordinate system: 

A^yy g'^ (10.15) 

i=i j=i 

Using the signs of covariant derivatives Vi — d/dx^ we can write the Laplace 
operator (10.15) as follows: 



i=i j=i 

The equality (10.16) differs from (10.15) not only in special notations for the 
derivatives. The Laplace operator defined as Atp = div grad <p can be applied only 
to a scalar field (p. The formula (10.16) extends it, and now we can apply the 
Laplace operator to any twice continuously diffcrcntiable tensor field F of any 
type (r, s). Due to this formula AF is the result of contracting the tensor product 
g (8) VVF with respect to two pairs of indices: AF = C(g (g) VVF). The resulting 
field AF has the same type (r, s) as the original field F. The laplace operator in 
the form of (10.16) is sometimes called the Laplace- Beltrami operator. 

§ 11. Potential and vorticular vector fields. 

Definition 11.1. A differentiable vector field F in the space E is called a 
potential field if rotF = 0. 

Definition 11.2. A differentiable vector field F in the space E is called a 
vorticular field if divF = 0. 

The theorem 10.1 yields some examples of potential vector fields, while the 
theorem 10.2 yields the examples of vorticular fields. Indeed, any field of the form 
grad</J is a potential field, and any field of the form rot F is a vorticular one. As it 
appears, the theorems 10.1 and 10.2 can be strengthened. 

Theorem 11.1. Any potential vector Geld F in the space E is a gradient of 
some scalar field (p, i. e. F = giadip. 

Proof. Let's choose a rectangular Cartesian coordinate system O, ei, 82, 63 
with orthonormal right-oriented basis ei, 62, 33. In this coordinate system the 
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potentiality condition rot F = for the vector field F is equivalent tho the following 
three relationships for its components: 

(11.1) 
(11.2) 
(11.3) 

The relationships (11.1), (11-2), and (11.3) are easily derived from (10.7) or from 
(10.9). Let's define the fimction <^{x^ x'^, x^) as the sum of three integrals: 

xi 

</3(a;^x^a;3) = c+ / F\x\0,0)dx^ + 



dF\x\x^,x^) 


dF^{x\x^ 


x^) 


dx'^ 


dx^ 




dF\x\x^,x^) 


dF^{x\x^ 


x') 


dx^ 


dx"^ 




dF^{x\x^,x'^) 


dF^{x\x^ 


x') 


dx^ 


dx^ 





(11.4) 

+ / F'^{x\x\0)dx^ + I F^{x^,x^,x^)dx^. 





Here c is an arbitrary constant. Now we only have to check up that the function 
(11.4) is that very scalar field for which F = gradiyS. 

Let's differentiate the function ip with respect to the variable x^. The constant 
c and the first two integrals in (11.4) do not depend on x^. Therefore, we have 

^ = ^1 F^{x\x\x')dx^=F\x\x\x% (11.5) 



In deriving the relationship (11.5) wc used the rule of differentiation of an integral 
with variable upper limit (sec [2]). 

Now let's differentiate the function (y9 with respect to x^ . The constant c and 
the first integral in (11.4) docs not depend on x^ . Differentiating the rest two 
integrals, we get the following expression: 

g = ^^^(.^x^o) + A j F\x\x\x^)dx\ 



The operations of differentiation with respect to x'^ and integration with respect 
to x^ in the above formula are commutative (see [2]). Therefore, we have 

^ = F^.'.^,,),J ^''%f^''\ .'. (11.6) 



In order to transform the expression being integrated in (11.6) we use the formula 
(11.2). This leads to the following result: 

(11.7) 

X=X3 

= F^{x\x^,0) + F'^{x\x^,x) = F^{x\x'^,x^). 
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In calculating the derivative dip/dx^ we use that same tricks as in the case of the 
other two derivatives dip/dx^ and dip/dx^: 

XI X2 

^ = A| Fi(.i,0,0)dx^ + A J F\x\x^O)dx^+ 





X2 



dx^ J dx^ 





To transform the last two integrals we use the relationships (11.1) and (11.3) 

X=X2 

+ 

x=0 



^ = F\x\0,0)+F\x\x,0) 
+ F'{x\x'',x) 



X=X3 

1^-1-2^^ ^F\x\x^x'). 



(11. 



x=0 



The relationships (11.5), (11.7), and (11.8) show that grad(y3 = F for the function 
(p(x^ , x'^ , x^) given by the formula (11.4). The theorem is proved. □ 

Theorem 11.2. Any vorticular vector field F in the space E is the rotor of 
some other vector field A, i. e. F = rot A. 

Proof. We perform the proof of this theorem in some rectangular Cartesian 
coordinate system with the orthonormal basis ei, 62, 63. The condition of vorticity 
for the field F in such a coordinate system is expressed by a single equation: 

dF\^) ^ dF\^) ^ dF\^) 
dx^ dx"^ dx^ 

Let's construct the vector field A defining its components in the chosen coordinate 
system by the following three formulas: 

X3 X2 

„2 



= J F'(x^x^x3)dT3- y F\x\x'',0)dx'' 



= 0. 

Let's show that the field A with components (11.10) is that very field for which 
rot A = F. We shall do it calculating directly the components of the rotor in the 
chosen coordinate system. For the first component we have 
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Here we used the rule of differentiation of an integral with variable upper limit. In 
calculating the second component we take into account that the second integral in 
the expression for the component in (11.10) does not depend on x^: 

'^-li^-ir, l^v.^-yxi^' -Fv,^'yy 



And finally, for the third components of the rotor we derive 



dA^ dA^ 



dx^ dx^ 
d 



fdF\x\x^,x^) dF^{x\x^,x^) 



dx^ 



X2 X3 



F%x\x,0)dx^ 



dx"^ 



dF^{x\x^x^) 



dx-^ 



dx^ 



dx-^ + 



F\x\x\0) = F\x\x\x) 



+ F%x\x^,0) = F%x\x^,x-^) 



In these calculations we used the relationship (11.9) in order to replace the sum of 
two partial derivatives dF^/dx^ +dF^/dx^ by —dF^/dx^. Now, bringing together 
the results of calculating all three components of the rotor, we see that rot A = F. 
Hence, the required field A can indeed be chosen in the form of (11.10). □ 



CHAPTER III 



CURVILINEAR COORDINATES 



§ 1. Some examples of curvilinear coordinate systems. 

The main purpose of Cartesian coordinate systems is the numeric representation 
of points: each point of the space E is represented by some unique triple of numbers 
Curvihnear coordinate systems serve for the same purpose. We begin 
considering such coordinate systems with some examples. 

Polar coordinates. Let's consider a plane, choose some point O on it (this will 
be the pole) and some ray OX coming out from this point. For an arbitrary 
point ^ ^ O of that plane its position is determined by two parameters: the 




Fig. 1.1 Fig. 1.2 

length of its radius- vector p = \oA\ and the value of the angle (p between the ray 
OX and the radius-vector of the point A. Certainly, one should also choose a 
positive (counterclockwise) direction to which the angle is laid (this is equivalent 
to choosing a preferable orientation on the plane). Angles laid to the opposite 
direction are understood as negative angles. The numbers p and (p are called the 
polar coordinates of the point A. 

Let's associate some Cartesian coordinate system with the polar coordinates as 
shown on Fig. 1.2. We choose the point O as an origin, then direct the abscissa 
axis along the ray OX and get the ordinate axis from the abscissa axis rotating it 
by 90°. Then the Cartesian coordinates of the point A are derived from its polar 
coordinates by means of the formulas 



= p cos{ip), 
= p sin((p). 



(1.1) 
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Conversely, one can express p and ip through and as follows: 



P^^{x^f + [x^)\ 
If = arctan(x^/a;"'^). 



(1.2) 



The pair of numbers (p, (p) can be treated as an element of the two-dimensional 



space 



In order to express M visually we represent this space as a coordi- 
nate plane. The coordinate plane (p, ip) 
has not its own geometric interpretation, 
it is called the map of the polar coordi- 
nate system. Not all points of the map 
correspond to the real geometric points. 
The condition p excludes the whole 
half-plane of the map. The sine and 
cosine both are periodic functions with 
the period 2tt — 360°. Therefore there 
are different points of the map that re- 
present the same geometric point. Thus, 
the mapping {p,(p) — > {x^,x'^) given by 
the formulas (1.1) is not injcctive. Let 
U be the unbounded domain highlighted 
with the light blue color on Fig. 1.3 (the 
points of the boundary are not included). Denote by V the image of the domain 
U under the mapping (1.1). It is easy to understand that V is the set of all points 
of the (x^, x"^) plane except for those lying on the ray OX. If wc restrict the 
mapping (1.1) to the domain U , we get the bijective mapping m : U V . 

Note that the formula (1.2) is not an exact expression for the inverse mapping 
TO~^ : V U. The matter is that the values of tan((^) at the points (a;^,a;^) and 
{—x^, —x^) do coincide. In order to express exactly it would be better to use 
the tangent of the half angle: 




Fig. 1.3 



tan((^/2) = . 

x^ + J [x^f + [x^f 



However, we prefer the not absolutely exact expression for ip from (1.2) since it is 
relatively simple. 

Let's draw the series of equidistant straight lines parallel to the axes on the 
map of the polar coordinate system (see Fig. 1.4 below). The mapping (1.1) 
takes them to the series of rays and concentric circles on the (.t^, x^) plane. 
The straight lines on Fig. 1.4 and the rays and circles on Fig. 1.5 compose 
the coordinate network of the polar coordinate system. By reducing the intervals 
between the lines one can obtain a more dense coordinate network. This procedure 
can be repeated infinitely many times producing more and more dense networks 
in each step. Ultimately (in the continuum limit), one can think the coordinate 
network to be maximally dense. Such a network consist of two families of lines: the 
first family is given by the condition ip = const, the second one — by the similar 
condition p = const. 



§ 1. SOME EXAMPLES OF CURVILINEAR COORDINATE SYSTEMS. 



47 



On Fig. 1.4 exactly two coordinate lines pass through each point of the map: 
one is from the first family and the other is from the second family. On the 
{x^,x'^) plane this condition is fulfilled at all points except for the origin O. Here 
all coordinate lines of the first family are crossed. The origin O is the only singular 
point of the polar coordinate system. 




Fig. 1.4 



Fig. 1.5 



The cylindrical coordinate system in the space E is obtained from the polar 
coordinates on a plane by adding the third coordinate h. As in the case of 

polar coordinate system, we associate some 
Cartesian coordinate system with the cylin- 
drical coordinate system (see Fig. 1.6). Then 




Fig. 1.6 



x^ — p cos{(p), 
x"^ ^ p sin((/?), 
x'^ — h. 



(1.3) 



Conversely, one can pass from Cartesian to 
cylindrical coordinates by means of the formu- 
la analogous to (1.2): 



' p^^{x^f + {x^)\ 
if ~ arctan(x^/a;"'^), 
h - x\ 



(1.4) 



The coordinate network of the cylindrical coordinate system consists of three 
families of lines. These are the horizontal rays coming out from the points of the 
vertical axis Ox^ , the horizontal circles with the centers at the points of the axis 
Ox^ , and the vertical straight lines parallel to the axis Ox^ . The singular points 
of the cylindrical coordinate system fill the axis Ox^ . Exactly three coordinate 
lines (one from each family) pass through each regular point of the space E, i. e. 
through each point that does not lie on the axis Ox^ . 
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A 



The spherical coordinate system in the spa- 
ce E is obtained by shght modification of 
the cyhndrical coordinates. The coordinate 
h is replaced by the angular coordinate i?, 
while the quantity p in spherical coordinates 
denotes the length of the radius-vector of the 
point A (see Fig. 1.7). Then 



X 



1 



p sin(i9) cos((y5), 
p sin('(?) sin((^), 
p cos{-d). 



X' 



.2 



(1.5) 



X- 



.3 



Fig. 1.7 



The spherical coordinates of a point are usual- 
ly written in the following order: p is the first 
coordinate, i9 is the second one, and Lp is 



the third coordinate. The converse transition 
from Cartesian coordinates to these quantities is given by the formula: 



Coordinate lines of spherical coordinates form three families. The first family is 
composed of the rays coming out from the point O; the second family is formed 
by circles that lie in various vertical planes passing through the axix Ox^\ and 
the third family consists of horizontal circles whose centers are on the axis Ox^ . 
Exactly three coordinate lines pass through each regular point of the space E, one 
line from each family. 

The condition p — const specifies the sphere of the radius p in the space E. The 
coordinate lines of the second and third families define the network of meridians 
and parallels on this sphere exactly the same as used in geography to define the 
coordinates on the Earth surface. 



Let D be some domain in the space E. Saying domain, one usually understand 
a connected open set. An open set D means that along with each its point 
A ^ D the set D comprises some spherical neighborhood 0(A) of this point. A 
connected set D means that any two points of this set can be connected by a 
smooth curve lying within D. See more details in [2]. Let's consider three numeric 
functions w^(x), u^(x), and it'^(x) defined in the domain D. Generally speaking, 
their domains could be wider, but we need them only within D. The values of 
three functions m^, u^, at each point form a triple of numbers, they can be 
interpreted as a point of the space . Then the triple of functions u^, u^, v?" 
define a mapping u : D . 

Definition 2.1. A triple of differentiable functions it^, u^, u'^ is called regular 
at a point A of the space E if the gradients of these functions gradu^, gradu^, and 
gradw'^ are linearly independent at the point A. 




(1.6) 



§ 2. Moving frame of a curvilinear coordinate system. 



§ 2. MOVING FRAME OF A CURVILINEAR COORDINATE SYSTEM. 
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Let's choose some Cartesian coordinate system in E, in this coordinate system 
the above functions u^, u^, are represented by the functions — 
of Cartesian coordinates of a point. The gradients of the differentiable functions 
w^, u^, form the triple of covectorial fields whose components are given by the 
partial derivatives of u^, u^, u'^ with respect to 



The matrix J of the form (2.2) is called the Jacobi matrix of the mapping 
u: D ^ given by the triple of the differentiable functions u^, v? , in the 
domain D. It is obvious that the regularity of the functions u^, v? , at a point 
is equivalent to the non-degeneracy of the Jacobi matrix at that point: det J 7^ 0. 

Theorem 2.1. If continuously differentiable functions u^, u^, with the do- 
main D are regular at a point A, then there exists some neighborhood 0{A) of the 
point A and some neighborhood 0{u{A)) of the point u{A) in the space M.^ such 
that the following conditions are fulhlled: 

(1) the mapping u: 0{A) 0{u{A)) is bijective; 

(2) the inverse mapping : 0{u{A)) 0{A) is continuously differentiable. 

The theorem 2.1 or propositions equivalent to it are usually proved in the course 
of mathematical analysis (see [2]). They are known as the theorems on implicit 
functions. 

Definition 2.2. Say that an ordered triple of continuously differentiable func- 
tions li^, u^, with the domain D (Z E define a curvilinear coordinate system in 
D if it is regular at all points of D and if the mapping u determined by them is a 
bijective mapping from D to some domain t/ C M'^. 

The cylindrical coordinate system is given by three functions ~ p{^)i 
= 'p(x), and = /i(x) from (1.4), while the spherical coordinate system is 
given by the fimctions (1.6). However, the triples of functions (1.4) and (1.6) 
satisfy the conditions from the definition 2.2 only after reducing somewhat their 
domains. Upon proper choice of a domain D for (1.4) and (1.6) the inverse 
mappings are given by the formulas (1.3) and (1.5). 

Suppose that in a domain Z) C E a curvilinear coordinate system u^, u^, is 
given. Let's choose an auxiliary Cartesian coordinate system in E. Then u^, u^, 
is a triple of functions defining a map u from D onto some domain [/ C M'^: 




(2.1) 



Let's compose a matrix of the gradients (2.1): 



dx^ dx"^ 

dv? du^ 

dx^ dx"^ 

dx^ dx"^ 



dv? 



(2.2) 




(2.3) 
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The domain D is called the domain being mapped, the domain L/ C K'^ is called 
the map or the chart, while u^^ : U D is called the chart mapping. The chart 
mapping is given by the following three functions: 

1 1/1 2 3\ 

X ~ X (U , U , U ), 

x^ = x'^{u^,u'^,u^), (2.4) 
a;3 — x'^{u^ ,u'^ ,u^). 



Denote by r the radius- vector r of the point with Cartesian coordinates x^, x^ , x^ . 
Then instead of three scalar functions (2.4) we can use one vectorial function 



3 

V{U^,U^,U^) =^X''{U^,U^,U^) ■ Bq. (2.5) 
9=1 



Let's fix some two of three coordinates u^, , and let's vary the third of them. 
Thus we get three families of straight lines within the domain [/ C M^^: 




Here c^, c^, are constants. The straight lines (2.6) form a rectangular coordinate 
network within the chart U . Exactly one straight line from each of the families 
(2.6) passes through each point of the chart. Substituting (2.6) into (2.5) we map 
the rectangular network from U onto a curvilinear network in the domain D <Z E. 
Such a network is called the coordinate network of a curvilinear coordinate system. 

The coordinate network of a curvilinear coordinate system on the domain D 
consists of three families of lines. Due to the bijectivity of the mapping u : D ^ U 
exactly three coordinate lines pass through each point of the domain D — one line 
from each family. Each coordinate line has its canonical parametrization: t = 
is the parameter for the lines of the first family, t = is the parameter for the 
lines of the second family, and finally, t = u'^ is the parameter for the lines of 
the third family. At each point of the domain D we have three tangent vectors, 
they are tangent to the coordinate lines of the three families passing through that 
point. Let's denote them Ei, E2, E3. The vectors Ei, E2, E3 are obtained by 
differentiating the radius-vector r(u^,u^,u'^) with respect to the parameters u^, 
M^, of coordinate lines. Therefore, we can write 



Ej(M\u^it3) = 



(2.7) 



Let's substitute (2.5) into (2.7). The basis vectors ei, 62, 63 do not depend on the 
variables u^, m^, u^, hence, we get 

E,(.^.^.3)^^M!W).^^. (2.8) 

q=l 
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The formula (2.8) determines the expansion of the vectors Ei, E2, E3 in the basis 
ei, 62, 63. The column- vectors of the coordinates of Ei, E2, and E3 can be 
concatenated into the following matrix: 







dx^ 


du^ 


du^ 


du^ 


dx'^ 


dx'^ 


dx'^ 


du^ 


du^ 


du^ 


dx^ 


dx^ 


dx^ 









(2.9) 



Comparing (2.9) and (2.2). we see that (2.9) is the Jacobi matrix for the mapping 
: U ^ D given by the functions (2.4). Let's substitute (2.4) into (2.3): 

u'{x^{u^,u^,u^),x^{u^,u^,u^),x^{u\u^,u^)) = u\ (2.10) 

The identity (2.10) follows from the fact that the functions (2.3) and (2.4) define 
two mutually inverse mappings u and u~^. Let's differentiate the identity (2.10) 
with respect to the variable u-': 



dxi 



(2.11) 



Here we used the chain rule for differentiating the composite function in (2.10). 
The relationship (2.11) shows that the matrices (2.2) and (2.9) are inverse to each 
other. More precisely, we have the following relationship 



/(w\u^u3) = J(a;\ 



X , X 



(2.12) 



by means of (2.4), or 
x^, x^ 



where ^ should be expressed through u^, 

conversely, u^, u^, should be expressed through x^, x^, x'^ by means of (2.3). 
The arguments shown in the relationship (2.12) are the natural arguments for the 
components of the Jacobi matrices / and J. However, one can pass to any required 
set of variables by means of (2.3) or (2.4) whenever it is necessary. 

The regularity of the triple of functions (2.3) defining a curvilinear coordinate 
system in the domain D means that the matrix (2.2) is non-degenerate. Then, due 
to (2.12), the inverse matrix (2.9) is also non-degenerate. Therefore, the vectors 
El, E2, E3 given by the formula (2.8) are linearly independent at any point of 
the domain D. Due to the linear independence of the coordinate tangent vectors 
El, E2, E3 they form a moving frame which is usually called the coordinate frame 
of the curvilinear coordinate system. The formula (2.8) now can be interpreted 
as the transition formula for passing from the basis of the auxiliary Cartesian 
coordinate system to the basis formed by the vectors of the frame Ei, E2, E3: 



E, 



3 

E 

9=1 



S%u\u\u'') 



(2.13) 



The transition matrix S in the formula (2.13) coincides with the Jacobi matrix 
(2.9), therefore its components depend on u^, u^, . These are the natural 
variables for the components of S. 
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The inverse transition from the basis Ei, E2, E3 to the basis of the Cartesian 
coordinate system is given by the inverse matrix T = S^^. Due to (2.12) 
the inverse transition matrix coincides with the Jacobi matrix (2.2). Therefore, 
x^, x^, x'^ are the natural variables for the components of the matrix T: 

3 

e, = ^T^(x\x2,x3).E,. (2.14) 

i=l 

The vectors Ei, E2, E3 of the moving frame depend on a point of the domain 
D G E. Since in a curvilinear coordinate system such a point is represented by its 
coordinates u^, u^, u^, these variables are that very arguments which are natural 
for the vectors of the moving frame: E^ = Ei(u^, m^, u^). 

§ 3. Change of curvilinear coordinates. 

Let u^, V? , V? be some curvilinear coordinates in some domain Di and let 
v} , , be some other curvilinear coordinates in some other domain D2. If 
the domains Di and D2 do intersect, then in the domain D = Di O D2 we have 
two coordinate systems. We denote by U and U the preimages of the domain 




Fig. 3.1 

D in the maps Ui and U2; i.e. we denote U = u{Di n D2) and we denote 
U = u{Di O 02). Due to the chart mappings the points of the domain D are in 
one-to-one correspondence with the points of the domains U and U. As for the 
chart mappings u^^ and u^^, they are given by the following hmctions: 

= x^{u^,u^,u'^), ( x^ = x^{u^,u'^, it^), 

x^ ^ x^ {u^ , u'^ , u^) , < x^ = x^ {u^ , u'^ , u^) , (3.1) 

— x^{u^,u'^, u'^), [ x'^ = x^{u^,u'^, u^). 

The mappings u and u inverse to the chart mappings are given similarly: 

= u^{x^,x'^,x^), ( = (x^ , x'^ , x^) , 

^ u'^{x^,x'^,x^), < u"^ = u'^{x^,x^,x^), (3.2) 

u'^(x^,x'^,x^), [ u'^ = u'^{x^,x'^,x^). 

Let's substitute the first set of functions (3.1) into the arguments of the second set 
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of functions (3.2). Similarly, we substitute the second set of functions (3.1) into 
the arguments of the first set of functions in (3.2). As a result we get the functions 



{x^ {u^ , u'^ , u^), x'^{v} ,v? ,u^), x^{u^ ,u'^ ,u^)) 

{x^ {u^ , ii'^ , u^), x^ {u^ , u'^ , u^) , x^{u^ ,u'^ ,u^)) 
u^(x^(u\ u^, u^), x^(u\ u^, u^), a;^(u\ u^, u^)) 
u^(x^({t\ u^, u^), x^(u\ u^, u^), x^(u\ u^, u^)) 



(3.3) 



(3.4) 



which define the pair of mutually inverse mappings u » u ^ and u » u ^ . For the 
sake of brevity we write these sets of functions as follows: 

= u'^ {u^ , u'^ , u^) , < u^(u\ u^, u^), (3.5) 

The formulas (3.5) express the coordinates of a point from the domain D in some 
curvilinear coordinate system through its coordinates in some other coordinate 
system. These formulas are called the transformation formulas or the formulas for 
changing the curvilinear coordinates. 

Each of the two curvilinear coordinate systems has its own moving frame within 
the domain D = Di nD2- Let's denote by S and T the transition matrices relating 
these two moving frames. Then we can write 



E,=^S;;.E„ E, = ^7f.Efc. (3.6) 

i=l fc=l 

Theorem 3.1. The components of the transition matrices S and T for the 
moving frames of two curviUnear coordinate system in (3.6) are determined by the 
partial derivatives of the functions (3.5): 

= Ttin\u^u^) = ^. (3.7) 

Proof. We shall prove only the first formula in (3.7). The proof of the second 
formula is absolutely analogous to the proof of the first one. Let's choose some 
auxiliary Cartesian coordinate system and then write the formula (2.8) applied to 
the frame vectors of the second curvilinear coordinate system: 

E,(.\.^.3)^^MW!).^^_ (3 3) 

9=1 

Applying the formula (2.14), we express ei, 62, 63 through Ei, E2, E3. Remember 
that the matrix T in (2.14) coincides with the Jacobi matrix J{x^,x^,x^) from 
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(2.2). Therefore, we can write the foUowmg formula: 



dx'i 



E,. (3.9) 



Now let's substitute (3.9) into (3.8). As a result we get the formula relating the 
frame vectors of two curvilinear coordinate systems: 

E, = y f y 9uK-'^^'^^') dx^{u\u^u^)\ 
^ ^ \ ^ dxi dui j 

1=1 \ q=l / 

Comparing (3.10) and (3.6), from this comparison for the components of S we get 

du'{x^,x^,x^) 9x'?(m\m^u3) 

Remember that the Cartesian coordinates x^ , x^ , x^ in the above formula (3.11) 
are related to the curvilinear coordinates v}, u^, by means of (3.1). Hence, the 
sum in right hand side of (3.11) can be transformed to the partial derivative of the 
composite function m*((x"'^(m"'^, u^, w'^), x^(i2^, u^, w'^), x^({t^, u^, u'^)) from (3.4): 

= ^ 

^ du^ ■ 

Note that the functions (3.4) written in the form of (3.5) are that very functions 
relating {t^, {t^, and m^, m^, u^, and their derivatives are in formula (3.7). The 
theorem is proved. □ 

A remark on the orientation. From the definition 2.2 we derive that the 
functions (2.3) are continuously differentiable. Due to the theorem 2.1 the func- 
tions (2.4) representing the inverse mappings are also continuously differentiable. 
Then the components of the matrix S in the formula (2.13) coinciding with the 
components of the Jacobi matrix (2.9) are continuous functions within the domain 
U. The same is true for the determinant of the matrix 5*: the determinant 
dct S {u^ , , u^) is a continuous function in the domain U which is nonzero at 
all points of this domain. A nonzero continuous real function in a connected set 
U cannot take the values of different signs in U. This means that det 5 > or 
det S < 0. This means that the orientation of the triple of vectors forming the 
moving frame of a curvilinear coordinate system is the same for all points of a 
domain where it is defined. Since the space E is equipped with the preferable 
orientation, we can subdivide all curvilinear coordinates in E into right-oriented 
and left-oriented coordinate systems. 

A remark on the smoothness. The definition 2.2 yields the concept of a 
continuously differentiable curvilinear coordinate system. However, the functions 
(2.3) could belong to a higher smoothness class C"\ In this case we say that we 
have a curvilinear coordinate system of the smoothness class C™ . The components 
of the Jacobi matrix (2.2) for such a coordinate system are the functions of the 
class C™~^. Due to the relationship (2.12) the components of the Jacobi matrix 
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(2.9) belong to the same smoothness class C™^^. Hence, the fmictions (2.4) 
belong to the smoothness class C™. 

If we have two curvilinear coordinate systems of the smoothness class C™, then, 
according to the above considerations, the transformation functions (3.5) belong 
to the class C™, while the components of the transition matrices S and T given 
by the formulas (3.7) belong to the smoothness class C™~^. 

§ 4. Vectorial and tensorial fields 
in curvilinear coordinates. 

Let M^, u^, be some curvilinear coordinate system in some domain Z? C E 
and let F be some vector field defined at the points of the domain D. Then 
at a point with coordinates u^, u^, wc have the field vector F(u^,u^,u^) and 
the triple of the frame vectors Ei(w-^, u^, u'^), £2(11-^, u^, u^), ^^{u^ , v? , u^) . Let's 
expand the field vector F in the basis formed by the frame vectors: 

3 

¥{u\u^,u^)^^F\u\u^,u^)-^,{u^,u^,u^). (4.1) 

i=l 

The quantities F^{u^, u^, u'^) in such expansion are naturally called the components 
of the vector field F in the curvilinear coordinates u^, v? , v? . If we have another 
curvilinear coordinate system u^, u^, in the domain D, then we have the other 
expansion of the form (4.1): 

3 

F(u\u^u3) = ^F\u^ ,u'^ ■%{u^ .u"^ ,u^). (4.2) 
1=1 

By means of the formulas (3.6) one can easily derive the relationships binding the 
components of the field F in the expansions (4.1) and (4.2): 

F\n) = j2S]{n)F\n), 

U (4-3) 

u' = u\u^,u^,u^). 

The relationships (4.3) are naturally interpreted as the generalizations for the 
relationships (1.2) from Chapter II for the case of curvilinear coordinates. 

Note that Cartesian coordinate systems can be treated as a special case of 
curvilinear coordinates. The transition functions = u'({t^, -u^, u'^) in the case of 
a pair of Cartesian coordinate systems are linear, therefore the matrix S calculated 
according to the theorem 3.1 in this case is a constant matrix. 

Now let F be either a field of covectors, a field of linear operators, or a field 
of bilinear forms. In any case the components of the field F at some point are 
determined by fixing some basis attached to that point. The vectors of the moving 
frame of a curvilinear coordinate system at a point with coordinates , v?, 
provide the required basis. The components of the field F determined by this 
basis are called the components of the field F in that curvilinear coordinates. The 
transformation rules for the components of the fields listed above under a change 
of curvilinear coordinates generalize the formulas (1.3), (1.4), and (1.5) from 
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Chapter II. For a covectorial field F the transformation rule for its components 
under a change of coordinates looks like 

U (4.4) 



The transformation rule for the components of an operator field F is written as 

p=lq=l (4 



In the case of a field of bilinear (quadratic) forms the generalization of the formula 
(1.5) from Chapter II looks like 

3 3 

p=lg=l (4.6) 
= M*({t\-U^,U^). 

Let F be a tensor field of the type (r, s). In contrast to a vectorial field, 
the value of such a tensorial field at a point have no visual embodiment in 
form of an arrowhead segment. Moreover, in general case there is no visually 
explicit way of finding the numerical values for the components of such a field 
in a given basis. However, according to the definition 1.1 from Chapter II, a 
tensor is a geometric object that for each basis has an array of components 
associated with this basis. Let's denote by F{u^ ,u^,u^) the value of the field F 
at the point with coordinates u^, v?^ v!^ . This is a tensor whose components in 
the basis Ei(u^, m^, u^), E2(u^, u^, u'^), 'Ei^{u^ , v? , v?) arc called components of the 
field F in a given curvilinear coordinate system. The transformation rules for the 
components of a tensor field under a change of a coordinate system follow from 
the formula (1.6) in Chapter II. For a tensorial field of the type (r, s) it looks like 

^1 :::;:(")= E si\{u) . . . s^^;{u)x 

^■■■.1: ^ j.q^^ F^i-::i^iu), (4.7) 



The formula (4.7) has two important differences as compared to the corresponding 
formula (1.7) in Chapter II. In the case of curvilinear coordinates 

(1) the transition functions m*(m^, m^, u^) should not be linear functions; 

(2) the transition matrices ^(u) and r(u) are not necessarily constant matrices. 

Note that these differences do not affect the algebraic operations with tensorial 
fields. The operations of addition, tensor product, contraction, index permutation, 
symmetrization, and alternation are implemented by the same formulas as in 
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Cartesian coordinates. The differences (1) and (2) reveal only in the operation of 
covariant differentiation of tensor fields. 

Any curvilinear coordinate system is naturally equipped with the the metric 
tensor g. This is a tensor whose components are given by mutual scalar products 
of the frame vectors for a given coordinate system: 



The components of the inverse metric tensor g are obtained by inverting the 
matrix g. In a curvilinear coordinates the quantities gij and g^^ are not necessarily 
constants any more. 

We already know that the metric tensor g defines the volume pseudotensor 
w. As before, in curvilinear coordinates its components are given by the formula 
(6.11) from Chapter II. Since the space E has the preferable orientation, the 
volume pseudotensor can be transformed to the volume tensor lj. The formula 
(8.1) from Chapter II for the components of this tensor remains valid in a 
curvilinear coordinate system either. 



Let A be a differentiable tensor field of the type (r, s). In § 5 of Chapter II we 
have defined the concept of covariant differential. The covariant differential VA 
of a field A is a tensorial field of the type (r, s + 1). In an arbitrary Cartesian 
coordinate system the components of the field VA are obtained by differentiating 
the components of the original field A with respect to x^, x'^, and x^. The use of 
curvilinear coordinates does not annul the operation of covariant differentiation. 
However, the procedure of deriving the components of the field VA from the 
components of A in curvilinear coordinates is more complicated. 

Let u^, u^, be some curvilinear coordinate system in a domain Z) C E. 
Let's derive the rule for covariant differentiation of tensor fields in a curvilinear 
coordinate system. We consider a vectorial field A to begin with. This is a field 
whose components are specified by one upper index: {u^ , , u^) . In order to 
calculate the components of the field B = VA we choose some auxiliary Cartesian 
coordinate system . Then we need to do the following maneuver: first 

we transform the components of A from curvilinear coordinates to Cartesian ones, 
then calculate the components of the field B = VA by means of the formula (5.1) 
from Chapter II, and finally, we transform the components of VA from Cartesian 
coordinates back to the original curvilinear coordinates. 

The Cartesian coordinates and the curvilinear coordinates u^, m^, 

are related by the following transition functions: 



5., = (E,(u)|E,(u)). 



(4.8) 



§ 5. Differentiation of tensor fields 
in curvilinear coordinates. 





(5.1) 



The components of the corresponding transition matrices are calculated according 
to the formula (3.7). When applied to (5.1), this formula yields 
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Denote by A''{x^,x'^,x^) the components of the vector field A in the Cartesian 
coordinate system x^, x^, x^. Then we get 

3 

p=i 

For the components of the field B — VA in these Cartesian coordinates, applying 
the formula (5.1) from Chapter II, we get 

p=i 

Now we perform the inverse transformation of the components of B from the 
Cartesian coordinates ^ back to the curvilinear coordinates w^, u^, u^: 

3 3 

V,A^ = B]{n) W ^'(") ^9 • (5.4) 

fe=l g=l 

Let's apply the Leibniz rule for calculating the partial derivative in (5.3). As a 
result we get two sums. Then, substituting these sums into (5.4), we obtain 

5=lp=lVfe=l / 

+EfEE^^w^»^) ^^(")- 

p=l \g=lp=l / 

Note that the matrices S and T are inverse to each other. Therefore, we can 
calculate the sums over k and p in the first summand. Moreover, we replace Tjixi) 
by the derivatives dx'' /du^ due to the formula (5.2), and we get 

^ 3^ ' dxi ^ dui dxi du3 

g=l g=l 

Taking into account all the above arguments, we transform the formula for the 
covariant derivative Vj A* into the following one: 

p=i \fe=i / 

We introduce the special notation for the sum enclosed into the round brackets in 
the above formula, we denote it by r*p: 

T]^[u) = Y^Sl{S.)^y. (5.5) 

k=l 
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Taking into account the notations (5.5), now we can write the rule of covariant 
differentiation of a vector field in curvilinear coordinates as follows: 

The quantities calculated according to (5.5) are called the connection compo- 
nents or the Christoffel symbols. These quantities are some inner characteristics of 
a curvilinear coordinate system. This fact is supported by the following lemma. 

Lemma 5.1. The connection components of a curvilinear coordinate system 
u^, u^, given by the formula (5.5) do not depend on the choice of an auxihary 
Cartesian coordinate system 

Proof. Let's multiply both sides of the equality (5.5) by the frame vector 
and then sum over the index i: 

3 3 3 f)rpk( \ 

^r;.^(u)E,(u) = EE (5.7) 

i=l 1=1 fc=l 

The sum over i in right hand side of the equality (5.7) can be calculated explicitly 
due to the first of the following two formulas: 



gfc = E^feE,, Ep^E^P^'^- (5.8) 

i=l fc=l 

These formulas (5.8) relate the frame vectors Ei, E2, E3 and the basis vectors 
ei, §2, 63 of the auxiliary Cartesian coordinate system. Now (5.7) is written as: 

2=1 fe=l k=l 

The basis vector does not depend on m^, m^, u^. Therefore, it is brought into 
the brackets under the differentiation with respect to . The sum over k in right 
hand side of the above formula is calculated explicitly due to the second formula 
(5.8). As a result the relationship (5.7) is transformed to the following one: 

if = Er}.-E. (5.9) 

i—l 

The formula (5.9) expresses the partial derivatives of the frame vectors back 
through these vectors. It can be understood as another one way for calculating 
the connection components Fjp. This formula comprises nothing related to the 
auxiliary Cartesian coordinates The vector 'Ep{u^ , , u^) is determined 

by the choice of curvilinear coordinates u^, u^, in the domain D. It is sufficient 
to differentiate this vector with respect to and expand the resulting vector in 
the basis of the frame vectors Ei, E2, E3. Then the coefficients of this expansion 
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yield the required values for Fjp. It is obvious that these values do not depend on 
the choice of the auxiliary Cartesian coordinates x^, x'^, above. □ 

Now let's proceed with deriving the rule for covariant differentiation of an 
arbitrary tensor field A of the type (r, s) in curvilinear coordinates. For this 
purpose we need another one expression for the connection components. It is 
derived from (5.5). Let's transform the formula (5.5) as follows: 

fe=l fc=l 

The matrices S and T are inverse to each other. Therefore, upon performing the 
summation over k in the first term we find that it vanishes. Hence, we get 

r;» = -E^.^Xu)^. (5.10) 

fc=l 

Let " be the components of a tensor field A of the type (r, s) in curvilinear 
coordinates. In order to calculate the components of B = VA we do the same 
maneuver as above. First of all we transform the components of A to some 
auxiliary Cartesian coordinate system: 



q\...qs / J v\ ■ ■ ■ Vr qi ' ' ' qs 



wi . . . Ws 



Wl... Wg 



Then we calculate the components of the field B in this auxiliary Cartesian 
coordinate system simply by differentiating: 

f}pi...p-r _ \^ \ "1 • • • 91 •■ 9s ^wi...wj 



Then we perform the inverse transformations of the components of B from the 
Cartesian coordinates back to the original curvilinear coordinate system: 



31---33 + 1 '^Pl 
Pi.-.Pt 
qi... 9s + l 



Qir- T^9l rpqs + 1 p>Pl...Pr 



E 

P1---Pt Vi...Vr 

gi-.-gs + i WX...W, 



E si\ 

■■Pr 

9s + l 

Ertil air rpqi rTn9s + l 

"Spi • ■ • Jp, J-ji ■ ■ ■ 



^ V ^1 ■ ■ ■ Vj. qi ■ ■ ■ qs wi...Ws 

dx'>-+^ 



(5.11) 



Applying the Leibniz rule for differentiating in (5.11), as a result we get three 
groups of summands. The summands of the first group correspond to differenti- 
ating the components of the matrix T, the summands of the second group arise 
when we differentiate the components of the matrix S in (5.11), and finally, the 
unique summand in the third group is produced by differentiating In 
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any one of these summands if the term T^™ or the term S^^ is not differentiated, 
then this term is built into a sum that can be evaluated explicitly: 

3 3 

E nz E ^1: = ^1: ■ 

Therefore, one can evaluate explicitly the most part of the sums in the formula 
(5.11). Moreover, we have the following equality: 



js+1 fl^q, + i 



Taking into accoimt all the above facts, we can bring (5.11) to 



m — lVm—^ \Pm — 1 



"^js+i^t-'-X - E E ( E -^p™ g:;^ 1 



s 3/3 r)Qw„ \ r)4*i - *'- 



n—\ Wji — \ \ Qn—l J 

Due to the formulas (5.5) and (5.10) one can express the sums enclosed into round 
brackets in the above equality through the Christoffcl symbols. Ultimately, the 
formula (5.11) is brought to the following form: 



^ /in---»r „ 3\---}s I 



E\ ^ -pim j\il...V^... ir _ \ ^ \ ^ plUn All-.-ir- 



(5.12) 



m— IDtti — 1 n—lWn — l 



The formula (5.12) is the rule for covariant differentiation of a tcnsorial field A of 
the type (r, s) in an arbitrary curvilinear coordinate system. This formula can be 
commented as follows: the covariant derivative Vj^^^ is obtained from the partial 
derivative d/duj^_^-^ by adding r + s terms — one per each index in the components 
of the field A. The terms associated with the upper indices enter with the positive 
sign, the other terms associated with the lower indices enter with the negative 
sign. In such additional terms each of the upper indices i,„ and each of the lower 
indices j„ arc sequentially moved to the Christoffcl symbol, while in its place we 
write the summation index Vm or Wn- The lower index js+i added as a result 
of covariant differentiation is always written as the first lower index in Christoffcl 
symbols. The position of the summation indices v,„ and Wn in Christoffcl symbols 
is always complementary to their positions in the components of the field A so 
that they always form a pair of upper and lower indices. Though the formula 
(5.12) is rather huge, we hope that due to the above comments one can easily 
remember it and reproduce it in any particular case. 
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§ 6. Transformation of the connection components 
under a change of a coordinate system. 

In deriving the formula for covariant differentiation of tcnsorial fields in curvi- 
linear coordinates we discovered a new type of indexed objects — these are 
Christoffel symbols. The quantities T^j are enumerated by one upper index and 
two lower indices, and their values are determined by the choice of a coordinate 
system. However, they are not the components of a tensorial fields of the type 
(1,2). Indeed, the values of all F^^- in a Cartesian coordinate system are identically 
zero (this follows from the comparison of (5.12) with the formula (5.1) in Chapter 
II). But a tensorial field with purely zero components in some coordinate system 
cannot have nonzero components in any other coordinate system. Therefore, 
Christoffel symbols are the components of a non-tensorial geometric object which 
is called a connection field or simply a connection. 

Theorem 6.1. Let u^, u^, and u^, u^, be two coordinate systems in a 
domain Z? C E. Then the connection components in these two coordinate systems 
are related to each other by means of the foUowing equality: 

= X! X! X! ^/ + X! 'q:[J' ■ 

m— 1 p—1 q—1 m—l 

Here S and T are the transition matrices given by the formulas (3.7). 

A remark on the smoothness. The derivatives of the components of T in 
(6.1) and the formulas (3.7), where the components of T are defined as the partial 
derivatives of the transition functions (3.5), show that the connection components 
can be correctly defined only for coordinate systems of the smoothness class not 
lower than C^. The same conclusion follows from the formula (5.5) for F* . 

Proof. In order to prove the theorem 6.1 we apply the formula (5.9). Let's 
write it for the frame vectors Ei, E2, E3, then apply the formula (3.6) for to 
express Ej through the vectors Ei, E2, and E3: 

k=l m=l 

Applying the Leibniz rule to the right hand side of (6.2), we get two terms: 

k=l m=l q=l 

In the first term in the right hand side of (6.3) we express E„i through the vectors 
El, E2, and E3. In the second term we apply the chain rule and express the 
derivative with respect to through the derivatives with respect to m^, m^, H^: 

2. r., - 2^ 2. ^™ ^ + 2^ 2. a-F • 

k=l k=lm=l 5=1 p=l 



§ 7. CONCORDANCE OF METRIC AND CONNECTION. 



63 



Now let's replace dvP/du^ by Tf relying upon the formulas (3.7) and then apply 
the relationship (5.9) once more in the form of 



9E„ ' 



m — 1 

As a result of the above transformations we can write the equahty (6.3) as fohows: 

3 3 3 p\rpm 3 3 3 

E E . = ;^ ^ ^ E , + 5: 5: 5: Tf f ™ . 

k—\ fc— 1 m— 1 q— 1 p— 1 m— 1 

Now we need only to express E„i through the frame vectors Ei, E2, E3 and collect 
the similar terms in the above formula: 

3 / 3 3 3 3 \ 

E - E EE^^^f ^^/r" - E ^™ ^ . 0. 

k—\ \ m—1 q—1 p—1 m—1 / 

Since the frame vectors Ei, E2, E3 are linearly independent, the expression en- 
closed into round brackets should vanish. As a result we get the equality exactly 
equivalent to the relationship (6.1) that we needed to prove. □ 

§ 7. Concordance of metric and connection. 
Another formula for Christoffel symbols. 

Let's consider the metric tensor g. The covariant differential Vg of the field g 
is equal to zero (sec formulas (6.7) in Chapter II). This is because in any Cartesian 
coordinates x^, x^, in E the components gij of the metric tensor do not depend 
on x^, x'^, x^. In a curvilinear coordinate system the components of the metric 
tensor gij[u^ ,v? ,u^) usually are not constants. However, being equal to zero in 
Cartesian coordinates, the tensor Vg remains zero in any other coordinates: 

Vfc.9y = 0. (7.1) 

The relationship (7.1) is known as the concordance condition for a metric and a 
connection. Taking into account (5.12), we can rewrite this condition as 



3 3 



ErL.9.,-Eri..9- = o. (7.2) 



The formula (7.2) relates the connection components r*"^ and the components 
of the metric tensor gtj. Due to this relationship we can express through 
the components of the metric tensor provided we remember the following very 
important property of the connection components (5.5). 

Theorem 7.1. The connection given by the formula (5.5) is a symmetric con- 
nection, i. e. =Tj^. 

Proof. From (5.2) and (5.5) for T^^ we derive the following expression: 

r^(u)-V5'=^-V^^^ (73) 

9=1 q=l 
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For the functions of the smoothness class the mixed second order partial 
derivatives do not depend on the order of differentiation: 



This fact immediately proves the symmetry of the Christoffel symbols given by 
the formula (7.3). Thus, the proof is over. □ 

Now, returning back to the formula (7.2) relating Tfj and gij, we introduce the 
following notations that simplify the further calculations: 

3 

r^,k^Y.^'^,9kr. (7.4) 

r=l 

It is clear that the quantities Tijk in (7.4) are produced from T^^ by means of 
index lowering procedure described in Chapter II. Therefore, conversely, F^^- are 
obtained from F.yfe according to the following formula: 

r^—E/T,,- (7.5) 

From the symmetry of F*^^ it follows that the quantities Tjjk in (7.4) are also 
symmetric with respect to the indices i and j, i.e. Tijk = ^jik- Using the 
notations (7.4) and the symmetry of the metric tensor, the relationship (7.2) can 
be rewritten in the following way: 

^-n.-n,.=0. (7.6) 

Let's complete (7.6) with two similar relationships applying two cyclic transposi- 
tions of the indices i ^ j ^ k ^ i to the formula (7.6). As a result we obtain 

dgjk 



du'- 
dgkt 



r.yfe - F,fc, = 0, (7.7) 

r^fei — Tjik = 0. 



Let's add the last two relationships (7.7) and subtract the first one from the sum. 
Taking into account the symmetry of Fyfe with respect to i and j, we get 

dgjk _^ dgkz^ _ %j_ _ 2F ■ . = 
du'^ du^ du^ 

Using this equality, one can easily express F^jfe through the components of the 
metric tensor. Then one can substitute this expression into (7.5) and derive 

r—1 ^ ^ 
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The relationship (7.8) is another formula for the Christoffel symbols Tfj, it follows 
from the symmetry of T^j and from the concordance condition for the metric 
and connection. It is different from (5.5) and (5.10). The relationship (7.8) has 
the important advantage as compared to (5.5): one should not use an auxiliary 
Cartesian coordinate system for to apply it. As compared to (5.9), in (7.8) one 
should not deal with vector-functions Ei{u^ , , u^) . All calculations in (7.8) are 
performed within a fixed curvilinear coordinate system provided the components 
of the metric tensor in this coordinate system are known. 




§ 8. Parallel translation. The equation 
of a straight line in curvilinear coordinates. 

Let a be a nonzero vector attached to some point A in the space E. In a 
Euclidean space there is a procedure of parallel translation; applying this pro- 
cedure one can bring the vector a from 
the point A to some other point B. This 
procedure does change neither the modu- 
lus nor the direction of the vector a being 
translated. In a Cartesian coordinate sys- 
tem the procedure of parallel translation 
is described in the most simple way: the 
original vector a at the point A and the 
translated vector a at the point B have 
the equal coordinates. In a curvilinear 
coordinate system the frame vectors at 
Fig. 8.1 the point A and the frame vector at the 

point B form two different bases. There- 
fore, the components of the vector a in the following two expansions 

a = a\A) ■ Ei(A) + a\A) ■ E^iA) + a'{A) ■ Es{A), 
a = a\B) ■ Ei(B) + a\B) ■ E2(B) + a^iB) ■ E,3{B) 

in general case are different. If the points A and B are closed to each other, 
then the triples of vectors Ei(A), E2(A), E3(^) and Ei(S), E2{B), E3(B) are 
approximately the same. Hence, in this case the components of the vector a in 
the expansions (8.1) are slightly different from each other. This consideration 
shows that in curvilinear coordinates the parallel translation should be performed 
gradually: one should first converge the point B with the point A, then slowly 
move the point B toward its ultimate position and record the coordinates of the 
vector a in the second expansion (8.1) at each intermediate position of the point 
B. The most simple way to implement this plan is to link A and B with some 
smooth parametric curve r — r(t), where t S [0, 1]. In a curvilinear coordinate 
system a parametric curve is given by three functions u^(t), u^(i), u^(t) that for 
each t G [0, 1] yield the coordinates of the corresponding point on the curve. 

Theorem 8.1. For a parametric curve given by tlircc functions u^{t), w^(t), 
and u^(t) in some curvilinear coordinate system the components of tlie tangent 
vector T(t) in the moving frame of that coordinate system are determined by the 
derivatives ii^ (t) , 'u?{t), v?{t). 
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Proof. Curvilinear coordinates u^, , determine the position of a point 
in the space by means of the vector- function r ~ r{u^, v?, u^), where r is the 
radius-vector of that point in some auxihary Cartesian coordinate system (see 
formulas (2.4) and (2.5)). Therefore, the vectorial-parametric equation of the 
curve is represented in the following way: 

v = v{u\t),u\t),u^{t)). (8.1) 

Applying the chain rule to the function r(t) in (8.1), we get 

Remember that due to the formula (2.7) the partial derivatives in (8.2) coincide 
with the frame vectors of the curvilinear coordinate system. Therefore the formula 
(8.2) itself can be rewritten as follows: 

T{t) =Y.^^{t) -E.^u^tWitW^t)). (8.3) 

It is easy to see that (8.3) is the expansion of the tangent vector r(t) in the basis 
formed by the frame vectors of the curvilinear coordinate system. The components 
of the vector r(i) in the expansion (8.3) are the derivatives v}{t), u^{t), u^{t). The 
theorem is proved. □ 

Let's apply the procedure of parallel translation to the vector a and translate 
this vector to all points of the curve linking the points A and B (see Fig. 8.1). 
Then we can write the following expansion for this vector 

3 

a = ^a'(0 • E,{u\t),u^{t),u^t)). (8.4) 

i=l 

This expansion is analogous to (8.3). Let's differentiate the relationship (8.4) with 
respect to the parameter t and take into account that a = const: 

Now let's use the formula (5.9) in order to differentiate the frame vectors of the 
curvilinear coordinate system. As a result we derive 

3 / 3 3 \ 

4=1 \ j = l fe=l / 

Since the frame vectors Ei, E2, E3 are linearly independent, we obtain 

«^+EEr;'.:"'«' = o. (8.5) 

j=i k=i 
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The equation (8.5) is called the differential equation of the parallel translation of 
a vector along a curve. This is the system of three linear differential equations of 
the first order with respect to the components of the vector a. Actually, in order 
to perform the parallel translation of a vector a from the point A to the point B 
in curvilinear coordinates one should set the initial data for the components of the 
vector a at the point A (i. e. for t = 0) and then solve the Cauchy problem for the 
equations (8.5). 

The procedure of the parallel translation of vectors along curves leads us to the 
situation where at each point of a curve in E we have some vector attached to 
that point. The same situation arises in considering the vectors x, n, and b that 
form the Frcnet frame of a curve in E (see Chapter I). Generalizing this situation 
one can consider the set of tensors of the type (r, s) attached to the points of some 
curve. Defining such a set of tensors differs from defining a tensorial field in E 
since in order to define a tensor field in E one should attach a tensor to each point 
of the space, not only to the points of a curve. In the case, where the tensors 
of the type (r, s) are defined only at the points of a curve, we say that a tensor 
field of the type (r, s) on a curve is given. In order to write the components of 
such a tensor field A we can use the moving frame Ei, E2, E3 of some curvilinear 
coordinate system in some neighborhood of the curve. These components form a 
set of functions of the scalar parameter t specifying the points of the curve: 

^::X=^::XW- (8-6) 

Under a change of curvilinear coordinate system the quantities (8.6) are transfor- 
med according to the standard rule 



4l::lW= E s;\{t) . . . s;i{t)x 

Pl-.-Pr (8.7) 

x-T]i{t)...T]:{t) AiY::i:{t) 



where S{t) and T{t) are the values of the transition matrices at the points of the 
curve. They are given by the following formulas: 

S{t) = S{u\t),u-'{t),u^{t)), 

T{t)^T{u^{t),u^{t),u^{t)). ■ ' 

We cannot use the formula (5.12) for differentiating the field A on the curve 
since the only argument, which the functions (8.6) depend on, is the parameter t. 
Therefore, we need to modify the formula (5.12) as follows: 



Zl . . . Ij. 



dA 



dt 

r 3 3 s 3 3 (8.9) 

■EE E rX^'^r^:-"-EE E r;X^^4i:::X...... 

m— 1 q—1 — 1 n—1 q—1 w^—l 



The formula (8.9) expresses the rule for covariant differentiation of a tensor field A 
with respect to the parameter t along a parametric curve in curvilinear coordinates 
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u^,u^,u^. Unlike (5.12), the index t beside the nabla sign is not an additional 
index. It is set only for to denote the variable t with respect to which the 
differentiation in the formula (8.9) is performed. 

Theorem 8.2. Under a change of coordinates u^, u^, for other coordinates 
u^, V?, the quantities B]^\"'^j = ^t^}i j"^ calculated by means of the formula 
(8.9) arc transformed according to the rule (8.7) and define a tensor held B ~ VtA 
of the type (r, s) which is called the covariant derivative of the field A with respect 
to the parameter t along a curve. 

Proof. The proof of this theorem is pure calculations. Let's begin with the 
first term in (8.9). Let's express through the components of the field A in 

the other coordinates v} , , v? by means of (8.7). In calculating dA^^^"^-'Jdt this 
is equivalent to differentiating both sides of (8.7) with respect to t: 

qi--- q3 



m=lPl---Pr 

qi--- qs 

s 

I \ ^ \ ^ on qir rpqi rpq„ rpq^ 'ipi-.-Pr 

"T '-'pi ■ ■ ■ '-'pr ^ji ■ ■ ■ ■ ■ ■ ^js ^qi...qs ■ 

n=lPl--Pr 

gi--- qs 

For to calculate the derivatives S"*^ and Tj" in (8.10) we use the fact that the 
transition matrices S and T are inverse to each other: 

3 3 3 3 7 / rrZ„i rpk \ 

~ 2^ 2^ ^v„, '^p,„ - ^-f. ^p^ 

k = l Vm = i "771 = 1 fc=l 

3 3 rfrpk 3 / 3 

•U77,=l fc=l 

3 3 3 / 3 ^rpfe \ 

= E E St - E E ^ T^:- 

k=lWn = l W„ = l \ k=l I 

In order to transform further the above formulas for the derivatives .S^™ and T^*" 
we use the second formula in (8.8): 

3/33 k \ 

^t = -E EE^r^"^p;:' (8-11) 

i;,„ = l \ fc=l q=\ I 
3/33 dT^ \ 

Tl:-Y. EE^r^^M?^^:- (8.12) 

tu„ = l \ k=l q=l / 




Let's substitute (8.11) and (8.12) into (8.10). Then, taking into account the 
relationship (8.7), we can perform the summation over pi, ... , pr and . . . , Qs 
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in the second and the third terms in (8.10) thus transforming (8.10) to 



qi---ql 

r 3 3/3 p.rpk \ 

m=lq=lvm = l \k=l / 

S 3 3 / 3^ Cirpk \ 



n—1 q—1 Wn—1 \ k—1 



JX...Wn---3s 



The second and the third terms in (8.9) and (8.13) are similar in their structure. 
Therefore, one can collect the similar terms upon substituting (8.13) into (8.9). 
Collecting these similar terms, we get the following two expressions 

dui ' ^ dui 

k=l k=l 

as the coefficients. Let's apply (6.1) to these expressions: 

3 flT^ 3 3 3 

•pirn _ \ ^ qim ^7n _ \ ^ \ ^ \ ^ ci„j pr Jip rpk 

qvm, 2^'^k Cj q ^Z^/^/^^r pfc g ^ -Urn J 



k—1 r—1 p—1 k—1 



3 dT^ 3 3 3 

-pWn \ ^ QWn Jri \ ^ \ ^ \ ^ QWn J^T rpp rpk 

^qjr^ Z^'^k "HTT - -^r- ^pk-Lq-l-j^ 



.14) 



dui 

k=l r=lp=lk=l 



If wc take into account (8.14) when substituting (8.13) into (8.9), then the equality 
(8.9) is written in the following form: 



Pl...pr 

qi---qs 



dAPl^P^ 
dt 



+ EEE E EE^:r-i-it^^^-:---- 

m—l q—1 Vm—1 pjn — 1 p—1 k—1 

-EE E EEE^rf^,„^T^r"^^}i:::X...... 

n—1 q—1 Wn—l qn — 1 p—1 k—1 

In order to transform it further we express ^j^''' jj'"' and AJ-_^'\\^^^ j^ through 
the components of the field A in the other coordinate system by means of the 
formula (8.7). Moreover, we take into account that TP ■ii'' upon summing up over 
q yields uP. As a result we obtain: 



31---JS 

Pi...p 
qi...q: 



Eoii air rpqi rpqs 



^ /IPl - Pr 



-EE E fX^'^?J:::;'r-'^- EE E fpg"„^^i?j;::S: 

m—l p—1 Vjj-^ — l n—1 p—1 Wn — 1 



dt 

(8.15) 
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Note that the expression enclosed into round brackets in (8.15) is VtA^i -Pj 
exactly. Therefore, the formula (8.15) means that the components of the field 
VtA on a curve calculated according to the formula (8.9) obey the transformation 
rule (8.7). Thus, the theorem 8.2 is proved. □ 

Now let's return to the formula (8.5). The left hand side of this formula 
coincides with the expression (8.9) for the covariant derivative of the vector field 
a with respect to the parameter t. Therefore, the equation of parallel translation 
can be written as Vta = 0. In this form, the equation of parallel translation can 
be easily generalized for the case of an arbitrary tensor A: 



VtA = 0. 



(8.16) 



The equation (8.16) cannot be derived directly since the procedure of parallel 
translation for arbitrary tensors has no visual representation like Fig. 8.1. 

Let's consider a segment of a straight line given parametrically by the functions 
u^{t), u^it), u^{t) in a curvilinear coordinates. Let t = s he the natural parameter 
on this straight line. Then the tangent vector r(t) is a vector of the unit length at 
all points of the line. Its direction is also unchanged. Therefore, its components li* 
satisfy the equation of parallel translation. Substituting a* = ?i* into (8.5), we get 



3 = 1 k=l 



0. 



.17) 



The equation (8.17) is the differential equation of a straight line in cmvilinear 
coordinates (written for the natural parametrization t ~ s). 

§ 9. Some calculations in polar, 
cylindrical, and spherical coordinates. 

Let's consider the polar coordinate system on a plane. It is given by formulas 
(1.1). Differentiating the expressions (1-1), we find the components of the frame 
vectors for the polar coordinate system: 



El 



cos{ip) 
sm{(p) 



-psm{(p) 
pcos{ip) 



(9.1) 



The column-vectors (9.1) are composed by the coordinates of the vectors Ei and 
E2 in the orthonormal basis. Therefore, we can calculate their scalar products and 
thus find the components of direct and inverse metric tensors g and g: 





1 





9'' = 


1 















p-' 



(9.2) 



Once the components of g and g are known, we can calculate the Christoffel 
symbols. For this purpose we apply the formula (7.8): 



r}i = o. 



r 



11 



0, 



-•^ 12 ^ ^ 21 ^ ^^ 



p2 _ p2 _ -1 



rL = -p. 



22 

-i2 
22 



(9.3) 



0. 
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Let's apply the connection components (9.3) in order to calculate the Laplace 
operator A in polar coordinates. Let '0 be some scalar field: "0 = ip{p, f). Then 



i=i j=i 



k=l 



(9.4) 



The formula (9.4) is a two-dimensional version of the formula (10.15) from 
Chapter II applied to a scalar field. Substituting (9.3) into (9.4), we get 



dp^ p dp p^ dip"^ 



(9.5) 



Now let's consider the cylindrical coordinate system. For the components of 
metric tensors g and g in this case we have 



(9.6) 





1 










1 


















9'' = 





p^' 













1 










1 



From (9.6) by means of (7.8) we derive the connection components: 



^ 11 


= 0, 


^ 12 


= 0, 


^ 21 


= 0, 


^ 13 


= 0, 


^ 31 


= 0, 


^ 22 


= -p, 


^ 23 


= 0, 


^ 32 


= 0, 


^ 33 




T.2 

^ 11 


= 0, 


p2 
^ 12 


= p-\ 


p2 
^ 21 


= p-' 


T.2 

^ 13 


= 0, 


p2 
^ 31 


= 0, 


p2 
22 




T.2 

^ 23 


= 0, 


p2 
^ 32 


= 0, 


p2 
^ 33 


= 0, 


T.3 

^ 11 


= 0, 


T.3 

^ 12 


= 0, 


T.3 

21 


= 0, 


T.3 

^ 13 


= 0, 


p3 
^ 31 


= 0, 


p3 
^ 22 


= 0, 


p3 
^ 23 


= 0, 


p3 
^ 32 


= 0, 


p3 
^ 33 


^ 0. 



(9.7) 



(9.8) 



(9.9) 



Let's rewrite in the dimension 3 the relationship (9.4) for the Laplace operator 
applied to a scalar field ip: 



Substituting (9.7), (9.. 



3 3 



9u' du^ 



k=l 



(9.10) 



and (9.9) into the formula (9.10), we get 



dp"^ p dp p^ dip'^ ^'■^ 



(9.11) 
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Now we derive the formula for the components of rotor in cyUndrical coordina- 
tes. Let A be a vector field and let A^, A^, be its components in cylindrical 
coordinates. In order to calculate the components of the field F = rot A we use 
the formula (10.5) from Chapter II. This formula comprises the volume tensor 
whose components are calculated by formula (8.1) from Chapter II. The sign 
factor in this formula is determined by the orientation of a coordinate system. 
The cylindrical coordinate system can be either right-oriented or left-oriented. It 
depends on the orientation of the auxiliary Cartesian coordinate system 
which is related to the cylindrical coordinates by means of the relationships (1.3). 
For the sake of certainty wc assume that the right-oriented cylindrical coordinates 
are chosen. Then = 1 and for the components of the rotor F = rot A we derive 



3 3 3 3 



= v/d^ E E E E 5"" ^^3k 9=' y,F\ 

i=l j = l k=l q=l 

Taking into account (9.7), (9.8), (9.9), (9.6) and using (9.12), we get 

1 1 dA^ dA^ 
-F = - — p 



(9.12) 



p dif dh 

1 dA^ 1 dA^ 

p dh p dp ' 

dA^ 1 dA^ 



(9.13) 



2 A' 



dp p dtp 

The relationships (9.13) can be written in form of the determinant: 



rot A = 





El 


E2 


E3 


1 


d 


d 


d 


p 


dp 


dip 


dh 




A^ 


p^A^ 


A^ 



(9.14) 



Here Ei, E2, E3 are the frame vectors of the cylindrical coordinates. 

In the case of spherical coordinates, we begin the calculations by deriving the 
formula for the components of the metric tensor g: 



9ij 



1 

p2 . (9.15) 

p^sin^{§) 

Then we calculate the connection components and write then in form of the array: 



r}i = 0, 
0, 



^ 13 



^ 31 



r. 



23 



0, 



r^2 = 0, 



r^i = 0, 

rl, = -psm'{d), 



(9.16) 



§ 9. SOME CALCULATIONS 



73 



r?i = 0, 
r?3 = o> 



r 3 = 0, 



r?2 = p- 



r^2 = 0, 



= p- 



^ 33 ~ 



sin(2i9) 



(9.17) 



p3 
11 


= 0, 


p3 
^ 12 


= 0, 


p3 
21 


-0, 


p3 
^ 13 




p3 
^ 31 


=p~\ 


p3 
^ 22 


= 0, 


p3 
^ 23 


= cot(??), 


p3 
-•^ 32 


^ cot(i9), 


p3 
^ 33 


= 0. 



(9.18) 



Substituting (9.16), (9.17), and (9.18) into the relationship (9.10), we get 



_ _^ 2 d^j ^1 d^tp ^ cot{-&) dip 
dp^ p dp p^ d-d"^ d-d 



1 

p2 sin2(,9) 9^' 



(9.19) 



Let A be a vector field with the components A^, A^, A'^ in the right-oriented 
spherical coordinates. Denote F = rot A. Then for the components of F we get 



^ sin(?9) 



9^3 



1 



d-d sin(t?) dip 
1 dA^ 



- 2 cos(?9) A^ 



p^ sin(i?) dip 
1 dA^ 



. .^.dA^ 2 sin I? ^3 

- sm(i?) — ^ A^, 

dp p 

1 dA^ 2 ,2 



(9.20) 



sin('!9) dp p2 sin(t?) psiii(i^) 
Like (9.13), the formulas (9.20) can be written in form of the determinant: 



rot A = 





El 


E2 


E3 


P^' 


d 


d 


d 


sin(79) 


dp 


dd 


dip 




A^ 


p'A^ 


p2 sin2(i5)A3 



(9.21) 



The formulas (9.5), (9.11), and (9.19) for the Laplace operator and the formulas 
(9.14) and (9.21) for the rotor is the main goal of the calculations performed just 
above in this section. They are often used in applications and can be found in 
some reference books for engineering computations. 

The matrices g in all of the above coordinate systems are diagonal. Such 
coordinate systems are called orthogonal, while the quantities Hi ~ ^fgii are called 
the Lame coefficients of orthogonal coordinates. Note that there is no orthonormal 
curvilinear coordinate system. All such systems are necessarily Cartesian, this fact 
follows from (7.8) and (5.9). 



CHAPTER IV 



GEOMETRY OF SURFACES. 



§ 1. Parametric surfaces. 
Curvilinear coordinates on a surface. 

A surface is a two-dimensional spatially extended geometric object. There 
are several ways for expressing quantitatively (mathematically) this fact of two- 
dimensionality of surfaces. In the three-dimensional Euclidean space E the choice 
of an arbitrary point implies three degrees of freedom: a point is determined by 
three coordinates. In order to decrease this extent of arbitrariness we can bind 
three coordinates of a point by an equation: 

F{x\x'^,x^) = 0. (1.1) 

Then the choice of two coordinates determines the third coordinate of a point. 
This means that we can define a surface by means of an equation in some 
coordinate system (for the sake of simplicity we can choose a Cartesian coordinate 
system). We have already used this method of defining surfaces (see formula (1.2) 
in Chapter I) when we defined a curve as an intersection of two surfaces. 

Another way of defining a surface is the parametric method. Unlike curves, 
surfaces are parameterized by two parameters. Let's denote them and u^: 





x^{u^ 




r = r(u\ u^) = 


x^{u^ 






x^{u^ 





The formula (1.2) expresses the radius-vector of the points of a surface in some 
Cartesian coordinate system as a function of two parameters u^, u^. Usually, only 
a part of a surface is represented in parametric form. Therefore, considering the 
pair of numbers (u^, u^) as a point of M^, we can assume that the point {u^, v?) 
runs over some domain [/ C M^. Let's denote by D the image of the domain U 
under the mapping (1.2). Then D is the domain being mapped^ U is the map or 
the chart, and (1.2) is the chart mapping: it maps U onto D. 

The smoothness class of the surface D is determined by the smoothness class 
of the fimctions x^{u^,v?'), x^(u^,u^), and x^{u^,v?') in formula (1.2). In what 
fallows we shall consider only those surfaces for which these functions are at least 
continuously differentiable. Then, differentiating these functions, we can arrange 
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their derivatives into the Jacobi matrix: 





dx^ 


du^ 


du^ 


dx'^ 


dx'^ 


du^ 


dv? 


dx^ 


dx^ 







(1.3) 



Definition 1.1. A continuously difFerentiable mapping (1.2) is called regular 
at a point (u^^u^) if the rank of the Jacobi matrix (1.3) calculated at that point 
is equal to 2. 

Definition 1.2. A set D is called a regular fragment of a continuously differ- 
entiable surface if there is a mapping u : D ^ U from D to some domain U cM.^ 
and the following conditions arc fulfilled: 

(1) the mapping u : D ^ U is bijective; 

(2) the inverse mapping u^^ : U ~* D given by three continuously difFerentiable 
functions (1.2) is regular at all points of the domain U. 

The Jacobi matrix (1.3) has three minors of the order 2. These are the 
determinants of the following 2x2 matrices: 



dx^ 


dx^ 




dx'^ 


dx'^ 






dx^ 


du^ 


du^ 




du^ 


du^ 




du^ 


du^ 


dx"^ 


dx"^ 




dx^ 


dx^ 




dx^ 


dx^ 


du^ 






du^ 






du^ 





(1.4) 



In the case of regularity of the mapping (1.2) at least one of the determinants (1.4) 
is nonzero. At the expense of renaming the variables x^, x^, a;'^ we always can do 
so that the first determinant will be nonzero: 



9a; ^ 


dx^ 






dx^ 


dx^ 




dv? 



In this case we consider the first two functions x^iiP'^iP') and x^iu^^ii}) in (1.2) 
as a mapping and write them as follows: 

x^ = x^iu^ ^ u^), 
x'^ = X [u ,u''). 

Due to (1.5) the mapping (1.6) is locally invertible. Upon restricting (1.6) to some 
sufficiently small neighborhood of an arbitrary preliminarily chosen point one can 
construct two continuously difFerentiable functions 

u^^u\x\x^), 
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that implement the inverse mapping for (1.6). This fact is well-known, it is a 
version of the theorem on implicit fmictions (see [2], see also the theorem 2.1 
in Chapter III). Let's substitute and from (1.7) into the arguments of the 
third function x^iu^jU^) in the formula (1.2). As a result we obtain the function 
F{x^,x'^) = x^{u^{x'^,x'^),u'^{x'^,x'^)) such that each regular fragment of a surface 
can locally (i. e. in some neighborhood of each its point) be presented as a graph 
of a continuously differentiable function of two variables; 

x^ = F{x\x^). (1.8) 

A remark on singular points. If we give up the regularity condition from 
the definition 1.2, this may cause the appearance of singular points on a surface. 
As an example we consider two surfaces given by smooth functions: 

x^ = (ui)3, r x^ = (wi)3, 

x^ = iu^f, I x^ = iu^f, (1.9) 

X3 = (ul)2 + (u2)2^ [ X^ = (U^)^ + {U^)\ 

In both cases the regularity condition breaks at the point v}- = = 0. As a result 
the first surface (1.9) gains the singularity at the origin. The second surface is 
non-singular despite to the breakage of the regularity condition. 

Marking a regular fragment D on a surface and defining a chart mapping 
:[/—>£) can be treated as introducing a curvilinear coordinate system on 




Fig. 1.1 Fig. 1.2 

the surface. The conditions — const and — const determine two families 
of coordinate lines on the plane of parameters u^, u^. They form the coordinate 
network in U. The mapping (1.2) maps it onto the coordinate network on the 
surface D (see Fig. 1.1 and Fig. 1.2). Let's consider the vectors Ei and E2 tangent 
to the lines of the coordinate network on the surface D: 



(LIO) 
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The formula (1.10) defines a pair of tangent vectors Ei and E2 attached to each 
point of the surface D. 

The vector-function r(it^,u^) which defines the mapping (1.2) can be written in 
form of the expansion in the basis of the auxihary Cartesian coordinate system: 

3 

r{u\u^)="^x''{u\u^)-eg. (1.11) 

q=l 

Substituting the expansion (1.11) into (1.10) we can express the tangent vectors 
El and E2 through the basis vectors ei, 62, 63: 

g=l 

Let's consider the column- vectors composed by the Cartesian coordinates of the 
tangent vectors Ei and E2: 



El 



dx'^ 
dx^ 



E2 



dx^ 
dx^ 



(1.13) 



Note that the column-vectors (1.13) coincide with the columns in the Jacobi 
matrix (1.3). However, from the regularity condition (see the definition 1.1) it 
follows that the column of the Jacobi matrix (1.3) are linearly independent. This 
consideration proves the following proposition. 

Theorem 1.1. The tangent vectors Ei and E2 arc linearly independent at each 
point of a surface. Therefore, they form the frame of the tangent vector fields in D. 

The frame vectors Ei and E2 attached to some point of a surface D define the 
tangent plane at this point. Any vector tangent to the surface at this point lies 
in the tangent plane, it can be expanded in the basis formed by the vectors Ei 
and E2. Let's consider some arbitrary curve 7 lying completely on the surface (see 
Fig. 1.1 and Fig. 1.2). In parametric form such a curve is given by two functions 
of a parameter t. They define the curve as follows: 



= ui(i), 

u2 = u'^{t). 



(1.14) 



By substituting (1.14) into (1.11) or into (1.2) we find the radius-vector of a point 
of the curve in the auxiliary Cartesian coordinate system r{t) = r{u^{t),u'^{t)). 
Let's differentiate r(<) with respect to t and find the tangent vector of the curve 
given by the above two functions (1.14): 



Tit) 



dr 
'dt 



E 



dr 



du' 
~dt 
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Comparing this expression with (1.10), we find that r is expressed as follows: 

2 

T(t) =^7i*-E,. (1.15) 

1=1 

Due to (1.15) the vector r is a hncar combination of the vectors Ei and E2 
forming the tangent frame. Hence, if a curve 7 hes completely on the surface, its 
tangent vector lies 011 the tangent plane to this surface, while the derivatives of 
the functions (1.14) are the components of the vector r expanded in the basis of 
the frame vectors Ei and E2. 

§ 2. Change of curvilinear coordinates on a surface. 

Let's consider two regular fragments Di and D2 on some surface, each equipped 
with with its own curvilinear coordinate system. Assume that their intersection 
D = Di n D2 is not empty. Then in D we have two curvilinear coordinate systems 
u^, and u^, v? . Denote by U and U the preimages of D imder the corresponding 
chart mappings (see Fig. 3.1 in Chapter III). Due to the bijectivity of the chart 
mappings (see definition 1.2) we can construct two mappings 

VLovT^-.U^U, uou^i;J7^[/. (2.1) 

The mappings uou~^ and uou"^ in (2.1) are also bijective, they can be repre- 
sented by the following pairs of functions: 



(2.2) 



Theorem 2.1. Tiie functions (2.2) representing uou ^ and uou ^ are conti- 
nuously difFerentiable. 

Proof. We shall prove the continuous differentiability of the second pair of 
functions (2.2). For the first pair the proof is analogous. Let's choose some point 
on the chart U and map it to D. Then we choose a suitable Cartesian coordinate 
system in E such that the condition (1.5) is fulfilled and in some neighborhood of 
the mapped point there exists the mapping (1.7) inverse for (1.6). The mapping 
(1.7) is continuously differentiable. 

The other curvilinear coordinate system in D induces the other pair of functions 
that plays the same role as the functions (1.6): 



(2.3) 



These are two of three functions that determine the mapping u~^ in form of (1.2). 
The functions = u^{v},v?') and v?' = v?'{v},v?) that determine the mapping 
uou~^ in (2.2) are obtained by substituting (2.3) into the arguments of (1.7): 

u^^u\x\u\u^),x^{U\u'')), ^^^^ 
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The compositions of continuously difFerentiable functions in (2.4) are continuously 
differentiable functions. This fact completes the proof of the theorem. □ 

The functions (2.2), whose continuous differentiability was proved just above, 
perform the transformation or the change of curvilinear coordinates on a surface. 
They are analogous to the functions (3.5) in Chapter III. 

A remark on the smoothness. If the functions (1.2) of both coordinate 
systems u^,u^ and belong to the smoothness class C™, then the transition 

functions (2.2) also belong to the smoothness class C™. 

Let r(u^,u^) and r{u^,u^) be two vector-functions of the form (1.2) for two 
curvilinear coordinate systems in D. They define the mappings and 
acting from the charts U and U to D. Due to the identity = <= (uou^^) 
the function r({t^,{t^) is obtained by substituting the corresponding transition 
functions (2.2) into the arguments of r(it^,u^): 

v{u\u^)=r{u\u\u^),u\u\ii^)). (2.5) 

Let's differentiate (2.5) with respect to and take into account the chain rule 
and the formula (1.10) for the vectors of the tangent frame: 

- dr dr du' du^ 

i—l i—1 

Differentiating the identity r{u^,u'^) = r{u^{u^,u^),iP{u^,u'^)), we derive the 
analogous relationship inverse to the previous one: 

dr dr du'' ^ du^ ~ 

sill = r = ^ J- • r = y r • XLj/i;. 

du^ Qyi Q^i 

fe=i 1=1 

It is clear that the above relationships describe the direct and inverse transitions 
from some tangent frame to another. Let's write them as 

2 2 

Ej=^5]-E„ E, = 5^i;'=.Efe, (2.6) 

i=i fc=i 

where the components of the matrices S and T are given by the formulas 



~2n 9'"' rrk/ 1 2\ '^^^ /r, r.\ 



From (2.7), we see that the transition matrices S and T arc the Jacobi matrices 
for the mappings given by the transition functions (2.2). They are non-degenerate 
and are inverse to each other. 

The transformations (2.2) and the transition matrices S and T related to them 
are used in order to construct the theory of tensors and tensor fields analogous to 
that which we considered in Chapter II and Chapter III. Tensors and tensor fields 
defined through the transformations (2.2) and transition matrices (2.7) are called 
inner tensors and inner tensor fields on a surface: 

^t-X = E ^l---^;^---^ (2-8) 

qi---qs 
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Definition 2.1. An inner tensor of the type {r,s) on a surface is a geometric 
object F whose components in an arbitrary curvihnear coordinate system on that 
surface are enumerated by (r + s) indices and under a change of coordinate system 
are transformed according to the rule (2.8). 

The formula (2.8) differs from the formula (1.6) in Chapter II only in the range 
of indices. Each index here runs over the range of two values 1 and 2. By 
setting the sign factor (—1)'^ = sign(detS') = ±1 into the formula (2.8) we get the 
definition of an inner pseudotensor 

^^v.-;: = E (-1)^ ^1 • • • ^ ^ ■ • ■ ^ (2-9) 

Pl...Pr 

Definition 2.3. An inner pseudotensor of the type (r, s) on a surface is 
a geometric object F whose components in an arbitrary curvilinear coordinate 
system on that surface are enumerated by (r + s) indices and under a change of 
coordinate system are transformed according to the rule (2.9). 

Inner tensorial and pseudotensorial fields are obtained by defining an inner 
tensor or pseudotensor at each point of a surface. The operations of addition, 
tensor product, contraction, transposition of indices, symmetrization and alter- 
nation for such fields are defined in a way similar to that of the case of the fields 
in the space E (see Chapter II). All properties of these operations are preserved. 

A remark on the differentiation. The operation of covariant differentiation 
of tensor fields in the space E was first introduced for Cartesian coordinate 
systems. Then it was extended to the case of curvilinear coordinates. On surfaces, 
as a rule, there is no Cartesian coordinate system at all. Therefore, the operation 
of covariant differentiation for inner tensor fields on a surface should be defined in 
a somewhat different way. 

§ 3. The metric tensor and the area tensor. 

The choice of parameters u^, on a surface determines the tangent frame 
El, E2 on that surface. Let's consider the scalar products of the vectors Ei, E2 
forming the tangent frame of the surface: 

5y = (E.|E,). (3.1) 

They compose the 2x2 Gram matrix g which is symmetric, non-degenerate, and 
positive. Therefore, we have the inequality 

detg>0. (3.2) 

Substituting (2.6) into (3.1), we find that under a change of a coordinate system 
the quantities (3.1) are transformed as the components of an inner tensorial field of 
the type (0, 2). The tensor g with the components (3.1) is called the metric tensor 
of the surface. Note that the components of the metric tensor arc determined by 
means of the scalar product in the outer space E. Therefore, we say that the 
tensor field g is induced by the outer scalar product. For this reason the tensor g 
is called the metric tensor of the induced metric. 
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Symmetric tensors of the type (0, 2) are related to quadratic forms. This fact 
yields another title for the tensor g. It is called the first quadratic form of a 
surface. Sometimes, for the components of the first quadratic form the special 
notations are used: gn — E, gi2 = (?2i = F, (722 = G. These notations are 
especially popular in the earlier publications on the differential geometry: 



(3.3) 



Since the Gram matrix g is non-degenerate, we can define the inverse matrix 
g = g~^. The components of such inverse matrix are denoted by g*-' , setting the 
indices i and j to the upper position: 





E 


F 


9^j = 


F 


G 



3 

E 



g'^ 9jk 



(3.4) 



For the matrix g the proposition analogous to the theorem 6.1 from Chapter II 
is valid. The components of this matrix define an inner tensor field of the type 
(2, 0) on a surface, this field is called the inverse metric tensor or the dual metric 
tensor. The proof of this proposition is completely analogous to the proof of the 
theorem 6.1 in Chapter II. Therefore, here we do not give this proof. 

From the symmetry of the matrix g and from the relationships (3.4) it follows 
that the components of the inverse matrix g are symmetric. The direct and inverse 
metric tensors are used in order to lower and raise indices of tensor fields. These 
operations are defined by the formulas analogous to (9.1) and (9.2) in Chapter II: 



B 



Jl---Js + 1 



-1 277 
~1 J77 + 



■3b + i 



A] 



■ tr + l 
■js-1 



k=l 
2 

E nil... 1777- 1 + qi. 



(3.5) 



k=l 



The only difference of the formulas (3.5) here is that the summation index k runs 
over the range of two numbers 1 and 2. Due to (3.4) the operations of raising and 
lowering indices (3.5) are inverse to each other. 

In order to define the area tensor (or the area pseudotensor) we need the 
following skew-symmetric 2x2 matrix: 






1 


1—' 






(3.6) 



The quantities (3.6) form the two-dimensional analog of the Levi-Civita symbol 
(see formula (6.8) in Chapter II). These quantities satisfy the relationship 



2 2 

p=l q 



(3.7) 



where M is some arbitrary square 2x2 matrix. The formula (3.7) is an analog of 
the formula (6.10) from Chapter II (see proof in [4]). 
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Using the quantities dij and the matrix of the metric tensor g in some curvihn- 
ear coordinate system, we construct the following quantities: 



v/d^dy. (3.8) 



From (3.7) one can derive the following relationship linking the quantities LOij and 
uipq defined according to the formula (3.8) in two different coordinate systems: 



3 3 



= sign(det5) ^^TfT^^ip,. (3.9) 

p=l q=l 



Due to (3.9) the quantities (3.8) define a skew-symmetric inner pseudotcnsorial 
field of the type (0, 2). It is called the area pseudotensor. If on a surface D one of 
the two possible orientations is marked, then the formula 



^^i] = in ydctg dij (3.10) 

defines a tensorial field of the type (0, 2). It is called the area tensor. The formula 
(3.10) differs from (3.8) only in sign factor which is the unitary pseudoscalar 
field defining the orientation (compare with the formula (8.1) in Chapter II). Here 
one should note that not any surface admits some preferable orientation globally. 
The Mobius strip is a well-known example of a non-orientablc surface. 

§ 4. Moving frame of a surface. 
Veingarten's derivational formulas. 

Each choice of a curvilinear coordinate system on a surface determines some 
frame of two tangent vectorial fields Ei, E2 on it. The vectors of such a 
frame define the tangent plane at each point of the surface. However, they are 
insufficient for to expand an arbitrary vector of the space E at that point of the 
surface. Therefore, they are usually completed by a vector that does not belong to 
the tangent plane. 

Definition 4.1. A unit normal vector n to a surface _D at a point A is a 
vector of the unit length attached to the point A and perpendicular to all vectors 
of the tangent plane to D at that point. 

The definition 4.1 fixes the unit normal vector n only up to the sign: at each 
point there are two opposite unit vectors perpendicular to the tangent plane. One 
of the ways to fix n uniquely is due to the vector product: 

^ [El, E2] 

|[Ei,E2]| ^ ' 

The vector n determined by the formula (4.1) depends on the choice of a curvi- 
linear coordinate system. Therefore, under a change of coordinate system it can 
change its direction. Indeed, the relation of the frame vectors Ei,E2 and Ei,E2 
is given by the formula (2.6). Therefore, we write 



Ei=T1-Ei+T2.E2, 



Erril -[-1 I rTi2 TTi 
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83 



Substituting these expressions into the vector product [Ei, E2], we obtain 

[El, E2] = (T^ Ti - Tf Ti) ■ [El, E2] = detr • [Ei, E2]. 

Now we easily derive the transformation rule for the normal vector n: 

n=(-l)^-n. (4.2) 

The sign factor (— 1)"^ = sign(det S) ~ ±1 here is the same as in the formula (2.8). 

Another way of choosing the normal vector is possible if there is a preferable 
orientation on a surface. Suppose that this orientation on D is given by the 
unitary pseudoscalar field ^d- Then n is given by the formula 

In this case the transformation rule for the normal vector simplifies substantially: 

n = n. (4.4) 

Definition 4.2. The tangent frame Ei, E2 of a curvilinear coordinate system 
u^, iL^ on a surface completed by the unit normal vector n is called the moving 
frame or the escort frame of this surface. 

If the normal vector is chosen according to the formula (4.1), the escort frame 
El, E2, n is always right-oriented. Therefore, in this case if we change the 
orientation of the tangent frame Ei, E2, the direction of the normal vector n is 
changed immediately. In the other case, if n is determined by the formula (4.3), 
then its direction does not depend on the choice of the tangent frame Ei, E2. This 
fact means that the choice of the orientation on a surface is equivalent to choosing 
the normal vector independent on the choice tangent vectors Ei,E2. 

There is a special case, when such an independent choice of the normal vector 
does exist. Let D be the boundary of a three-dimensional domain. Then one of 
two opposite normal vectors is the inner normal, the other is the outer normal 
vector. Thus, we conclude that the boundary of a three-dimensional domain in 
the space E is always orientable. 

Let D be some fragment of a surface of the smoothness class . The vectors of 
the escort frame of such a surface are continuously differentiable vector-functions 
of curvilinear coordinates: Ei(m^,u^), E2(m^,u^), and n(u^,M^). The derivatives 
of such vectors are associated with the same point on the surface. Hence, we can 
write the following expansions for them: 

ff=^r^.E, + 6,.n. (4.5) 

fc=i 

The derivatives of the unit vector n are perpendicular to this vector (see lemma 3.1 
in Chapter I). Hence, we have the equality 

S-Eof'B. (4.0) 

fc=l 
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Let's consider the scalar product of (4.5) and the vector n. We also consider the 
scalar product of (4.6) and the vector Ej. Due to (Efe| n) = we get 



(9Ej/auMn) = 6„(n|n) = 6„-, (4.7) 

2 

{E,\dn/du^) = Y,c'ig,,. (4.8) 

k=l 

Let's add the left hand sides of the above formulas (4.7) and (4.8). Upon rather 
easy calculations we find that the sum is equal to zero: 

{dEj/du' I n) + (Ej | dn/du') = a(Ej | n)/du' = 0. 

From this equality we derive the relations of bij and c*' in (4.5) and (4.6): 

2 

k=l 

By means of the matrix of the inverse metric tensor g we can invert this relation- 
ship. Let's introduce the following quite natural notation: 

2 

^f = E^^^-5''- (4-9) 

Then the coefficients in (4.6) can be expressed through the coefficients 6y in 
(4.5) by means of the following formula: 

c- = -bt (4.f0) 
Taking into account (4.10), we can rewrite (4.5) and (4.6) as follows: 

(4.11) 

k=l 

The expansions (4.5) and (4.6) written in form of (4.11) are called the Veingarten's 
derivational formulas. They determine the dynamics of the moving frame and play 
the central role in the theory of surfaces. 

§ 5. Christoffel symbols 
and the second quadratic form. 

Let's study the first Veingarten's derivational formula in two different coordinate 
systems u^, and w^, v? on a surface. In the coordinates w^, v?' it is written as 

ff = Er^ + (5.1) 

fc=i 
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In the other coordinates u^, this formula is rewritten as 

^=^f-.E„. + 5,,.fi. (5.2) 

m— 1 

Let's express the vector E^- in the left hand side of the formula (5.1) through the 
frame vectors of the second coordinate system. For this purpose we use (2.6): 

du' ^ du' ^ du' 9 ^ J Qyt ■ 

q=l q=l q=l 

For the further transformation of the above expression we use the chain rule for 
differentiating the composite function: 

du' ^ du' 2^L^\ i J duP ■ 

m— 1 1 p—l 

The values of the partial derivatives di^q/dvP are determined by the formula (5.2). 
Moreover, we should take into account (2.7) in form of the equality dvP jdu^ = Tf : 

n-p _2- r)T™ 

S - E ^ • + E E E f r,) • e™+ 

rn— 1 9—1 p—l m— 1 



2 2 2 2 

2 2 2 2 2 2 



'r)T" 

+EE(^;^fw-n=EE(idr^^^O-E^ 

(j=l p=l m=l fe=l 



+ E E E E(^' f ™ 5,^„) . E, + E E(^' ^..) • 

g— 1 p—l m— 1 k—1 q—1 p—l 

The unit normal vectors n and n can differ only in sign: n = ±n. Hence, the 
above expansion for dEj/du'^ and the expansion (5.1) both are the expansions in 
the same basis Ei, E2, n. Therefore, we have 

2 no^m 2 2 2 

= E ^™ ^ + E E E ^™ ^1 (5-3) 

m— 1 m— 1 p—l q—1 

h,^±jzjlT!Tfb,,. (5.4) 

p=l q=l 

The formulas (5.3) and (5.4) express the transformation rules for the coefficients 
and hij in Veingarten's derivational formulas under a change of curvilinear 
coordinates on a surface. 

In order to make certain the choice of the sign in (5.4) one should fix some rule 
for choosing the unit normal vector. If we choose n according to the formula (4.1), 
then under a change of coordinate system it obeys the transformation formula 
(4.2). In this case the equality (5.4) is written as 

K=jljl{-^)'T!T]b,,. (5.5) 
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It is easy to see that in this case btj are the components of an inner pseudotensorial 
field of the type (0, 2) on a surface. 

Otherwise, if we use the formula (4.3) for choosing the normal vector n, then n 
does not depend on the choice of a curvilinear coordinate system on a surface (see 
formula (4.4)). In this case 6^ arc transformed as the components of a tensorial 
field of the type (0, 2). The formula (5.4) then takes the form 

p=l q=l 

Tensors of the type (0, 2) correspond to bilinear and quadratic forms. Pseudo- 
tensors of the type (0, 2) have no such interpretation. Despite to this fact the 
quantities bij in Vcingarten's derivational formulas are called the components of 
the second quadratic form b of a surface. The following theorem supports this 
interpretation. 

Theorem 5.1. The quantities F^^ and bij in Veingarten's derivational formulas 
(4.11) arc symmetric with respect to the lower indices i and j. 

Proof. In order to prove the theorem we apply the formula (1.10). Let's write 
this formula in the following way: 

) = d^j • ^^-^^ 

Then let's substitute (5.7) into the first formula (4.11): 

i9^r(u^, w^) 



Y.^'^r'Ek + hj-n. (5.8) 



k=l 

The values of the mixed partial derivatives do not depend on the order of 
differentiation. Therefore, the left hand side of (5.8) is a vector that does not 
change if we transpose indices i and j. Hence, the coefficients T^j and bij of its 
expansion in the basis Ei, E2, n do not change imder the transposition of the 
indices i and j. The theorem is proved. □ 

Sometimes, for the components of the matrix of the second quadratic form the 
notations similar to (3.3) are used: 



L 


M 


M 


N 



(5.9) 



These notations are especially popular in the earlier publications. 

The tensor fields g and b define a pair of quadratic forms at each point of a 
surface. This fact explains in part their titles — the first and the second quadratic 
forms. The first quadratic form is non-degenerate and positive. This situation is 
well-known in linear algebra (see [1]). Two forms, one of which is positive, can be 
brought to the diagonal form simultaneously, the matrix of the positive form being 
brought to the unit matrix. The diagonal elements of the second quadratic form 
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upon such diagonalization are called the invariants of a pair of forms. In order to 
calculate these invariants we consider the following contraction: 



(5.10) 



The quantities b!^ enter the second Veingarten's derivational formula (4.11). They 
define a tensor field (or a pseudotensorial field) of the type (1, 1), i. c. an operator 
field. The operator with the matrix (5.10) is called the Veingarten operator. The 
matrix of this operator is diagonalized simultaneously with the matrices of the 
first and the second quadratic forms, and its eigenvalues are exactly the invariants 
of that pair of forms. Let's denote them by fci and ^2- 

Definition 5.1. The eigenvalues fci(u^,u^) and k2(u^,v?) for the matrix of 
the Veingarten operator are called the principal curvatures of a surface at its point 
with the coordinates , u^. 

From the computational point of view the other two invariants are more 
convenient. These are the following ones: 



H = 



ki + k2 



K = kik2. 



(5.11) 



The invariants (5.11) can be calculated without knowing the eigenvalues of the 
matrix 6f. It is sufficient to find the trace for the matrix of the Veingarten 
operator and the determinant of this matrix: 



H 



2tr(6f), 



K = det{b^). 



(5.12) 



The quantity H in the formulas (5.11) and (5.12) is called the mean curvature, 
while the quantity K is called the Gaussian curvature. There are formulas, 
expressing the invariants H and K through the components of the first and the 
second quadratic forms (3.3) and (5.9): 



H 



1 EN + GL-2 FM 

2 EG-F^ ' 



LN - AP 



Let v{u^,u'^) and w(w^,u^) be the vectors of the basis in which the matrix of 
the first quadratic form is equal to the unit matrix, while the matrix of the second 
quadratic form is a diagonal matrix: 



(5.14) 



Then v and w are the eigenvectors of the Veingarten operator. The vectors (5.14) 
have their three-dimensional realization in the space E: 



El 



w 



El 



w 



E2 



(5.15) 



This is the pair of the unit vectors lying on the tangent plane and being per- 
pendicular to each other. The directions given by the vectors (5.15) are called 
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the principal directions on a surface at the point with coordinates u^, u^. If the 
principal curvatures at this point are not equal to each other: ki 7^ fc2, then the 
principal directions are determined uniquely. Otherwise, if ki = fc2, then any 
two mutually perpendicular directions on the tangent plane can be taken for the 
principal directions. A point of a surface where the principal curvatures are equal 
to each other (fci = fc2) is called an umbilical point. 

A remark on the sign. Remember that depending on the way how we choose 
the normal vector the second quadratic form is either a tensorial field or a 
pseudotensorial field. The same is true for the Veingarten operator. Therefore, 
in general, the principal curvatures ki and k2 are determined only up to the sign. 
The mean curvature H is also determined only up to the sign. As for the Gaussian 
curvature, it has no uncertainty in sign. Moreover, the sign of the Gaussian 
curvature divides the points of a surface into three subsets: for any point of a 
surface if K > 0, the point is called an elliptic point; if K < 0, the point is called 
a hyperbolic point; and finally, if K ^ 0, the point is called a parabolic point. 



§ 6. Covariant differentiation 
of inner tensorial fields of a surface. 

Let's consider the formula (5.3) and compare it with the formula (6.1) in 
Chapter III. These two formulas differ only in the ranges over which the indices 
run. Therefore the quantities T^j, which appear as coefficients in the Veingarten's 
derivational formula, define a geometric object on a surface that is called a 
connection. The connection components are called the Christoffel symbols. 

The main purpose of the Christoffel symbols V^^ is their usage for the covariant 
differentiation of tensor fields. Let's reproduce here the formula (5.12) from 
Chapter III for the covariant derivative modifying it for the two-dimensional case: 



V7 ■ 4«i---''- 



\^ \^ pJm Ail-.-Vr, 



s 2 



(6.1) 



Theorem 6.1. The formula (6.1) correctly defines the covariant differentiation 
of inner tensor fields on a surface that transforms a field of the type (r, s) into a 
field of the type (r, s + 1) if and only if the quantities T'^j obey the transformation 
rule (5.3) under a change of curvilinear coordinates on a surface. 

Proof. Let's begin with proving the necessity. For this purpose we choose 
some arbitrary vector field A and produce the tensor field B = VA of the type 
(1, 1) by means of the formula (6.1). The correctness of the formula (6.1) means 
that the components of the field B are transformed according to the formula (2.8). 
From this condition we should derive the transformation formula (5.3) for the 
quantities F*^ in (6.1). Let's write the formula (2.8) for the field B = VA: 



1^ + E - E E s^^^ f w + E r;-, A . 
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Then we expand the brackets in the right hand side of this relationship. In the 
first summand we replace T[ by dvP /du^ according to the formula (2.7) and we 
express A™ through A' according to the transformation rule for a vector field: 



^ ai" A dvP dA^ 



duP ^ du} duP du^ 

2 (6.2) 

\k=\ / i=i j=i 

Taking into account (6.2), we can cancel the partial derivatives in the previous 
equality and bring it to the following form: 

2 2 2 ciT^TTi 2 2 2 

J — 1 J — 1 m—l m— 1 p— 1 1 

Now we need only to express A'^ through A^ applying the transformation rule for 
the components of a vectorial field and then extract the common factor A^ : 

2 r)7^^" 2 2 2 \ 

Eqk j _ \^ \^ \^ qk rj.p rj.q pm 1 _4J _ n 
g^t 2^ 2^ Z^^rn-Li -Lj ^ pq \ ^ - U. 

m—l m—l p—1 q—1 / 

Since A is an arbitrary vector field, each coefhcicnt enclosed into round brackets 
in the above sum should vanish separately. This condition coincides exactly with 
the transformation rule (5.3). Thus, the necessity is proved. 

Let's proceed with proving the sufficiency. Suppose that the condition (5.3) 
is fulfilled. Let's choose some tensorial field A of the type (r, s) and prove that 
the quantities determined by the formula (6.1) are transformed as 

the components of a tensorial field of the type (r, s + 1). Looking at the formula 
(6.1) we see that it contains the partial derivative dA^^_^"^^J du^"'^^ and other r + s 
terms. Let's denote these terms by A^ '':'' ["] andA^ ' !" [°] . Then 

(6.3) 




EAil...i^ rm] <U...jr rO] 



The tensorial nature of A means that its components are transformed according 
to the formula (2.8). Therefore, in the first term of (6.3) we get: 

Q^il...ir 

h---js _ " I \ ^ oil qir rpqi rpqs API 

+ i + i ^ ^Pl---^P. Ji--- 3. 91... 9. 




Pl---Pr gs + i = l 
qi---qs 



1- a^j.+i 
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m=lPl.--Pr 
131... qs 

~^ '^Pl ■ ■ ■ "^Pr ^Ji ■ ■ ■ QiiJs + l ■ ' ■ ■ 

n=lPl...Pr 

gi... gs 

Here we used the Leibniz rule for differentiating the product of multiple functions 
and the chain rule in order to express the derivatives with respect to u^''+^ through 
the derivatives with respect to . For the further transformation of the above 
formula we replace i9u'^=+i/i9u^''+i by T^'"^^^ according to (2.7) and we use the 
following identities based on the fact that S and T are mutually inverse matrices: 



Pm _ _ ci„, '-'-^k qk 



Qu3. + i ^ ^ + i P" 

9 9 6.4 

ar*" 2 2 ^„g„ 

^ Jr. = 7-t"n Qk "^k 

tu„ = l k=\ 

Upon substituting (6.4) into the preceding equality it is convenient to transpose 

the summation indices: p,„ with and with Wn- Then we get 



Qjji\...i^ 

J'-'' = y s;\... s;^ Ti'... X 

Pl...Pr 

qi... gs 

xf E ^:.^§^-E^[o]+E^ra)' 

\ q^^i — l m— 1 71—1 / 

where the following notations are used for the sake of simplicity: 

^[0]= EE^LSTr^"--."r---^^ 
^ra= EE^lS77<^■■■S:■■■- 

io„=l fc=l 



(6.5) 



(6.6) 



No let's consider " [™] in (6.3). They are obviously defined as follows: 

2 

Aii...ir rm"| _ \ ^ pim Aii...Vm---ii- 

l',„ = l 

Applying the transformation rule (2.8) to the components of the field A, we get: 

2 

Aii-.- iT fm] _ \ ^ \ ^ oil qvm Q^r 
^Ji...J. + iLoJ - 2^ -^Pi • ■ ■ "^Pr^ ■ ■ ■ '='p. ^ 



^,;:::?:"'"=i (6.7) 



rjiqi rpqs p4m 4Pl---Pr 
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For the further transformation of the above expression we use (5.3) written as 

2 fiT'^ 2 2 2 

_ \ " Qim " I \ ' \ " \ " rjiq. + l rpq pfc 

k=l k=lq=lqs + l=i 

Immediately after substituting this expression into (6.7) we perform the cychc 
transposition of the summation indices: r —f pm f- Some sums in the 

resulting expression are evaluated explicitly if we take into account the fact that 
the transition matrices S and T arc inverse to each other: 

Pl...pj. 

Qi ... Os 

/ 2 X (6-8) 

\9s + l = l / 

By means of the analogous calculations one can derive the following formula: 

P1---PT 

Qi ... Qs 

, . (6.9) 

\qs + l = l J 

Now we substitute (6.5), (6.8), and (6.9) into the formula (6.3). Then the entries 
of V^['q] and do cancel each other. A a residue, upon collecting the similar 

terms and cancellations, wc get the formula expressing the transformation rule 
(2.8) applied to the components of the field VA. The theorem 6.1 is proved. □ 

The theorem 6.1 yields a universal mechanism for constructing the covari- 
ant differentiation. It is sufficient to have a connection whose components are 
transformed according to the formula (5.3). We can compare two connections: 
the Euclidean connection in the space E constructed by means of the Cartesian 
coordinates and a connection on a surface whose components are given by the 
Veingarten's derivational formulas. Despite to the different origin of these two 
connections, the covariant derivatives defined by them have many common prop- 
erties. It is convenient to formulate these properties using covariant derivatives 
along vector fields. Let X be a vector field on a surface. For any tensor field A of 
the type (r, s) wc define the tensor field B = VxA of the same type (r, s) given 
by the following formula: 

2 

9=1 

Theorem 6.2. The operation of covariant difFerentiation of tensor Gelds posses- 
ses the foUowing properties: 

(1) Vx(A + B) = VxA + VxB; 

(2) Vx+yA = VxA + VyA; 

(3) V^.xA = e • VxA; 

(4) Vx(A (g) B) = VxA (g) B + A (8) VxB; 

(5) VxC(A) = C(VxA); 
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where A and B are arbitrary differentiable tensor fields, while X and Y are arbi- 
trary vector fields and ^ is an arbitrary scalar field. 

Looking attentively at the theorem 6.2 and at the formula (6.10), we see that 
the theorem 6.2 is a copy of the theorem 5.2 from Chapter II. while the formula 
(6.10) is a two-dimensional analog of the formula (5.5) from the same Chapter II. 
However, the proof there is for the case of the Euclidean connection in the space 
E. Therefore we need to give another proof. 

Proof. Let's choose some arbitrary curvilinear coordinate system on a surface 
and prove the theorem by means of direct calculations in coordinates. Denote 
C = A + B, where A and B are two tensorial fields of the type (r, s). Then 

+^Jl...J.• 
Substituting Cj^'"^"^ into (6.1), upon rather simple calculations we get 

The rest is to multiply both sides of the above equality by X^^+^ and perform 
summation over the index js+i- Applying (6.10), we derive the formula of the 
item (1) in the theorem. 

Note that the quantities -S'J '*'' in the formula (6.10) are obtained as the linear 
combinations of the components of X. The items (2) and (3) of the theorem follow 
immediately from this fact. 

Let's proceed with the item (4). Denote C = A ig) B. Then for the components 
of the tensor field C we have the equality 

Let's substitute the quantities Cj^^.'.'j^ll^ from (6.11) into the formula (6.1) defining 
the covariant derivative. As a result for D = VC we derive 



js + q+l'Vm. jl---jB j 



m— 1 Vjri—^ 

r+p 2 

E\ " ^ii...Vpi„ j^ir + l—V„,...ir+p 

s 2 



n— 1 Wn—1 
s+q 2 



U + q 
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Note that upon collecting the similar terms the above huge formula can be 
transformed to the following shorter one: 



ll... Ir^ 

(6.12) 

Now in order to prove the fourth item of the theorem it is sufficient to multiply 
(6.12) by X^'+i^^ and sum up over the index js+5+1. 

Proceeding with the last fifth item of the theorem, we consider two tensorial 
fields A and B one of which is the contraction of another: 

k=l 

Substituting (6.13) into (6.1), for j^+iBj'^'''^ we obtain 

2 ^ lii . . . ip — i k ip . . . 

^ 3l---3q-lkjg...]s , 



k=l 



m— 1 k—1 Vjn — 1 
s 2 2 



n— 1 k—1 Wn — 1 



The index Vm in (6.14) can be either to the left of the index k or to the right of it. 
The same is true for the index u)„. However, the formula (6.14) does not comprise 
the terms where Vm or w„ replaces the index k. Such terms would have the form: 



tt^L.^^:::i^Vcl' (6-15) 

k=l v=l 
k=l w=l 



It is easy to note that (6.15) and (6.16) differ only in sign. It is sufficient to 
rename fc to w and u; to fc in the formula (6.16). Adding both (6.15) and (6.16) to 
(6.14) would not break the equality. But upon adding them one can rewrite the 
equality (6.14) in the following form: 

k=l 

No in order to complete the proof of the item (5), and thus prove the theorem in 
whole, it is sufficient to multiply the equality (6.17) by X^^^^ and sum up over the 
index js+i. □ 
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Among the five properties of the covariant derivative hsted in the theorem 6.2 
the fourth property written as (6.12) is most often used in calculations. Let's 
rewrite the equality (6.12) as follows: 

The formula (6.18) is produced from (6.12) simply by renaming the indices; 
however, it is more convenient for reception. 

§ 7. Concordance of metric and connection on a surface. 

Earlier we have noted that the covariant differential of the metric tensor g 
in the Euclidean space E is equal to zero (see formula (6.7) in Chapter II). 
This property was called the concordance of the metric and connection. Upon 
passing to curvilinear coordinates we used this property in order to express the 
components of the Euclidean connection through the components of the metric 
tensor (see formula (7.8) in Chapter III). Let's study whether the metric and the 
connection on surfaces are concordant or not. The answer here is also positive. It 
is given by the following theorem. 

Theorem 7.1. The components of the metric tensor g^^ and the connection 
components T^j in arbitrary coordinates on a surface are related by the equality 

^^5- = Is - E 9rj E r'., 3., = (7.1) 

q=l q=l 

which expresses the concordance condition for the metric and connection. 

Proof. Let's consider the first Veingarten's derivational formula in (4.11) and 
let's rewrite it renaming some indices: 

||=ErL-E, + 6..-n. (7.2) 

9=1 

Let's take the scalar products of both sides of (7.2) by Ej and remember that the 
vectors Ej and n are perpendicular. The scalar product of Eq and Ej in the right 
hand side yields the element gqj of the Gram matrix: 

2 

idE./du''\-E,)^J2^l9qj- (7.3) 

9=1 

Now let's transpose the indices i and j in (7.3) and take into account the symmetry 
of the Gram matrix. As a result we obtain 

2 

{E.,\^■E,/^u'')^Y.^l^9^q■ (7.4) 

9=1 
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Then let's add (7.3) with (7.4) and remember the Leibniz rule as applied to the 
differentiation of the scalar product in the space E: 

{d^^/du'' I Ej) + (E, I d^j/du'') ^ 9(E, | Ej)/du'' = 

2 2 

= dg^./du" = J2 9,j + E 9^,■ 

9=1 9=1 

Now it is easy to see that the equality just obtained coincides in essential with 
(7.1). The theorem is proved. □ 

As an immediate consequence of the theorems 7.1 and 5.1 we get the following 
formula for the connection components: 

r— 1 

We do not give its prove here since, in essential, it is the same as in the case of 
the formula (7.8) in Chapter III. 

From the condition V qgij = one can easily derive that the covariant deriva- 
tives of the inverse metric tensor are also equal to zero. For this purpose one 
should apply the formula (3.4). The covariant derivatives of the identical operator 
field with the components 5\ are equal to zero. Indeed, we have 

^"^^ - ^ + E Vk jl - 0. (7.6) 

r— 1 r—1 

Let's differentiate both sides of (3.4) and take into account (7.6): 



V, J2 9'' 9jk - Y.^^i9"' 9jk + g'' y.gjk) = 
\ j=l I J=l 

2 

= T.'^,9''9jk = y,Sl = 0. (7.7) 

In deriving (7.7) we used the items (4) and (5) from the theorem 6.2. The 
procedure of lowering j by means of the contraction with the metric tensor gjk is 
an invertiblc operation. Therefore, (7.7) yields V^g'-' = 0. Now the concordance 
condition for the metric and connection is written as a pair of two relationships 

Vg = 0, Vg = 0, (7.8) 

which look exactly like the relationships (6.7) in Chapter II for the case of metric 
tensors in the space E. 

Another consequence of the theorem 7.1 is that the index raising and the index 
lowering operations performed by means of contractions with the metric tensors 
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g and g commute with the operations of covariant differentiations. This fact is 
presented by the following two formulas: 



\ k=l / k=l I— „N 

/ 2 \ 2 ^^-^^ 



\ fe=i / fc=i 

The relationship (7.9) is easily derived from (7.8) using the items (4) and (5) in 
the theorem 6.2. 

Theorem 7.2. The covariant differential of the area pscudotensor (3.8) on any 
surface is equal to zero: Va; = 0. 

In order to prove this theorem we need two auxiliary propositions which are 
formulated as the following lemmas. 

Lemma 7.1. For any square matrix M whose components are different iable 
functions of some parameter x the equality 

d(lndetM) 



dx 



ir{M-'M') (7.10) 



is fulfilled, where M' is the square matrix composed by the derivatives of the 
corresponding components of the matrix M. 

Lemma 7.2. For any square 2x2 matrix M the equality 

2 

Y{M,qdq.j+Mjqd,q) = trMdy (7.11) 

9=1 

is fulfilled, where dij are the components of the skew- symmetric matrix determined 
by the relationship (3.6). 

The proof of these two lemmas 7.1 and 7.2 as well as the proof of the above 
formula (3.7) from §3 can be found in [4]. 

Let's apply the lemma 7.1 to the matrix of the metric tensor. Let x ^ . Then 
we rewrite the relationship (7.10) as follows: 

1 aydet g _ 1 op dg^p 

VdHt^ du'' 2 ^ ^ ^ Su*^ ■ ^ ' ' 

^ ° q—l p—1 

Note that in (7.11) any array of four numbers enumerated with two indices can 
play the role of the matrix ]\I . Having fixed the index fc, one can use the 
connection components r{.j as such an array. Then we obtain 

2 2 

Y.i^ld,,+n,d^,)=Y.n,d^,. (7.13) 

q=l 9=1 
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Proof for the theorem 7.2. The components of the area pseudotensor uj 
are determmed by the formula (3.8). In order to find the components of the 
pseudotensor Vuj we apply the formula (6.1). It yields 



(^^/det^ 




q=l / 



ydctg ^ 

For the further transformation of this expression we apply (7.12) and (7.13): 

V..., = Vd^(ltt9'''^-t d.,. (7.14) 

\ q=lp=l q=l I 

Now let's express F'.^ through the components of the metric tensor by means of 
the formula (7.5). Taking into account the symmetry of g^', we get 

2 2 2 2 2 

q—1 <?— 1 p— 1 9—1 p— 1 

Substituting this expression into the formula (7.14), we find that it vanishes. 
Hence, VfcWy = 0. The theorem is proved. □ 

A remark on the sign. The area tensor differs from the area pseudotensor 
only by the scalar sign factor Therefore, the proposition of the theorem 7.2 
for the area tensor of an arbitrary surface is also valid. 

A remark on the dimension. For the volume tensor (and for the volume 
pseudotensor) in the Euclidean space E we have the analogous proposition: it 
states that = 0. The proof of this proposition is even more simple than 
the proof of the theorem 7.2. The components of the field a; in any Cartesian 
coordinate system in E are constants. Hence, their derivatives are zero. 

§ 8. Curvature tensor. 

The covariant derivatives in the Euclidean space E are reduced to the partial 
derivatives in any Cartesian coordinates. Is there such a coordinate system for 
covariant derivatives on a surface ? The answer to this question is related to the 
commutators. Let's choose a vector field X and calculate the tensor field Y of the 
type (1,2) with the following components: 

= W^WjX'' - V.X'^. (8.1) 

In order to calculate the components of the field Y we apply the formula (6.1): 

V. V.X'^ . + ± Ff,V,X^ - ± ^VqX\ 

q=l 9=1 
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Let's subtract the second relationship (8.2) from the first one. Then the last terms 



in them do cancel each other due to the symmetry of 



5=1 \ r=l J q=l \ r=l / 

Upon expanding the brackets and some cancellations we get 

= E f # - U + E^r], - r,^,n.)) x^- (8.3) 

r=l \ g=l / 

It is important to note that the formula (8.3) does not contain the derivatives of 
the components of X — they are canceled. Let's denote 

r,k _ jr _ t^J^ ir , ST^ pfc y.q _ ST^ pfc pg a\ 
rij Q^i Q^j ^ Z_^^ iq^ jr Z_^^iq^ir- \°-^) 

g=l q=\ 

The formula (8.3) for the components of the field (8.1) then can be written as 

2 

(V.V, - V, V,)X^ = ^r^, X\ (8.5) 

Let's replace the vector field X by a covector field. Performing the similar 
calculations, in this case we obtain 

2 

(V,V, - V, VOXfc = -Y, Rl^J Xr. (8.6) 

r=l 

The formulas (8.5) and(8.6) can be generalized for the case of an arbitrary tensor 
field X of the type (r, s): 



. r"-' (8.7) 



n—1 Wn — ^ 



Comparing (8.5), (8.6), and (8.7), we see that all of them contain the quantities 



given by the formula (8.4). 



Theorem 8.1. The quantities R'^^j introduced by the formula (8.4) deSne a 
tensor field of the type (1,3). This tensor held is called the curvature tensor or the 
Riemann tensor. 

The theorem 8.1 can be proved directly on the base of the formula (5.3). 
However, we give another proof which is more simple. 



CopyRight © Sharipov R.A., 1996, 2004. 



§ 8. CURVATURE TENSOR. 



99 



Lemma 8.1. Let R be a geometric object which is presented by a four-dimensio- 
nal array R^^j in coordinates. If the contraction of R with an arbitrary vector X 

2 

Y^'=T.^'^J^' (8.8) 

9=1 

is a tensor of the type (1,2), then the object R itself is a tensor of the type (1,3). 

Proof of the lemma. Let u^, and m^, v? be two curvilinear coordinate 
systems on a surface. Let's fix some numeric value of the index r (r = 1 or r = 2). 
Since X is an arbitrary vector, we choose this vector so that its r-th component 
in the coordinate system m^, is equal to unity, while all other components are 
equal to zero. Then for Y^^ in this coordinate system we get 

2 
9=1 

For the components of the vector X in the other coordinate system we derive 

2 

9=1 

then we apply the formula (8.8) on order to calculate the components of the tensor 
Y in the second coordinate system: 

2 2 
m— 1 m=l 

The rest is to relate the quantities Y^'j from (8.9) and the quantities 1^'^ from 
(8.10). From the statement of the theorem we know that these quantities are the 
components of the same tensor in two different coordinate systems. Hence, we get 

n—l p—1 9=1 

Substituting (8.9) and (8.10) into the formula (8.11), we find 

2 2 2 2 
nfe _ qk rpmrpp rpq jjn 

n—l m— 1 p—1 q—1 

This formula exactly coincides with the transformation rule for the components of 
a tensorial field of the type (1, 3) under a change of coordinates. Thus, the lemma 
is proved. □ 

The theorem 8.1 is an immediate consequence of the lemma 8.1. Indeed, the 
left hand side of the formula (8.6) defines a tensor of the type (1, 2) for any choice 
of the vector field X, while the right hand side is the contraction of R and X. 
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The components of the curvature tensor given by the formula (8.4) are enumer- 
ated by three lower indices and one upper index. Upon lowering by means of the 
metric tensor the upper index is usually written in the first position: 

2 

Rqrtj = ^ Rrij 9kq- (8-12) 
k=l 

The tensor of the type (0,4) given by the formula (8.12) is denoted by the same 
letter R. Another tensor is derived from (8.4) by raising the first lower index: 

2 

r=l 

The raised lower index is usually written as the second upper index. The tensors 
of the type (0,4) and (2,2) with the components (8.12) and (8.13) are denoted by 
the same letter R and called the curvature tensors. 

Theorem 8.2. The components of the curvature tensor R determined by the 
connection (7.5) according to the formula (8.4) satisfy the following relationships: 

/I N Tik Ttk 

(2) Rqrij — Rrqij ^ 

(3) Rqrij — Rijqrj 

(4) + + = 0. 

Proof. The first relationship is an immediate consequence of the formula (8.4) 
itself. When transposing the indices i and j the right hand side of (8.4) changes 
the sign. Hence, we get the identity (1) which means that the curvature tensor is 
skew-symmetric with respect to the last pair of its lower indices. 

In order to prove the identity in the item (2) we apply (8.7) to the metric 
tensor. As a result we get the following equality: 

2 

(V, V, - V, V,) = i^^m 9kr + 5,fe) • 
fc=i 

Taking into account (8.12), this equality can be rewritten as 

(V,Vj - VjV,;) ggr = Rrqtj + Rqri] ■ (8.14) 

Remember that due to the concordance of the metric and connection the covariant 
derivatives of the metric tensor are equal to zero (see formula (7.1)). Hence, the 
left hand side of the equality (8.14) is equal to zero, and as a consequence we get 
the identity from the item (2) of the theorem. 

Let's drop for a while the third item of the theorem and prove the fourth item 
by means of the direct calculations on the base of the formula (8.4). Let's write 
the relationship (8.4) and perform twice the cyclic transposition of the indices in 
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it: i j r i. As a result we get the following three equalities: 

q=l q=l 

2 2 
nfc _ rt _ '^^ j» I pfc p9 _ pfc pg 

rlur 2-^ 39 " 2-^ 3' 









du3 























11 



9=1 9=1 
2 2 



nfc _ »j _ I pfc p? _ pfc p 



9 

q=l 9=1 



Let's add all the three above equalities and take into account the symmetry of the 
Christoffer symbols with respect to their lower indices. It is easy to see that the 
sum in the right hand side will be zero. This proves the item (4) of the theorem. 

The third item of the theorem follows from the first, the second, and the third 
items. In the left hand side of the equality that we need to prove we have Rqrij- 
The simultaneous transposition of the indices q ^ r and i ^ j does not change 
this quantity, i. e. we have the equality 

-Rqrij — Rrqji- (8.15) 

This equality follows from the item (1) and the item (2). Let's apply the item (4) 
to the quantities in both sides of the equality (8.15): 



Rqrij Rqijr Rqjri , 

Rrqji Rrjiq Rriqj • 



(8.16) 



Now let's perform the analogous manipulations with the quantity Rijqr- 
To each quantity in (8.17) we apply the item (4) of the theorem: 



Rijqr " Riqrj Rirjq-, 
Rjirq ~ Rjrqi ~^ Rjqir- 



(8.18) 



Let's add the equalities (8.16) and subtract from the sum the equalities (8.18). It 
is easy to verify that due to the items (1) and (2) of the theorem the right hand 
side of the resulting equality is zero. Then, using (8.15) and (8.17), we get 

2 Rqj^ij 2 R^jq^ 0. 

Dividing by 2, now we get the identity that we needed to prove. Thus, the theorem 
is completely proved. □ 

The curvature tensor R given by its components (8.4) has the indices on both 
levels. Therefore, we can consider the contraction: 

2 

Rr] = ^ Rrky (8.19) 
fc=l 
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The formula (8.19) for Rrj can be transformed as follows: 

2 2 

i=l k=l 

From this equality due to the symmetry g^^ and due to the item (4) of the 
theorem 8.2 we derive the symmetry of the tensor Rrj'- 

Rrj = Rjr- (8.20) 

The symmetric tensor of the type (0,2) with the components (8.19) is called the 
Ricci tensor. It is denoted by the same letter R as the curvature tensor. 

Note that there are other two contractions of the curvature tensor. However, 
these contractions do not produce new tensors: 



k=l fc=l 



, Rkrj =0, ^ Rrik = ~Rr 



Using the Ricci tensor, one can construct a scalar field R by means of the formula 

2 2 

i? = ^^i?.,g'-^-. (8.21) 

r=l j = l 

The scalar R{u^,v?) defined by the formula (8.21) is called the scalar curvature of 
a surface at the point with the coordinats The scalar curvature is a result 

of total contraction of the curvature tensor R given by the formula (8.13): 

The formula (8.22) is easily derived from (8.21). Any other ways of contracting 
the curvature tensor do not give other scalars essentially different from (8.21). 

In general, passing from the components of the curvature tensor R^J to the 
scalar curvature, we should lose a substantial part of the information contained in 
the tensor R: this means that we replace 16 quantities by the only one. However, 
due to the theorem 8.2 in two-dimensional case we do not lose the information at 
all. Indeed, due to the theorem 8.2 the components of the curvature tensor R^J 
are skew-symmetric both with respect to upper and lower indices. If /c = r or 
i = j, they do vanish. Therefore, the only nonzero components are RI2, R\\, R^i, 
f?2ii a-iid they satisfy the equalities RI2 = i?2i = ^-^12 = ^^21- Hence, we get 

D C?12 I d21 (-) T>12 

n — + ^21 — ^ ^12- 
Let's consider the tensor D with the following components: 

j^kr ^ R /^k gr _gk gr\ 
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The tensor D is also skew-symmetric with respect to upper and lower indices and 
= i?i2- Hence, these tensors do coincide: D = R. In coordinates we have 

= § (S"^ ' SI) . (8.23) 
By lowering the upper index r, from (8.23) we derive 

Rr., = ^ {S^ 9rj - Or.) ■ (8.24) 

The formula (8.24) determines the components of the curvature tensor on an 
arbitrary surface. For the Ricci tensor this formula yields 



R^l = J9^^ (8.25) 



The Ricci tensor of an arbitrary surface is proportional to the metric tensor. 

A remark. The curvature tensor determined by the symmetric connection 
(7.5) possesses another one (fifth) property expressed by the identity 

VkR%, + + Vji?,\. = 0. (8.26) 

The relationship (8.23) is known as the Bianchi identity. However, in the case of 
surfaces (in the dimension 2) it appears to be a trivial consequence from the item 
(1) of the theorem 8.2. Therefore, we do not give it here. 

§ 9. Gauss equation and Peterson-Codazzi equation. 

Let's consider the Veingarten's derivational formulas (4.11). They can be 
treated as a system of nine vectorial equations with respect to three vector- 
functions Ei(u^,u^), E2(w^,u^), and n{u^,u^). So, the number of the equations is 
greater than the number functions. Such systems are said to be overdetermined. 
Overdetermined systems are somewhat superfiuous. One usually can derive new 
equations of the same or lower order from them. Such equations are called 
differential consequences or compatibility conditions of the original equations. 

As an example we consider the following system of two partial differential 
equations with respect to the function / = f{x,y): 

|^=a(x,y), ^ = b{x,y). (9.1) 

Let's differentiate the first equation (9.1) with respect to y and the second equation 
with respect to x. Then we subtract one from another: 

dO' db 



The equation (9.2) is a compatibility condition for the equations (9.1). It is a 
necessary condition for the existence of the function satisfying the equations (9.1). 
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Similarly, one can derive the compatibility conditions for the system of Vein- 
garten's derivational equations (4.11). Let's write the first of them as 

^=Erj.-E, + 6,.-n. (9.3) 

9=1 

Then we differentiate (9.3) with respect to and express the derivatives dEk/du^ 
and dn/du^ arising therein by means of the derivational formulas (4.11): 



9=1 / 



2 /^p9 2 \ (9.4) 

q=l \ s=l , 



Let's transpose indices i and j in the formula (9.4). The value of the second order 
mixed partial derivative does not depend on the order of differentiation. Therefore, 
the value of the left hand side of (9.4) does not change under the transposition of 
indices i and j . Let's subtract from (9.4) the relationship obtained by transposing 
the indices. As a result we get 

q=l \ ,5 = 1 s=l / 

\ 9=1 9=1 / 

The vectors Ei, E2, and n composing the moving frame are linearly independent. 
Therefore the above equality can be broken into two separate equalities 



s—l s—1 



q=l q=l 



Note that the left hand side of the first of these relationships coincides with the 
formula for the components of the curvature tensor (see (8.4)). Therefore, wc can 
rewrite the first relationship as follows: 

Rl^J^b,,b!~hkl/^. (9.5) 

The second relationship can also be simplified: 

V,b,k = V,6,fc. (9.6) 

It is easy to verify this fact immediately by transforming (9.6) back to the initial 
form applying the formula (6.1). 
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The equations (9.5) and (9.6) are differential consequences of tlie Veingarten's 
derivational formulas (4.11). The first of them is known as the Gauss equation 
and the second one is known as the Peterson-Codazzi equation. 

The tensorial Gauss equation (9.5) contains 16 separate equalities. However, 
due to the relationship (8.24) not all of them are independent. In order to simplify 
(9.5) let's raise the index k in it. As a result we get 

|(5f^-<5J<5f)=6f6^-6^6^ (9.7) 

The expression in right hand side of (9.7) is skew-symmetric both with respect to 
upper and lower pairs of indices and each index in (9.7) runs over only two values. 
Therefore the right hand side of the equation (9.7) can be transformed as 

6f6^-6]6f = S(<5?^^-<Sj5f). (9.8) 

Substituting q = I, k = 2, i = I, j = 2 into (9.8), for B in (9.8) we get 

B = blbl-bl fe? = det(6f ) = i^, 

where K is the Gaussian curvature of a surface (see formula (5.12)). The above 
considerations show that the Gauss equation (9.5) is equivalent to exactly one 
scalar equation which is written as follows: 

R = 2K. (9.9) 

This equation relates the scalar and Gaussian curvatures of a surface. It is also 
called the Gauss equation. 
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CURVES ON SURFACES 



§ 1. Parametric equations of a curve on a surface. 

Let r(t) be the vectorial-parametric equation of a differentiable curve all points 
of which lie on some differentiable surface. Suppose that a fragment D containing 
the points of the curve is charted, i. e. it is equipped with curvilinear coordinates 
u^, V? . This means that there is a bijective mapping u : D ^ U that maps 
the points of the curve to some domain i7 C M^. The curve in the chart U is 
represented not by three, but by two functions of the parameter t: 

(compare with the formulas (1.14) from Chapter IV). The inverse mapping u""'^ is 
represented by the vector-function 

r = r(u\u2). (1.2) 

It determines the radius-vectors of the points of the surface. Therefore, we have 

r{t) = r{u\t),u^{t)). (1.3) 

Differentiating (1.3) with respect to t, we obtain the vector r: 

2 

T(^)=^^^'.E,(^.l(^),^^'(^)). (1.4) 

1=1 

This is the tangent vector of the curve (compare with the formulas (1.15) in 
Chapter IV). The formula (1.4) shows that the vector r lies in the tangent plane 
of the surface. This is the consequence of the fact that the curve in whole lies on 
the surface. 

Under a change of curvilinear coordinates on the surface the derivatives are 
transformed as the components of a tensor of the type (1,0). They determine the 
inner (two-dimensional) representation of the vector r in the chart. The formula 
(1.4) is used to pass from inner to outer (tree-dimensional) representation of this 
vector. Our main goal in this chapter is to describe the geometry of curves lying 
on a surface in terms of its two-dimensional representation in the chart. 

The length integral is an important object in the theory of curves, see formula 
(2.3) in Chapter I. Substituting (1.4) into this formula, we get 
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The expression under integration in (1.5) is the length of the vector t in its inner 
representation. If s = s{t) is the natural parameter of the curve, then, denoting 
du' = dt^ we can write the following formula: 

2 2 

ds^ = ^^9i]du'du^. (1.6) 

1=1 j=i 

The formula (1.6) approves the title «the first quadratic form» for the metric 
tensor. Indeed, the square of the length differential ds^ is a quadratic form of 
differentials of the coordinate functions in the chart. If i = s is the natural 
parameter of the curve, then there is the equality 

2 2 

YT.9^,i^^ = l (1-7) 

1=1 J=l 

that expresses the fact that the length of the tangent vector r of a curve in the 
natural parametrization is equal to unity (see § 2 in Chapter I) . 

§ 2. Geodesic and normal curvatures of a curve. 

Let t = she the natural parameter of a parametric curve given by the equations 
(1.1) in curvilinear coordinates on some surface. Let's differentiate the tangent 
vector t{s) of this curve (1.4) with respect to the parameter s. The derivative of 
the left hand side of (1.4) is given by the Frenet formulas (3.8) from Chapter I: 

2 2 2 o-p 

k ■ neurv = E • + E E • g-T • - (2.1) 

k=l 1=1 j=l 

By ncurv we denote the unit normal vector of the curve in order to distinguish 
it from the unit normal vector n of the surface. For to calculate the derivatives 
dEi/du^ we apply the Veingarten's derivational formulas (4.11): 

2 2 \ / 2 2 \ 

E E r'. ■ Efe + E E • (2-2) 

j=l j=l / \ i=l 3 = 1 j 

Let's denote by fcnorm the coefficient of the vector n in the formula (2.2). This 
quantity is called the normal curvature of a curve: 

2 2 

fcnorm = E E ""^ ' '^2-3) 

1=1 3 = 1 

In contrast to the curvature k, which is always a non-negative quantity, the normal 
curvature of a curve (2.3) can be either positive, or zero, or negative. Taking into 
account (2.3), we can rewrite the relationship (2.2) itself as follows: 

k ■ ncurv - fcnorm ' " = E [ + E E ^i* J ' (2.4) 

k=l \ i=l 3 = 1 / 
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The vector in the right hand side of (2.4) is a hnear combination of the vectors Ei 
and E2 that compose a basis in the tangent plane. Therefore, this vector lies in 
the tangent plane. Its length is called the geodesic curvature of a curve: 



fe=i \ i=i j=i / 



(2.5) 



Due to the formula (2.5) the geodesic curvature of a curve is always non-negative. 
If fcgcod 7^ 0, then, taking into account the relationship (2.5), one can define the 
unit vector n inner and rewrite the formula (2.4) as follows: 



The unit vector n inner 
curve on a surface. 




— ^geod ' 1^ inner- C^-^) 

in the formula (2.6) is called the inner normal vector of a 

Due to (2.6) the vector n inner is a 
linear combination of the vectors ricurv 
and n which are perpendicular to the 
unit vector r lying in the tangent plane. 
Hence, n inner -L T- On the other hand, 
being a linear combination of the vectors 
El and E2, the vector n inner itself lies 
in the tangent plane. Therefore, it is 
determined up to the sign: 



Hi 



= ±k, nl. 



(2.7) 



Fig. 2.1 



Let's move the term k. 



Using the relationship (2.7), sometimes 
one can extend the definition of the vec- 
tor n inner evcn to tliosc poiuts of the 
curve where fcgcod = 0. 
n to the right hand side of the formula (2.6). Then 



this formula is rewritten as follows: 



k ■ THr 



'' geod ' ^ inner 



kr 



n. 



(2., 



The relationship (2.8) can be treated as an expansion of the vector k ■ ricurv as a 
sum of two mutually prpendicular components. Hence, we have 



(^geod)' ^" (^norm) 



(2.9) 



The formula (2.3) determines the value of the normal curvature of a curve in 
the natural parametrization t ~ s. Let's rewrite it as follows: 



2 2 

EE^'>'"' 

2 2 

E E 

i=\ j=i 



(2.10) 
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In the natural parametrization the formula (2.10) coincides with (2.3) because 
of (1.7). When passing to an arbitrary parametrization all derivatives are 
multiplied by the same factor. Indeed, we have 

~di ^ ~ds"dt' ^ ' ^ 

But the right hand side of (2.10) is insensitive to such a change of u\ Therefore, 
(2.10) is a valid formula for the normal curvature in any parametrization. 

The formula (2.10) shows that the normal curvature is a very rough characte- 
ristic of a curve. It is determined only by the direction of its tangent vector r in 
the tangent plane. The components of the matrices gij and bij characterize not 
the curve, but the point of the surface through which this curve passes. 

Let a be some vector tangent to the surface. In curvilinear coordinates w^, 
it is given by two numbers a^, — they are the coefficients in its expansion in 
the basis of two frame vectors Ei and E2. Let's consider the value of the second 
quadratic form of the surface on this vector: 

2 2 

b(a, a) = ^ ^ b,j a' a\ (2.12) 

i=i j=i 

Definition 2.1. The direction given by a nonzero vector a in the tangent 
plane is called an asymptotic direction if the value of the second quadratic form 
(2.12) on such vector is equal to zero. 

Note that asymptotic directions do exist only at those points of a surface where 
the second quadratic form is indefinite in sign or degenerate. In the first case 
the Gaussian curvature is negative: K < 0, in the second case it is equal to zero: 
K = 0. At those points where > there are no asymptotic directions. 

Definition 2.2. A curve on a surface whose tangent vector r at all its points 
lies in an asymptotic direction is called an asymptotic line. 

Due to (2.12) the equation of an asymptotic line has the following form: 

2 2 

h{T,T) = ^^biju'u^ ^0. (2.13) 

i=i j=i 

Comparing (2.13) and (2.10), we see that asymptotic lines are the lines with zero 
normal curvature: fcnorm = 0. On the surfaces with negative Gaussian curvature 
K < at each point there are two asymptotic directions. Therefore, on such 
surfaces always there are two families of asymptotic lines, they compose the 
asymptotic network of such a surface. On any surface of the negative Gaussian 
curvature there exists a curvilinear coordinate system u^, whose coordinate 
network coincides with the asymptotic network of this surface. However, we shall 
not prove this fact here. 

The curvature lines are defined by analogy with the asymptotic lines. These 
are the curves on a surface whose tangent vector lies in a principal direction at 
each point (see formulas (5.14) and (5.15) in § 5 of Chapter IV). The curvature 
lines do exist on any surface, there are no restrictions for the Gaussian curvature 
of a surface in this case. 
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Definition 2.3. A geodesic line on a surface is a curve whose geodesic curva- 
ture is identically equal to zero: fcgcod — 0. 

From the Frenet formula t = k ■ ricurv and from the relationship (2.8) for the 
geodesic lines we derive the following equality: 

dr 

-J-^fcnorm-n. (2.14) 

ds 

In the other words, the derivative of the unit normal vector on a geodesic line is 
directed along the unit normal vector of a surface. This is the external description 
of geodesic lines. The inner description is derived from the formula (2.5): 

2 2 

^' + EEr'>'^' =0- (2.15) 

i=i ]=i 

The equations (2.15) are the differential equations of geodesic lines in natural 
parametrization. One can pass from the natural parametrization to an arbitrary 
one by means of the formula (2.11). 

§ 3. Extremal property of geodesic lines. 

Let's compare the equations of geodesic lines (2.15) with the equations of 
straight lines in curvilinear coordinates (8.17) which we have derived in Chap- 
ter in. These equations have the similar structure. They differ only in the ranges 
of indices: in the case of geodesic lines on a surface they run over two numbers 
instead of three numbers in the case of straight lines. Therefore, geodesic lines 
are natural analogs of the straight lines in the inner geometry of surfaces. The 
following theorem strengthens this analogy. 

Theorem 3.1. A geodesic line connecting two given points A and B on a surface 
has the extremal length in the class of curves connecting these two points. 

It is known that in the Euclidean space E the shortest path from a point A 
to another point B is the segment of straight line connecting these points. The 
theorem 3.1 proclaims a very similar proposition for geodesic lines on a surface. 
Remember that real functions can have local maxima and minima — they are 
called extrema. Apart from maxima and minima, there are also conditional 
extrema (saddle points), for example, the point x = Q for the function y ~ . All 
those points are united by the common property — the linear part of the function 
increment at any such point is equal to zero: /(xq + h) = fix^) + 0{h?). 

In the case of a geodesic line connecting the points A and i? on a surface we 
should slightly deform (variate) this line keeping it to be a line on the surface 
connecting the same two points A and B. The deformed line is not a geodesic line. 
Its length differs from the length of the original line. The condition of the extremal 
length in the theorem 3.1 means that the linear part of the length increment is 
equal to zero. 

Let's specify the method of deforming the curve. For the sake of simplicity 
assume that the points A and -B, and the geodesic line connecting these points in 
whole lie within some charted fragment D of the surface. Then this geodesic line is 
given by two functions (1.1). Let's increment by one the number of arguments in 
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these functions. Then we shaU assume that these functions are sufficiently many 
times differentiable with respect to all their arguments: 

= u^(t, h), 
u = u [t, h). 

For each fixed h in (3.f) we have the functions of the parameter t, they define 
a curve on the surface. Changing the parameter ft,, we deform the curve so 
that in the process of this deformation its point are always on the surface. The 
differentiability of the functions (3.f) guarantees that small deformations of the 
curve correspond to small changes of the parameter h. 

Let's impose to the functions (3.1) a series of restrictions which are easy to 
satisfy. Assume that the length of the initial geodesic line is equal to a and let the 
parameter t run over the segment [0, a]. Let 

u^^XO,/!) = w''(0,0), v!'{a,h) ^u^{a,Q). (3.2) 

The condition (3.2) means that under a change of the parameter h the initial point 
A and the ending point B of the curve do not move. 

For the sake of brevity let's denote the partial derivatives of the functions 
u^{t,h) with respect to t by setting the dot. Then the quantities ?i* = du^/dt 
determine the inner representation of the tangent vector to the curve. 

Assume that the initial line correspond to the value ft = of the parameter ft. 
Assume also that for ft = the parameter t coincides with the natural parameter 
of the geodesic line. Then for ft = the functions (3.1) satisfy the equations (1.7) 
and (2.15) simultaneously. For ft 7^ the parameter t should not coincide with the 
natural parameter on the deformed curve, and the deformed curve itself should 
not be a geodesic line in this case. 

Let's calculate the lengths of the deformed curves. It is the function of the 
parameter ft determined by the length integral of the form (1.5): 



a 

m = / 



'^^'^^ Qij v} ii^ dt. (3.3) 

1=1 j=i 



For ft = we have L(0) = a. The proposition of the theorem 3.1 on the extremity 
of the length now is formulated as L{h) = a + 0{h'^) or, equivalently, as 



dL{h) 



dh 



= 0. (3.4) 

h=0 



Proof of the theorem 3.1. Let's prove the equality (3.4) for the length 
integral (3.3) under the deformations of the curve described just above. Denote by 
A(t, ft) the expression under the square root in the formula (3.3). Then by direct 
differentiation of (3.3) we obtain 

dm_[dX/dh^^^ (3.5) 



dh y 2 VA 
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Let's calculate the derivative in the numerator of the fraction (3.5): 

V i=i j=i / i=i j=i k=i 

k^l i^l j^l i^l j^l k^l 

2 2 2 Q ' ^ 2 2 2 ■ Q ' ^ 

+ E E E SI i'-'' ^ + E E E ^ K . 

k=l i=l j = l fe=l i=l 3 = 1 

Due to the Kronecker symbols 5\. and 5^ in the above expression we can perform 
explicitly the summation over k in the last two terms. Moreover, due to the 
symmetry of gij they are equal to each other: 



i=l j=l k=l 1=1 j=l 



Substituting this expression into (3.5), we get two integrals: 



i=i j=i k=i 



du^ 2y/X dh 



The integral (3.7) contain the second order mixed partial derivatives of (3.1): 

~dh ~ dtdh' 

In order to exclude such derivatives we integrate (3.7) by parts: 



gik dii^ _ gik ii' dv/^ 

Va a/i ^ ^ dh 



a a 

9 ( gik \ du^ 







at V VA / dh 



dt. 



Let's differentiate the equalities (3.2) with respect to h. As a result we find that 
the derivatives du^ /dh vanish at the ends of the integration segment over t. This 
means that non-integral terms in the above formula do vanish. Hence, for the 
integral I2 in (3.7) we obtain 
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Now let's add the integrals Ii and I2 from (3.6) and (3.8). As a result for the 
derivative dL/dh in (3.5) we derive the following equality: 

dL{h) 
dh ^ 

1=1 k=i ^ 

In this equality the only derivatives with respect to the parameter h are du^ /dh. 
For their values at /i = we introduce the following notations: 




(3.9) 



h=i) 



The quantities 5u^ = 5u^{t) in (3.9) are called the variations of the coordinates 
on the initial curve. Note that under a change of curvilinear coordinates these 
quantities arc transformed as the components of a vector (although this fact does 
not matter for proving the theorem). 

Let's substitute h = Q into the above formula for the derivative dL/dh. When 
substituted, the quantity A in the denominators of the fractions becomes equal to 
unity: A(t, 0) = 1. This fact follows from (1.7) since t coincides with the natural 
parameter on the initial geodesic line. Then 



dL{h) 



dh 



h=0 



i=l k 




Since the above equality does not depend on h any more, we replace the partial 
derivative with respect to t by d/dt. All of the further calculations in the right 
hand side are for the geodesic line where t is the natural parameter. 

Let's move the sums over i and k under the integration and let's calculate the 
coefficients of Su'' denoting these coefficients by Uk- 



%y u'-u-J d{giku') 



du'' 2 dt 
1=1 \j=i 



i:i:G&-|f)"'-i:.-- 

i=i j=i ^ ' i=\ 



(3.10) 



Due to the symmetry of 'u? the second term within round brackets in the 
formula (3.10) can be broken into two terms. This yields 

i=l j = l ^ ' i=l 

Let's raise the index k in Uk, i- c. consider the quantities given by the formula 



2 

fe=i 
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For these quantities from the previously derived formula we obtain 

^qk 



^ 2 \du' du3 du 



i=i j=i k=i 



^ 1 u'u'. 



Let's compare this formula with the formula (7.5) in Chapter IV that determines 
the connection components. As a result we get: 



2 2 

1=1 j=i 



(3.11) 



Now it is sufficient to compare (3.11) with the equation of geodesic lines (2.15) 
and derive — 0. The quantities Uk are obtained from If^ by lowering the index: 

2 

Uk = Y.9kq U". 

9=1 

Therefore, the quantities Uk are also equal to zero. From this fact we immediately 
derive the equality (3.4) which means exactly that the extremity condition for the 
geodesic lines is fulfilled. The theorem is proved. □ 

§ 4. Inner parallel translation on a surface. 

The equation of geodesic lines in the Euclidean space E in form of (8.17) was 
derived in Chapter III when considering the parallel translation of vectors in 
curvilinear coordinates. The differential equation of the parallel translation (8.5) 
can be rewritten now in the two-dimensional case: 



3 = 1 k=l 



jk 



0. 



(4.1) 



The equation (4.1) is called the equation of the inner parallel translation of vectors 
along curves on a surface. 

Suppose that we have a surface on some fragment of which the curvilinear 
coordinates u^, v? and a parametric curve (1.1) are given. Let's consider some 
tangent vector a to the surface at the initial point of the curve, i. e. at t = 0. 
The vector a has the inner representation in form of two numbers a^, a^, they are 
its components. Let's set the Cauchy problem for the differential equations (4.1) 
given by the following initial data at t = 0: 



a\t) 



a\t) 



(4.2) 



Solving the Cauchy problem (4.2), we get two functions a^{t) and a^{t) which 
determine the vectors a(t) at all points of the curve. The procedure described 
just above is called the inner parallel translation of the vector a along a curve on 
a surface. 
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Let's consider the inner parallel translation of the vector a from the outer point 
of view, i. e. as a process in outer (three-dimensional) geometry of the space E 
where the surface under consideration is embedded. The relation of inner and 
outer representations of tangent vectors of the surface is given by the formula: 



J2 a' ■ E,. (4.3) 



Let's differentiate the equality (4.3) with respect to t assuming that and 
depend on t as solutions of the differential equations (4.1): 

i—1 i—1 j—1 

The derivatives dEi/du^ are calculated according to Vcingarten's derivational 
formulas (see formulas (4.11) in Chapter IV). Then 



j=i fe=i / \j=i fc=i 



Since the functions a'(i) satisfy the differential equations (4.1), the coefficients at 
the vectors E.^ in this formula do vanish: 

'" = fEEv«^«''^V"- (4.5) 



dt 



\j=i k=i 



The coefficient at the normal vector n in the above formula (4.5) is determined by 
the second quadratic form of the surface. This is the value of the corresponding 
symmetric bilinear form on the pair of vectors a and t. Therefore, the formula 
(4.5) is rewritten in a vectorial form as follows: 

§=b(r,a).n. (4.6) 

The vectorial equation (4.6) is called the outer equation of the inner parallel 
translation on surfaces. 

The operation of parallel translation can be generalized to the case of inner 
tensors of the arbitrary type (r, s) . For this purpose we have introduced the 
operation of covariant differentiation of tensorial function with respect to the 
parameter t on curves (see formula (8.9) in Chapter III). Here is the two- 
dimensional version of this formula: 

(4 7) 

+ EE E rX^^4;:::r-"-EE E r^"^^}l:::X...... 

m— 1 q—1 i^TTi — 1 n— 1 q—1 Wn — ^ 
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In terms of the covariant derivative (4.7) the equation of the inner parallel 
translation for the tcnsorial field A is written as 

VtA = 0. (4.8) 

The consistence of defining the inner parallel translation by means of the equation 
(4.8) follows from the two-dimensional analog of the theorem 8.2 from Chapter III. 

Theorem 4.1. For any inner tensorial function A{t) determined at the points 
of a parametric curve on some surface the quantities Bj^^'Y^ = VtA^-^"*^ calculated 
according to the formula (4.7) define a tensorial function B(t) = VtA of the same 
type {r, s) as the original function A{t). 

The proof of this theorem almost literally coincides with the proof of the 
theorem 8.2 in Chapter III. Therefore, we do not give it here. 

The covariant differentiation Vt defined by the formula (4.7) possesses a series 
of properties similar to those of the covariant differentiation along a vector field 
Vx (see formula (6.10) and theorem 6.2 in Chapter IV). 

Theorem 4.2. The operation of covariant differentiation of tensor-valued func- 
tions with respect to the parameter t along a curve on a surface possesses the 
following properties: 

(1) Vt(A + B) = VtA + VtB; 

(2) Vf (A ® B) = A ® B + A ® VtB; 

(3) VtC(A) = C(VtA). 

Proof. Let's choose some curvilinear coordinate system and prove the theorem 
by means of direct calculations in coordinates. Let C = A + B. Then for the 
components of the tensor- valued function C{t) we have 

+^Jl...J.• 
Substituting C'^*j^ "jJ into (4.7), for the covariant derivative VtC we get 

This equality proves the first item of the theorem. 

Let's proceed with the item (2). Denote C = A B. Then for the components 
of the tensor- valued function C(t) we have 

^il...ir+p _ Aii...i^ -nir+l-.-ir+p (AQ^ 

Let's substitute the quantities Cj^'''^''^'' from (4.9) into the formula (4.8) for the 
covariant derivative. As a result for the components of VjC we derive 
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m— 1 q—1 Vm — 1 

+ E EE^r;:rX"^^;::::::;r;r^^^- 

m— r+l g— 1 v^n — l 
s 2 2 

n— 1 (?— 1 10,1 — 1 

^ E EE^r;:r-^^i?;:::::-r,.,,. 

Note that upon collecting the similar terms the above huge formula can be 
transformed to the following one: 

(4.10) 

x^;::::::;:::+^;i:::;:(v*i^;:::::;;:::)- 

Now it is easy to see that the formula (4.10) proves the second item of the theorem. 

Let's choose two tensor-valued functions A{t) and B(t) one of which is the 
contraction of another. In coordinates this fact looks like 

2 

_ \ " A.il...ip-ikip...ir 

- z^^ji...j,-ifcj,...j.- y^-^^> 

k=l 

Let's substitute (4.11) into the formula (4.7). For Vti?*J;;;]^ we derive 

2 ^j^ii---ip-i kip...ir 

T-T Tyil - ir \ ^ jl--- jq-l jq--- js 

- 2. jt + 

k=l 

+ E E E E (4-12) 

rn — l k—1 q—1 Vm — 1 

s 2 2 2 
\ " \ " \ " \ " pi«„ -a ^il---ip-lkip...i,. 

^ z^l^z^ 2^ ^ji...«,„...fe...j, ■ 

n— 1 k—1 q—1 Wji — l 

In the formula (4.12) the index sequentially occupies the positions to the left 
of the index k and to the right of it. The same is true for the index Wn- However, 
the formula (4.12) has no terms where w„i or Wn replaces the index k. Such terms, 
provided they would be present, according to (4.7), would have the form 

EEEr^"'4::::t^:ri' (4-i3) 

fe=l q=l v=l 

-EEEr^^"^4::::rit;i;;:X- (4.14) 

k=l q=l w=l 



118 



CHAPTER V. CURVES ON SURFACES. 



It is easy to note that (4.13) and (4.14) differ only in sign. Indeed, it is sufficient 
to rename k to v and w to k in the formula (4.14). If we add simultaneously 
(4.13) and (4.14) to (4.12), their contributions cancel each other thus keeping the 
equality valid. Therefore, (4.12) can be written as 



k=l 



The relationship (4.15) proves the third item of the theorem and completes the 
proof in whole. □ 

Under a reparametrization of a curve a new parameter t should be a strictly 
monotonic function of the old parameter t (see details in § 2 of Chapter I) . Under 
such a reparametrization and Vt are related to each other by the formula 



V,A=^.V,A (4.16) 



for any tensor-valued function A on a curve. This relationship is a simple 
consequence from (4.7) and from the chain rule for differentiating a composite 
function. It is an analog of the item (3) in the theorem 6.2 of Chapter IV. 

Let A be a tensor field of the type (r, s) on a surface. This means that at each 
point of the surface some tensor of the type (r, s) is given. If we mark only those 
points of the surface which belong to some curve, we get a tensor-valued function 
A(t) on that curve. In coordinates this is written as 

A^::X(t) = A^::Xiu'(t},v^(t}). (4.17) 

The function A{t) constructed in this way is called the restriction of a tensor field 
A to a curve. The specific feature of the restrictions of tensor fields on curves 
expressed by the formula (4.17) reveals in differentiating them: 

''Jhl_^^y''Jl2±^u'i_ (4.18) 
dt ^ dui ^ ' 

9=1 

Substituting (4.18) into the formula (4.7), we can extract the common factor in 
the sum over q. Upon extracting this common factor we find 



31- 



Y.i^v,A-^::X. (4.19) 



The formula (4.19) means that the covariant derivative of the restriction of a 
tensor field A to a curve is the contraction of the covariant differential VA with 
the tangent vector of the curve. 

Assume that VA = 0. Then due to (4.19) the restriction of the field A to any 
curve is a tensor- valued function satisfying the equation of the parallel translation 
(4.8). The values of such a field A at various points are related to each other by 
parallel translation along any curve connecting these points. 
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Definition 4.1. A tensor field A is called an autoparallel field or a covariantly 
constant field if its covariant differential is equal to zero identically: VA = 0. 

Some of the well-known tensor fields have identically zero covariant differentials: 
this is the metric tensor g, the inverse metric tensor g, and the area tensor 
(pseudotensor) x. The autoparallelism of these fields plays the important role for 
describing the inner parallel translation. 

Let a and b be two tangent vectors of the surface at the initial point of some 
curve. Their scalar product is calculated through their components: 

2 2 

1=1 j=i 

Let's perform the parallel translation of the vectors a and b along the curve 
solving the equation (4.8) and using the components of a and b as initial data in 
Cauchy problems. As a result we get two vector- valued functions a(i) and h{t) on 
the curve. Let's consider the function ^{t) equal to their scalar product: 

2 2 

m = (a(t) |b(t)) = YY.9^Ai) a\t)b>{t). (4.20) 

1=1 j=i 

According to the formula (4.7) the covariant derivative VtV' coincides with the 
regular derivative. Therefore, we have 

1=1 j=i 

Here we used the items (2) and (3) of the theorem 4.2. But Vta* = and 
WtV ~ since we a(i) and b(t) are obtained as a result of parallel translation of 
the vectors a and b. Moreover, Vtgij = due to autoparallelism of the metric 
tensor. For the scalar function ip{t) defined by (4.20) this yields di/j/dt = and 
tp{t) ~ (a I b) = const. As a result of these considerations we have proved the 
following theorem. 

Theorem 4.3. The operation of inner parallel translation of vectors along cur- 
ves preserves the scalar product of vectors. 

Preserving the scalar product, the operation of inner parallel translation pre- 
serves the length of vectors and the angles between them. 

From the autoparallelism of metric tensors g and g we derive the following 
formulas analogous to the formulas (7.9) in Chapter IV: 



\k=l / k=l 

V* I Yd'' A-.:., A = Yd'' A 



(4.21) 



... k 

k=l J k=l 



Then from the formulas (4.21) we derive the following fact. 

Theorem 4.4. The operation of inner parallel translation of tensors commutate 
with the operations of index raising and index lowering. 
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5. Integration on surfaces. Green's formula. 



Let's consider tlie two-dimensional space R^. Let's draw it as a coordinate 
plane u^, u^. Let's choose some simply connected domain Q, outlined by a closed 

piecewise continuously differentiable con- 
tour 7 on the coordinate plane u^, . 
Then we mark the direction (orientation) 
on the contour 7 so that when moving 
in this direction the domain lies to 
the left of the contour. On Fig. 5.1 this 
direction is marked by the arrow. In 
other words, we choose the orientation 
on 7 induced from the orientation of f2, 
i. e. 7 = 9ri. 

Let's consider a pair of continuously 
differentiablc functions on the coordinate 
plane: P{u^,u^) and Q{u^,u'^). Then, 
if all the above conditions are fulfilled, 
there is the following integral identity: 




Fig. 5.1 



Pdu'+Q du^) = 1 1 ^\ du^du\ 



(5.1) 



The identity (5.1) is known as Green's formula (see [2]). The equality (5.1) is an 
equality for a plane. We need its generalization for the case of an arbitrary surface 
in the space E. In such generalization the coordinate plane u^, or some its 
part plays the role of a chart, while the real geometric domain and its boundary 
contour should be placed on a surface. Therefore, the integrals in both parts 
of Green's formula should be transformed so that one can easily write them for 
any curvilinear coordinates on a surface and their values should not depend on a 
particular choice of such coordinate system. 

Let's begin with the integral in the left hand side of (5.1). Such integrals are 
called path integrals of the second kind. Let's rename P to vi and Q to V2. Then 
the integral in the left hand side of (5.1) is written as 



(5.2) 



In order to calculate the integral (5.2) practically the contour 7 should be 
parametrized, i.e. it should be represented as a parametric curve (1.1). Then 
the value of an integral of the second kind is calculated as follows: 




(5.3) 



This formula reducing the integral of the second kind to the regular integral over 
the segment [a, h] on the real axis can be taken for the definition of the integral 
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(5.2). The sign is chosen regarding to the direction of the contour on Fig. 5.1. If 
a < b and if when t changes from a to b the corresponding point on the contour 
moves along the arrow, we choose plus in (5.3). Otherwise, we choose minus. 
Changing the variable t = Lp{t) in the integral (5.3) and choosing the proper 
sign upon reparametrization of the contour, one can verify that the value of this 
integral does not depend on the choice of the parametrization on the contour. 

Now let's change the curvilinear coordinate system on the surface. The 
derivatives in the integral (5.3) under a change of curvilinear coordinates on the 
surface are transformed as follows: 

- du' ^ du^ j i 

—r- = > S]u^. (5.4) 
dt ^ du^ dt ^ ^ ^ ' 

Substituting (5.4) into the formula (5.3), for the integral / we derive: 

(5.5) 



Now let's write the relationship (5.3) in coordinates v}, iP. For this purpose we 
rename to and Vi to Vi in the formula (5.3): 

/ = ± J (Y.v^u'^ dt. (5.6) 

Comparing the formulas (5.5) and (5.6), we see that these formulas are similar 
in their structure. For the numeric values of the integrals (5.3) and (5.6) to be 
always equal (irrespective to the form of the contour 7 and its parametrization) 
the quantities Vi and Vi should be related as follows: 

2 2 

i=l i=l 

These formulas represent the transformation rule for the components of a covec- 
torial field. Thus, we conclude that any path integral of the second kind on a 
surface (5.2) is given by some inner covectorial field on this surface. 

Now let's proceed with the integral in the right hand side of the Green's formula 
(5.1). Distracting for a while from the particular integral in this formula, let's 
consider the following double integral: 

I = F du^du^. (5.7) 



A change of curvilinear coordinates can be interpreted as a change of variables in 
the integral (5.7). Remember that a change of variables in a multiple integral is 
performed according to the following formula (see [2]): 

Fdu^du^^ F\dctJ\du^dii^, (5.8) 
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where J is the Jacobi matrix determined by the change of variables: 



J 









dv? 




dv^ 




dv? 



(5.9) 



The Jacobi matrix (5.9) coincides with the transition matrix S (see formula (2.7) 
in Chapter IV). Therefore, the function F being integrated in the formula (5.7) 
should obey the transformation rule 



F = I det S*! 



(5.10) 



under a change of curvilinear coordinates on the surface. The quantity F has 
no indices. However, due to (5.10), this quantity is neither a scalar nor a 
pseudoscalar. In order to change this not very pleasant situation the integral (5.7) 
over a two-dimensional domain on a surface is usually written as 



\J det g / du^dv?, 



(5.11) 



where det g is the determinant of the first quadratic form. In this case the 
quantity / in the formula (5.11) is a scalar. This fact follows from the equality 
detg = (detT)^ detg that represent the transformation rule for the determinant 
of the metric tensor under a change of coordinate system. 

Returning back to the integral in the right hand side of (5.1), we transform it 
to the form (5.11). For this purpose we use the above notations P = ui, Q = U2, 
and remember that vi and V2 are the components of the covectorial field. Then 



dQ dP dv2 dvi 
du^ dv^ du^ du^ 



(5.12) 



The right hand side of (5.12) can be represented in form of the contraction with 
the unit skew-symmetric matrix d^^ (see formula (3.6) in Chapter IV): 



(~v r\ '2. '2, i~, 2, r\ / 2. 

i=i j=i i=i 



dv} d 



(5.13) 



Note that the quantities dij with lower indices enter the formula for the area 
tensor lj (see (3.7) in Chapter IV). Let's raise the indices of the area tensor by 
means of the inverse metric tensor: 



2 2 2 2 

J2 E 9" = E E 9'" 9" dp 

p—l q—l p—1 q—1 



Applying the formula (3.7) from Chapter IV, we can calculate the components of 
the area tensor a;*-' in the explicit form: 



Vdetg-id*^'. 



(5.14) 
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The formula (5.14) expresses w*-' through d*-'. Now we use (5.14) in order to 
express d*-' m the formula (5.13) back through the components of the area tensor: 

4=1 \J = 1 / 

In order to simplify the further calculations we denote 

2 

y' = ^t^'^\;,. (5.15) 

i=i 

Taking into account (5.15), the formula (5.13) can be written as follows: 



dv2 dvi ^ d{y/dctg 



du^ du^ ^ dv} 

(5.16) 

dy^ l9(lndetg) \ 



du^ 2 du'- 



The logarithmic derivative for the determinant of the metric tensor is calculated 
by means of the lemma 7.1 from Chapter IV. However, we need not repeat these 
calculations here, since this derivative is already calculated (see (7.12) and the 
proof of the theorem 7.2 in Chapter IV): 

With the use of (5.17) the formula (5.16) is transformed as follows: 

i=l \ 9=1 / 

In this formula one easily recognizes the contraction of the covariant differential of 
the vector field y. Indeed, we have 

i—l 

Using the formula (5.18), the notations (5.15), and the autoparallclism condition 
for the area tensor WqUj'^^ ~ 0, we can write the Green's formula as 

j^Vidu"^ ^ j I ^^w'^V.Wj yjdei^du^du^. (5.19) 



The sign factor in (5.19) should be especially commented. The condition that 
the domain O should lie to the left of the contour 7 when moving along the arrow 
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is not invariant under an arbitrary change of coordinates u^, li} by m^, m^. Indeed, 
if we set v} — ~v} and — u^, we would have the mirror image of the domain 
n and the contour 7 shown on Fig. 5.1. This means that the direction should be 
assigned to the geometric contour 7 lying on the surface, not to its image in a 
chart. Then the sign factor (^d in (5.19) can be omitted. 

The choice of the direction on a geometric contour outlining a domain on a 
surface is closely related to the choice of the normal vector on that surface. The 
normal vector n should be chosen so that when observing from the end of the 
vector n and moving in the direction of the arrow along the contour 7 the domain 
il should lie to the left of the contour. The choice of the normal vector n defines 
the orientation of the surface thus defining the unit pseudoscalar field ^^i. 

§ 6. Gauss-Bonnet theorem. 

Let's consider again the process of inner parallel translation of tangent vectors 
along curves on surfaces. The equation (4.6) shows that from the outer (three- 
dimensional) point of view this parallel translation differs substantially from the 
regular parallel translation: the vectors being translated do not remain parallel 
to the fixed direction in the space — they change. However, their lengths are 
preserved, and, if we translate several vectors along the same curve, the angles 
between vectors are preserved (see theorem 4.3). 

From the above description, we see that in the process of parallel translation, 
apart from the motion of the attachment point along the curve, the rotation of 
the vectors about the normal vector n occurs. Therefore, we have the natural 
problem — how to measure the angle of this rotation ? We consider this problem 
just below. 

Suppose that we have a surface equipped with the orientation. This means 
that the orientation field and the area tensor lj are defined (see formula (3.10) 
in Chapter IV). We already know that fixes one of the two possible normal 
vectors n at each point of the surface (see formula (4.3) in Chapter IV). 

Theorem 6.1. The inner tensor field Q of the type (1, 1) with the components 

2 

0] = Y,co,kg'' (6.1) 

k=l 

is an operator field describing the counterclockwise rotation in the tangent plane 
to the angle 7r/2 = 90° about the normal vector n. 

Proof. Let a be a tangent vector to the surface and let n be the unit normal 
vector at the point where a is attached. Then, in order to construct the vector 
b = 0(a) obtained by rotating a counterclockwise to the angle tt/2 ~ 90° about 
the vector n one can use the following vector product: 

b = e(a) = [n, a]. (6.2) 

Let's substitute the expression given by the formula (4.3) from Chapter IV for the 
vector n into (6.2). Then let's expand the vector a in the basis Ei, E2: 



a = • El + • E2. 



(6.3) 
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As a result for the vector b in the formula (6.2) we derive 

In order to calculate the denominator in the formula (6.4) we use the well-known 
formula from the analytical geometry (see [4]): 



[El, E2II' = det 



(El I El) (E1IE2) 
(E2 I El) (E2 I E2) 



detg. 



As for the numerator in the formula (6.4), here we use the not less known formula 
for the double vectorial product: 

[[El, E2], Ej] = E2 • (El I Ej) - El • (E, I E2). 

Taking into account these two formulas, we can write (6.4) as follows: 

gi,--E2__g2,.Ei ^^^. 
^ Vdetg 

Using the components of the area tensor (5.14), no we can rewrite (6.5) in a more 
compact and substantially more elegant form: 

i=l \ j = l k=l / 

From this formula it is easy to extract the formula (6.1) for the components of the 
linear operator O relating b and a. The theorem is proved. □ 

The operator field 8 is the contraction of the tensor product of two fields uj 
and g. The autoparallelism of the latter ones means that 8 is also an autoparallcl 
field, i.e. V8 = 0. 

We use the autoparallelism of 8 in the following way. Let's choose some 
parametric curve 7 on a surface and perform the parallel translation of some unit 
vector a along this curve. As a result we get the vector- valued function a.{t) on the 
curve satisfying the equation of parallel translation Vja = (see formula (4.8)). 
Then we define the vector-function h{t) on the curve as follows: 

b(<) - 8(a(t)). (6.6) 

From (6.6) we derive Vt(b) = Vt8(a) + 8(Vta) = 0. This means that the 
function (6.6) also satisfies the equation of parallel translation. It follows from 
the autoparallelism of 8 and from the items (2) and (3) in the theorem 4.2. The 
vector- functions a(t) and h(t) determine two mutually perpendicular unit vectors 
at each point of the curve. There are the following obvious relationships for them: 

8(a) = b, 8(b) = -a. (6.7) 
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Let's remember for the further use that a{t) and h{t) are obtained by parallel 
translation of the vectors a(0) and b(0) along the curve from its initial point. 

Now let's consider some inner vector field x on the surface (it is tangent to the 
surface in the outer representation). If the field vectors x(w^,u^) are nonzero at 
each point of the surface, they can be normalized to the unit length: x x/|x|. 
Therefore, we shall assume x to be a field of unit vectors: |x| = 1. At the points 
of the curve 7 this field can be expanded in the basis of the vectors a and b: 

X = cos{(p) ■ a + sin((p) • b. (6.8) 

The function ip{t) determines the angle between the vector a and the field vector x 
measured from a to x in the counterclockwise direction. The change of ip describes 
the rotation of the vectors during their parallel translation along the curve. 

Let's apply the covariant differentiation Vf to the relationship (6.8) and take 
into account that both vectors a and b satisfy the equation of parallel translation: 

Vtx = (— sin((^) • a + cos((^) • b) • (6.9) 

Here we used the fact that the covariant derivative Vt for the scalar coincides with 
the regular derivative with respect to t. In particular, we have Vtf ~ (p. Now we 
apply the operator 8 to both sides of (6.8) and take into account (6.7): 

8(x) = cos(i^) • b — sin((y9) • a. (6.10) 

Now we calculate the scalar product of 0(x) from (6.10) and VtX from (6.9). 
Remembering that a and b are two mutually perpendicular unit vectors, we get 

(e(x) I Vfx) = (cos2((^) + sin^iip)) ip^ip. (6.11) 

Let's write the equality (6.11) in coordinate form. The vector-function x(i) on the 
curve is the restriction of the vector field x, therefore, the covariant derivative VtX 
is the contraction of the covariant differential Vx with the tangent vector of the 
curve (see formula (4.19)). Hence, we have 

222 

q=l i=l 3=1 

Here in deriving (6.12) from (6.11) we used the formula (6.1) for the components 
of the operator field 8. 

Let's discuss the role of the field x in the construction described just above. 
The vector field x is chosen as a reference mark relative to which the rotation 
angle of the vector a is measured. This way of measuring the angle is relative. 
Changing the field x, we would change the value of the angle tp. We have to admit 
this inevitable fact since tangent planes to the surface at different points are not 
parallel to each other and we have no preferable direction relative to which we 
could measure the angles on all of them. 

There is a case where we can exclude the above uncertainty of the angle. Let's 
consider a closed parametric contour 7 on the surface. Let [0, 1] be the range over 

CopyRight © Sliaripov R.A., 1996, 2004. 
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which the parameter t runs on such contour. Then x(0) and x(l) do coincide. 
They represent the same field vector at the point with coordinates u^(0), m^(0): 

x(0) x(l) = x{u\0),u^{0)). 

Unhke x(t), the function a{t) is not the restriction of a vector field to a curve 7. 
Therefore, the vectors a(0) and a(l) can be different. This is an important feature 
of the inner parallel translation that differs it from the parallel translation in the 
Euclidean space E. 

In the case of a closed contour 7 the difference (p{l) — (p{0) characterizes the 
angle to which the vector a turns a as a result of parallel translation along the 
contour. Note that measuring the angle from x to a is opposite to measuring it 
from a to X in the formula (6.8). Therefore, taking for positive the angle measured 
from X in the counterclockwise direction, we should take for the increment of the 
angle gained during the parallel translation along 7 the following quantity: 

1 

Atp = (p{0) - ip{l) = - / (pdt. 



Let's substitute (6.12) for (p into this formula. As a result we get 



1 ' 2 2 2 



Aip^- ^^^(a;*t^y V,a;J)u« Ut. (6.13) 

\q=l 1=1 J=l / 







Comparing (6.13) with (5.3), we see that (6.13) now can be written in the form of 
a path integral of the second kind: 



Aip = 



« — 1 — 1 ^ — 1 



^ q=l i=l j=l 



Assume that the contour 7 outlines some connected and simply connected 
fragment fl on the surface. Then for this fragment O we can apply to (6.14) the 
Green's formula written in the form of (5.19): 

» » 2 2 2 2 

A^=-^D / / EEEE'^'' W^ixPujpgVjX") ^"d^du^du^. 
{/ i=li=lp=lg=l 

If the direction of the contour is in agreement with the orientation of the surface, 
then the sign factor can be omitted: 

2222 

= - // E E E E ^''^f ^ + 

O i=lj=lp=lq=l (6.15) 



ViX^ uj'-^ujpqVjx'^) ^/detgdu^du^. 
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Let's show that the term WiX^ uj^^ ujpqW jx"^ in (6.15) yields zero contribution to 
the value of the integral. This feature is specific to the two-dimensional case where 
we have the following relationship: 

u'^LUpg = d'' dp, =^Sl,Sl~ SI (5^. (6.16) 

The proof of the formula (6.16) is analogous to the proof of the formula (8.23) in 
Chapter IV. It is based on the skew-symmetry of d'^^ and dpq. 

Let's complete the inner vector field x of the surface by the other inner vector 
field y = 0(x). The vectors x and y form a pair of mutually perpendicular unit 
vectors in the tangents plane. For their components we have 



'^x''Xg = l, Xi = ^gikX^, y,=^g.^ky'', (6.17) 

9=1 k=l k=l 

2 2 2 

^VfcX^x, =0, % = ^a;p,xP, f=Y,^''xr (6.18) 

9=1 p=i i=i 

The first relationship (6.17) expresses the fact that |x| = 1, other two relationships 
(6.17) determine the covariant components Xi and Ui of x and y. The first 
relationship (6.18) is obtained by differentiating (6.17), the second and the third 
relationships (6.18) express the vectorial relationship y = S(x). 

Let's multiply (6.16) by Vkx'' Xj xP and then sum up over q, p, and j taking 
into account the relationships (6.17) and (6.18): 

Vfex' ^\Y.yi ^kx") f = y\ (6.19) 

. 9=1 / 



Using (6.19), i.e. substituting V^x^ = ZiyP and ^ jx'^ = Zj into (6.15), we see 
that the contribution of the second term in this formula is zero. Then, applying 
(6.16) to (6.15), for the increment Aip we derive 



ALp = - j j ^^a;'(V,VjX^ - VjVia;-'') Aet g du^ du^ . 

n '=1 J=i 

Now we apply the relationship (8.5) from Chapter IV to the field x. Moreover, we 
take into account the formulas (8.24) and (9.9) from Chapter IV: 



2 2 

// ^^{Kg.jx'x^) Act gdu^du^. 



i=l ] = ! 



Remember that the vector field x was chosen to be of the unit length from the 
very beginning. Therefore, upon summing up over the indices i and j we shall 
have only the Gaussian curvature under the integration: 



JJ K ^/dctgdu^du^. (6.20) 
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Now let's consider some surface on which a connected and simply connected 
domain Q outlined by a piecewise continuously differentiable contour 7 is given 

(see Fig. 6.1). In other words, we have 
a polygon with curvilinear sides on the 
surface. The Green's formula (5.1) is 
applicable to a a piecewise continuously 
differentiable contour, therefore, the for- 
mula (6.20) is valid in this case. The par- 
allel translation of the vector a along a 
piecewise continuously differentiable con- 
tour is performed step by step. The 
result of translating the vector a along 
a side of the curvilinear polygon 7 is 
used as the initial data for the equations 
of parallel translation on the succeeding 
side. Hence, ip(t) is a continuous func- 
tion, though its derivative can be discon- 
tinuous at the corners of the polygon. 

Let's introduce the natural paramct- 
rization t ~ s on the sides of the polygon 
7. Then we have the unit tangent vector r on them. The vector-function T(t) is a 
continuous function on the sides, except for the corners, where T(t) abruptly turns 
to the angles A'i/'i, A?/'2, ■ • ■ , Ai/in (see Fig. 6.1). Denote by tp{t) the angle between 
the vector T(t) and the vector a(t) being parallel translated along 7. We measure 
this angle from a to t taking for positive the counterclockwise direction. The 
finction iplt) is a continuously differentiable function on 7 except for the corners. 
At these points it has jump discontinuities with jumps at/ii, ATp2i ■ • ■ , ^i^n- 

Let's calculate the derivative of the function 'il'{t) out of its discontinuity points. 
Applying the considerations associated with the expansions (6.8) and (6.9) to the 
vector T{t), for such derivative we find: 

^ = (e(r)|VtT). 





^ — \ 
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Fig. 6.1 





(6.21) 



Then let's calculate the components of the vector VtT in the inner representation 
of the surface (i. e. in the basis of the frame vectors Ei and E2): 



i=l J=l 



(6.22) 



Keeping in mind that t = s is the natural parameter on the sides of the polygon 
7, we compare (6.22) with the formula (2.5) for the geodesic curvature and with 
the formula (2.4). As a result we get the equality 



Vtx — k ■ n 



curv norm 



n — fegeod ' Hi 



(6.23) 



But n inner IS a Unit vector in the tangent plane perpendicular to the vector r. 
The same is true for the vector 6(r) in the scalar product (6.21). Hence, the unit 
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vectors ninncr and Q{t) are collinear. Let's denote by e{t) the sign factor equal to 
the scalar product of these vectors: 

e= (e(x)|ni„ncr) = ±1. (6.24) 
Now from the formulas (6.23) and (6.24) we derive: 

'(/) = efcgood- (6.25) 

Let's find the increment of the function ■0(i) gained as a result of round trip along 
the whole contour. It is composed by two parts: the integral of (6.25) and the sum 
jumps at the corners of the polygon 7: 

/n 
e k good ds + ^ Alp,. (6.26) 

7 

The angle ALp is measured from x to a in the counterclockwise direction, while 
the angle ATp is measured from a to r in the same direction. Therefore, the sum 
Aif + Alp is the total increment of the angle between x and r. It is important 
to note that the initial value and the final value of the vector t upon round trip 
along the contour do coincide. The same is true for the vector x. Hence, the sum 
of increments ALp + Aip is an integer multiple of the angle 27r = 360°: 

ALp + Alp ^2Tir. (6.27) 

Practically, the value of the number r in the formula (6.27) is equal to unity. 
Let's prove this fact by means of the following considerations: we perform the 
continuous deformation of the surface on Fig. 6.1 flattening it to a plain, then 
we continuously deform the contour 7 to a circle. During such a continuous 
deformation the left hand side of the equality (6.27) changes continuously, while 
the right hand side can change only in discrete jumps. Therefore, under the above 
continuous deformation of the surface and the contour both sides of (6.27) do not 
change at all. On a circle the total angle of rotation of the unit tangent vector is 
calculated explicitly, it is equal to 27r. Hence, r = 1. We take into account this 
circumstance when substituting (6.20) and (6.26) into the formula (6.27): 

j j K y^det^dv}dv? + (j) e k g^od ds + ^ A^p, ^ 2Tr . (6.28) 



The formula (6.28) is the content of the following theorem which is known as the 
Gauss-Bonnet theorem. 

Theorem 6.2. The sum of the external angles of a curvilinear polygon on a 
surface is equal to 2tt minus two integrals: the area integral of the Gaussian cur- 
vature over the interior of the polygon and the integral of the geodesic curvature 
(taken with the sign factor e ) over its perimeter. 
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It is interesting to consider the case where the polygon is formed by geodesic 
hnes on a surface of the constant Gaussian curvature. The second integral in 
(6.28) then is equal to zero, while the first integral is easily calculated. For the 
sum of internal angles of a geodesic triangle in this case we derive 

ai + a2 + ce3 ^ TT + K S, 

where K S is the product of the Gaussian curvature of the surface and the area of 
the triangle. 

A philosophic remark. By measuring the sum of angles of some sufficiently 
big triangle we can decide whether our world is flat or it is equipped with the 
curvature. This is not a joke. The idea of a curved space became generally 
accepted in the modern notions on the structure of the world. 
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